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Canonical Light-Cone Commutators 
and Their Applications 
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I. Introduction 


I shall report here on recent work done by several colleagues! and me 
where we succeeded in exhibiting the canonical structure of the commu- 
tators of currents on the light cone, and applied it to several interesting 
physical problems. Two reasons may be advanced at the present time 
for studying this topic. 

(1) It has been appreciated already in the early days of current algebra, 
that the fixed mass sum rules, the most famous of which is due to Dashen, 
Fubini and Gell-Mann [1,2], which are usually derived with the help 
of the unreliable p— oo technique [3], are in fact equivalent to appropriate 
light-cone commutators [4]. 


* A. P. Sloan Fellow. 
1 This work was done in collaboration with J. M. Cornwall, D. Discus and V. Teplitz. 
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(2) Recently it has been shown that the dramatic MIT-SLAC electro- 
production experiments [5] measure directly the commutator of electro- 
magnetic currents on the light cone [6]. 

The following notation will be used. For the coordinate x” and for 
all other vector or tensor quantities, we define the + components by 


x = 27/2 (x + x’). (I-1) 


The remaining components will be denoted by i,i=1,2, or by the 

subscript L. The metric tensor now is g**=g~~ =g!?=g*!=0; 

gt =g  *t=-—g!!=-—g?*=1. Consequently x?=2x* x" — x?. 
Consider the commutator of two currents. 


Cao (x) = ER]. (1-2) 


In a given field theory, the usual canonical formalism permits one to 
evaluate the equal time commutator, C%}(x)|,o=0, Without solving the 
theory. I shall show how similarly one can compute the light cone 
commutator, C%}(x)|,+=09, Without unraveling the dynamics. The reason 
for the nomenclature can be given. Since C“}(x) vanishes for negative x”, 
and since x?=2x*x~ — x} =-— x? when x*=0, knowledge of C“}(x) 
at the point x*=0 is equivalent to the knowledge of C#}(x) on the 
light cone, x? = 0. The relation of the present results to other investiga- 
tions is as follows. Fritsch and Gell-Mann [7] independently and simul- 
taneously suggested that the most singular contribution to C4#}(x) near 

2 = () might be as in the free-field theory, where of course it is completely 
calculable. Subsequently Gross and Treiman [8] verified the Fritsch- 
Gell-Mann [7] conjecture in various interacting theories. For purposes 
of studying the deep-inelastic, “scaling” region [9], both our approach 
and the Fritsch-Gell-Mann [7] technique give identical results. However, 
as will be seen below, the determination of the fixed mass sum rules 
requires the full commutator at x* =0, not just the most singular part 
at x?=0. 


Il. Light Cone Current Commutators 


Since, as we have argued, light cone current commutators are of interest 
both for experimental and theoretical physics, we now show how one 
might calculate them [10]. Consider first a conserved, internal symmetry 
current Jë. The time independent charge is given by Q,= | d?xJ?(x), 
and the assumption that the current transforms in a known fashion in 
group space implies 


[Q., JEO] =i fasc (O) . (II-1) 
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It is not hard to show that for conserved currents an alternate formula 
for Q, may be given. 
Q= fP xidx Ji (x). (II-2) 


Of course Q, is independent of the unintegrated variable x*. From 
(II-1) and (II-2) it now follows that 


Ea 0), Jf OV] +0 = i fave JEO) 8? (x1) (x7) 


; (11-3) 
+0_ RR) + 6,575 (x). 


The terms in (II-3) which disappear upon integration over x” and x, 
are not unlike the Schwinger terms that are present in selected equal 
time commutators. Indeed by taking vacuum expectation values, one 
can show that they must be present, though of course the vacuum matrix 
element is not sensitive to their c or g number character. However in 
one important respect these structures are different from what has been 
encountered previously. S#,(x) and S} (x) are not, in general, local in x”, 
though they do possess this properly in x,. The reason for this is the 
following. The commutator must vanish for x?=2x*x~ — x? < 0. When 
x*=0, as it is in (II-3), x?=—x?, and by causality, the commutator 
must be local on x, (i.e. it has support only at x, = 0). However no con- 
straint on x” is imposed. Therefore we expect the objects on the right 
hand side of (II-3) to have an arbitrary dependence on x”. Such quantities 
are called bilocal operators. 

The bilocal operators are model dependent, and at the present time 
there is no a priori way of calculating them. However, as we shall de- 
monstrate repeatedly, they are measurable, physical quantities, and their 
form must be specified, if we wish to understand various and diverse 
physical phenomena. As will be demonstrated below, they are known 
to have non-zero connected matrix elements; hence they are q numbers. 
It will be shown presently that the poo technique of current algebra 
would predict their vanishing; hence we have the possibility of rectifying 
that unreliable method, if we can specify the bilocal operators. It is this 
task which we now embark upon. 


A. Light-Cone Quantization 


The reason it was possible to give a plausible model for equal time 
commutators without solving any of the dynamical equations of quantum 
field theory, was of course due to the fact that equal time commutators 
for fields form a boundary condition on the theory. They are assumed 
to be of a definite canonical form which can be given without solving 
the theory. Consequently for any function of the fields, like a current, 
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the appropriate commutators can be computed. It turns out that some- 
thing similar may be done for light cone commutators. Rather than 
quantizing the theory on a space like surface — which corresponds to 
the usual equal time quantization — one may equivalently quantize on a 
light-like surface. With the latter technique, light cone commutators 
emerge canonically [11]. 

In order to obtain familiarity with these ideas, let us begin by ex- 
amining a simple scalar theory given by the Lagrangian 


L£=50,089- Lee? tig. (11-4) 
In the usual way, this leads to the equation of motion 
ODgp+wWop=ip°. (11-5) 


However the quantum theory is not specified only by (II-5). One also 
postulates the following equal time commutation relations 


iLe(x), PO)],o-9=9, (II-6a) 
id PX), PO)] 0-9 = 8° (x), (II-6b) 
i[6o P(x), Cop (0)] 0-9 = 9 (II-6c) 


where the special role of 0,@ follows from the fact that it is the canonical 
momentum 6,9 =0£/50). This scheme is the quantum analogue of 
the classical (non-quantized) procedure of solving the partial differential 
(II-5) with the specification of the Cauchy initial value data: p(x) and 
ĉo p(x) on the space like surface x° =0. 

Let us now rewrite things in terms of the light cone variables. Of 
course the Lagrange function and the equation of motion are equivalent. 


L=08_98,9+39p dp -3Wp’+44p*, (11-7) 
28,6_9 +9; do +Wp=Ap°. (11-8) 


The new element comes in if we define the canonical momentum to be 
6 £/6 640 =0_@. Consequently we postulate that 


iL0_—(x), PO)].+=0 =24(x ) (x4). (II-9a) 

This may be integrated and gives the basic light cone commutator 
i[p X), 9PO)].- 20 =F E(x )67(x,) (H-9b) 
where e(x) = — e( — x) = 1 for x > 0. This procedure is the quantum analog 


of solving the classical theory by specifying the initial data on the light 
like surface x* =0. It can be shown that classically a unique solution 
is obtained by specifying only @. The knowledge of derivatives of ọ is 
unnecessary [12]. 
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Therefore (II-9b) is the fundamental light cone commutator in this 
theory, and light cone commutators of all other functions of p(x) may 
be now computed without solving the theory. There are peculiar features 
of this quantization scheme that should be commented upon. Note that 
the commutator between field and momentum (Il-9a) has an unexpected 
factor of 4. Furthermore according to (II-9b) fields do not commute 
with themselves, and it is also clear that canonical momenta are non- 
commuting. This difference from the state of affairs in equal time quanti- 
zation arises from the fact that the canonical momentum 7(x) is related 
by an equation of constraint to the canonical field: x(x) = 3- ọ(x). In 
the present scheme equations of motion involve + derivatives; the + co- 
ordinate plays the role of time. Consequently the canonical formalism 
is more complicated than at equal times, and careful investigation shows 
that the Ansatz (Il-9) is correct. A related phenomenon is that the light 
cone method has changed a second order differential equation in time, 
to a first order equation in +. Therefore the present scheme is somewhat 
analogous to the conventional quantization of Fermion theories, which 
are also first order. 

The question now arises whether or not the light cone theory is the 
same as the equal time theory, or whether it is different; i.e. is the S matrix 
the same for both schemes. This is an important question, for if the two 
theories are different, then we cannot assume the simultaneous validity 
of the conventional equal time commutators and the new light cone 
commutators. To answer this question with complete certainty would 
require solving the two theories, and comparing the results — a task 
impossible to carry out at present. Hence we must content ourselves 
with a partial answer which relies on various formal properties of the 
two theories. The following four facts are offered as evidence that the 
two methods of quantization result in the same physical theory. 

1. It can be shown that the free theories are the same [11]. 

2. If the light cone theory is developed to the point of computing 
S matrix elements, then one encounters instead of the usual Feynman- 
Dyson rules, the Weinberg “p-o” [13] rules, which are supposed to 
be equivalent [11]. 

3. In the light cone theory, Green’s functions are ordered along 
the + direction, rather than along the time direction. However causality 
requires the two to be the same, up to seagulls. To see this, consider 


T (x)= [A (x), B(0)] 0(x°) , (II-10a) 
T.(x)=[A(x), B(0)] 0(x*). (II-10b) 


The only place that (II-10a) and (II-10b) appear to be unequal is when 
x°>0 and x* <0, or x°<0 and x*>0. However in this region x? is 
spacelike, and the commutator function vanishes by causality [14]. 
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4. The vacuum expectation value of light cone commutators can be 
calculated in the usual framework. The results are consistent with the 
canonically postulated forms. For example, consider 


i<0| Le (x), o (0)] 10) = i| da?o(a?) A(xla?) (I-11) 
0 


where ọ(a°) is some spectral function, and A(x|a?) is the free-field com- 
mutator function. 


A (xla?) =(2n)~3 f d+ke(k®) ô(k?— ad)e x. (11-12) 


In the usual way, we have from the equal time commutators the result 
that 


fda? o(a?)=1. (11-13) 
0 


To calculate the light cone commutator, we observe that from (II-12) 
it follows that 
A(x|a7)|.+ =o = (—i/4) e(x7) 07 (x1). (I-14) 


Consequently, we have with the help of (II-13) 
i<0| Eo (x), PO] l0>lx+ =0=4 a(x”) 8’ (x1) (1-15) 


which reproduces the canonically postulated (II-9b), as far as vacuum 
expectation values are concerned. Other evaluations of vacuum expec- 
tation values always lead to the same result — the light cone method 
is consistent with and equivalent to the equal time method. 

Supported by the above four partial arguments, as well as by the 
fact that no conflict has been found between the two methods, we shall 
assume that the form of light cone commutators evaluated in the light 
cone quantization scheme is valid simultaneously with the equal time 
commutators. 

To give a form for the SU(3)x SU(3) current commutators, we turn 
to a model which has previously served as an inspiration for equal time 
commutators: the quark model. Since we wish to present results which 
are possibly true even in an interacting field theory, we shall assume that 
the quarks interact with a vector gluon — the end result does not depend 
on the nature of the gluon-quark interaction, as long as it is a non- 
derivative one. Thus we are led to the Lagrangian 


L =P) &,-gB,)y"— M]w 
ee ua En (Il-16) 
FH = "BY — 6° B". 
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We have set the gluon (B*) mass to zero in order that the computations 
be simple — the end results with which we deal are insensitive to this 
simplification. The quark field y has mass M. If we wish to introduce 
SU(3) symmetry breaking we would allow M to be a mass matrix, 
though for the present, we do not do this. The light cone quantization 
of this model may now be taken over from the published literature [11]. 
Various peculiar features of this scheme are the following. The 4 compo- 
nent Fermion field y cannot be viewed as a canonically independent 
field. Rather a projection of it y, =P, y, P} =4y7y" is independent; 
while the other projection y_=P_y, P_=%y*y” is a dependent field. 
The canonical dependence is expressed by the following equation of 
constraint, which may be shown to hold in this theory. 


(i6--gB*)y-=3[iö,-gB)y+M]y'y,.  (I-17a) 


Since B* is a massless vector meson field, there exists a gauge freedom, 
which may be exploited to set B* to zero [15]. Then (II-17a) is easily 
integrated, to give a formula, at fixed x*, for y_ in terms of the canonically 
independent fields B; and y+. 


p_(x)=(—i/4) | dée(x” — 2) 


: Í : R 4 (II-17b) 
. {Li 0,—gB;(x > é, x,)] y +M} ? p(x ’ Gs x) $ 
The canonical commutators are 
[B'(x), Bi(0)],.+ =0 = (i/4) g e(x7) ° (x1), (II-18) 
{p+ (x), pt O)}x.+ =0 = 27? P ô (x7) 6? (x,). (II-19) 


The operators p, and y¥ anticommute with themselves, and commute 
with the B’. Another commutator which we shall need is the one between 
pt and w_. We have from (II-17) and (II-19) 


{p_(x), yt (x) t= t 

= —(i/4 V/D e(x7 —x'7) {Li ð;—g B(x] 7 +M} 9* (x, — x). 
It is also possible to give other commutators, e.g. between y- and y*. 
These are extremely complicated, and fortunately we have no need of 


them. Note that the commutator (II-20) makes reference to the inter- 
action. 


(11-20) 


B. Derivation of Current Commutators 


The current commutators may now be calculated. For the present, 
considerations of internal symmetry are ignored, and the current is 
defined by 

t=pyy. (II-21a) 
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In terms of the fields y, and y_ the formula for the current is as follows 


Jt =2' pt p, (Il-21 b) 
J7 =2 P y*y_ (H-21c) 
Fa P yty yp, +2 yty yy. (11-21 d) 


The commutator of J* with itself is quite simple. Only the y, field is 
involved whose commutator, (II-19), contains no interaction terms. 


We find 
[J* (x), J*(y)].+ ay+ =9- (11-22) 


For the [J*,J~] commutator, we need the more complicated formula 
(II-20). A rather lengthy computation yields 


[J* (x), I~ Wer =y+ 

=(—i/2)e(x” —y~) B(x) Lid, +9 By] + M} P_w)5?(x1—y,) (1-23a) 
328" —y)ORLB(x) y' P_ py) yy} — hee. 

(The abbreviation h.c. stands for Hermitian conjugate which must be 

subtracted from the right hand side of (II-23a).) The commutator is 

not expressible in terms of the current itself. Like the commutator of 

the fields it appears to depend explicitly on the coupling to g B". However 


this dependence may be eliminated. By use of the equation (II-17a), 
(II-23a) may be reexpressed by 


IT e+ =y+ 
= {tex -y) P(x, -y) POY vy} 
+ 6F{ —fe(x” —y) 6?(x,-— yi) [ER yv) 
EB) yy wO] — he. 


We have used the relations y' P_=4'—e'/y;)s; 


(II-23 b) 


yayıyı'ıyy, H=-e'=1 for i=1,j=2. 
The remarkable feature of this final formula is its elegant simplicity. 
The light cone commutator of these two currents is expressible in terms 
of bilocal operators, which are straight-forward generalizations of the 
local currents. These bilocal operators are (x) y“p(y) and p(x) p“ys p(y). 
It is clear? that the operators enter with (x—y)?=0, since x* = y* and 


2 Our formulas for the bilocal currents are gauge dependent, since they were calculated 
in the special gauge B* =0. An explicitly gauge independent form is 


PO, x7, 0)»*yO) expig | dy” B+ (0, y7,0); 
0 


see Gross and Treiman [8]. 
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x,=y,. Consequently the non-locality is confined to the — direction. 
The bilocal terms do not contribute when an integration over x” and 
x, is performed. Thus the structure (II-23b) is of the type (II-3), which 
we arrived at by general considerations. 

Finally a completely similar argument gives 


[J* (x), FON 

= 0 {—ge(x” —y7) 67 (x,—y LP) y'vly) +e" BR) PO] 
+ OF {4 a(x” — yy) & (x. —y.) Lo B(x) »* vO) 
+e BR) yt ysy} — h.c. 


Again the light cone commutator is expressible in terms of bilocal 
generalizations of the vector and axial vector currents. 

Commutators of current components, not involving at least one + 
component, are considerably more complicated as they require the 
commutators between the y_ and yw* fields. We shall not concern 
ourselves with these here. 

There exists an alternate method for computing commutators, due 
to Schwinger [16]. Schwinger’s original arguments apply to ordinary 
quantization, but it is not hard to generalize them to light cone quanti- 
zation. This “action principle” procedure gives the following formula 
for the equal time commutator of the time component of a current, with 
the current itself. 


IR), J*(y)] x0 = 0 = — 10, 5 J*(x)/5.A,(y). (11-25) 


Here it is assumed that the current is conserved, and that it couples 
to an external field A”. At the end of the computation A” is set to zero. 
The variation in (Il-25) is performed with fixed canonical variables, 
ie. J* is first expressed in terms of canonical coordinates and momenta 
and the variation is performed only with respect to the residual depend- 
ence on A“. The analogous formula arising from light cone quantization 
is 


(11-24) 


[J* (x), MO] c+ =y = — 10, 6.I"(x)/5.A,(y). (11-26) 


In many models, J* has no dependence on external fields. (In our 
theory for example J* = 2 ył wy and y.. isa canonical variable which 
is held fixed in the explicit variation.) In that case (H-26) becomes 


II Wer == —10_ JT (X)/5A,(y) iS OAO). (27 


This formula is exactly of the general form (II-3). An abstract represen- 
tation for the bilocal operators is given by (II-27): —idJ~(x)/6A,(y) 
and —i ôJ'(x)/ A,(x), and it is seen that when an integration over x” 
and x, is performed, the bilocal terms cease to contribute. Moreover 
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we again recognize that these structures are light-cone generalizations 
of the ordinary Schwinger term. That object, it will be recalled, is given 
by öJ”(x)/6 A,(x), where the variation is performed at a common value 
of time [16]. 

A specific form for the bilocal operators is obtainable only by speci- 
fying a model, i.e. by specifying the dependence on A”. It is not hard 
to show that in the quark-gluon model this dependence is such that 
(11-22), (11-23) and (II-24) are true, regardless of the nature of the gluoan, 
provided it is not coupled derivatively. 


C. Non-Canonical Contributions 


Let us take the vacuum expectation value of the current commutator 
function 


<0 LI", PO] 10> = ("CO — 2") f d 220 (22) A(x). (1-28) 
0 


The representation is conserved, as a consequence of current conserva- 
tion, and a non-negative spectral function o(A”) has been introduced. 
Consider now the ++ components of (II-28) when xt =0. We have 
from (11-14), a nonvanishing result. 


<0] LJ * (x), J* (0) Od #=0 
= 00. (d920(22) Aw] c+ <0 (1-29) 
0 


= 8 ô [li/4) ex”) 52(x,)] Í d220(22) +0. 
0 


On the other hand the canonical light cone commutator (11-22) would 
lead us to expect zero for the right hand side of (TI-29). Of course, this 
is just the ancient problem that naive canonical commutators do not 
yield the conventional Schwinger term in Fermion theories. Recall that 
the Schwinger term is given by 


<O| EJ? (x), JECO [00-9 =i #8? (x) i dA?0(2?) (11-30) 
0 


which involves the same spectral integral as (II-29). Hence we learn that 
the light cone commutator (IJ-22) must be modified by a non-canonical 
contribution, just as the equal time commutator. We shall make the 
very important assumption that all these non-canonical additions are 
c numbers, and therefore given by their vacuum expectation value. This 
assumption is not true in perturbation theory, as will be seen in Section III 


Canonical Light-Cone Commutators and Their Applications 11 


below. Thus we set 
[J* (x), J Mer == (1/4) & & [e(x7 —y7) 6?(x,—-y,) S] 
S= [d220(22). (11-31) 
0 


It is instructive to examine this commutator in theories where the usual 
Schwinger term emerges canonically e.g. in scalar electro-dynamics. 
A straight forward calculation yields (IJ-31), except S becomes a bilocal 
operator S(x|y), with S(x|x) coinciding with the usual Schwinger term 
[17]. 

Similar modifications are present in other commutators. For the 
+ — components (II-28) gives 


<0| [J~ (x), J7 (0)] 10>, + =0 


i (11-32) 
= (i/8) [ 42202? — 6,4) a (4?) e(x7) 8°x). 
0 


Apart from the term not involving derivatives, which can be reproduced 
by the canonical commutators, we find in (I-32) a contribution propor- 
tional to Sô; é'67(x,). This second derivative structure has no analogue 
in the canonical result, hence the latter must be modified. 


[J* (x), JO] + =)+ 
= & {elx —y7)6?(x1—y) P(X) VO 
+g fiel — y7)67(x1— y LB) ' VY) — e P(X) 7750) 
—h.c. —(i/8) GF & {e(x~ — y~) 6?(x, — y1) S}. 
For the +i component (IL-28) gives 
<0| LJ* (x), JO] l0>x+ =0 = (i/4) 0_ Ge(x) 5°(x,)S+0 (I-34) 


while the canonical form gives zero. Hence this commutator is also 
modified 


[J* (x), TOJ- =yt 
= & {2887 y7) x y Pr wl) Pyp 
+ {elx =y) a ya) PE) 97 VY) 
+(x) ysp(y)]} — hc. + (i/4) & & {ex —y~)07(x,—y,)S}. 
In this discussion we have pretended that S is finite. If it is quadratically 
divergent, which is the scale invariant result, then in addition to the 
exhibited, quadratically divergent vacuum singularities, there are also 


well defined, non-infinite terms involving higher derivatives of 5(x_,). 
We shall not concern ourselves with these objects. 


(11-33) 


(11-35) 
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D. SU (3) Generalization 


The SU(3) generalization of the previous results is obtained easily. The 
current is given by 


V(x) = B(x) y3 Ava). (11-36) 
We shall also need the axial current 
Au) = iP (x) y"y53 4,0). (11-37) 
The internal symmetry matrices satisfy 
Aut = (i fare + dabe) Ae- (11-38) 


In deriving the current commutators, we allow for non-conservation 
of the currents by introducing the mass matrix M in the Lagrangian 
(II-16). The commutators are the following 


LV." (x), Vo" Wee =y+ fe ER =y) a y) 
—(i/4) 54, & [ex =y) P(x, ya) S], 


IR) Mer = + — i fare Vo (0) 5(x™ yT) y) 
= ~} {i fave + dano [Lex —y7) (x1 ya) Ke (Kly)] 
+3 elx —y7) 8y) (Vlxly) ie Aly 
+ gie — y) 8 (x ya) PO) y Y Aap) — h.c. 
310, {ex —y~) 6?(x,—y1) S}, (11-40) 
LV," (x), Vi s+ =y+ -ifane Vi(x) 8x" — y~) exay) 
= ~$ {i fave t dane} {Tex -y Ya y) (Villy) ie Aly] 
+ [a(x =y) y g V yti At ly 
+ sie =y) e xy) P) y* Ay) — he. 


(11-39) 


+30, tea —y) e (x y1) S}. (11-41) 
We have introduced the bilocal generalizations of the vector and axial 
vector current 

Vi(xly) = B(x) yp) (11-42) 
Ab(xly) = iP (x) y*y5 p(y). (II-43) 


These enter the commutation relations with (x— y} =0, x7 — y7 £0. 
The term involving A,, is an internal symmetry operator given by 
Aap = [M, A,] 4, Le. it probes the non conservation of the current VA 
It is not hard to see how the commutators with the axial current work 
out, but we do not present them here. In Section IV the model dependence 
of these formulas will be discussed. 
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E. BJL Theorem on the Light Cone 


We conclude this Section by presenting a method for computing light 
cone commutators directly from Green’s functions. The method is analo- 
gous to the usual BJL [18, 19] technique which allows one to determine 
an equal time commutator. Consider the x* ordered product of two 
operators 


T, (q) = | d*xe'**<a|T, O, (x) 0,(0) |B. (11-44) 
An integration by parts gives 
T, (a) = (i/q7) | dxe" 3, <a] T+ O, (x) 02(0) |B> 
=(i/q7) | d’x dx” e't ein al [O (x), 02(0)] IBdr+ =0 
+0 [(1/47)°]. (11-45) 


Not unexpectedly, the 1/9” term in a x* order product is the light-cone 
commutator, just as the 1/q° term in a time ordered product is the equal 
time commutator. The non-trivial aspect of the theorem emerges when 
we recall that a covariant Green’s function (T* product) is equal to 
the T, product, apart from seagulls. Consequently we have 


T*(q) = | d*xe'** <a] T* O, (x) 02(0) |B> 
q > polynomials +(i/q”) f dx dx ett = e”’a'* (11-46) 
-<al [O1 (x), 02(0)] 1Bdx+=0+ 0 


A more careful analysis indicates that q7 should become large away 
from the real axis. In the usual way, if the limit diverges, we interpret 
this as the statement that the commutator has divergent matrix elements. 

As an example, consider the free propagator of Boson fields. For 
large q7 it goes as i/2q* q7 and 1/2q* is indeed the Fournier transform 
of the canonical (free field) light cone commutator: 


—fie(x7)6?(x,). 


III. Applications to Deep Inelastic Processes 


The theory which we have developed will now be applied to the electro- 
production processes studied by the MIT-SLAC experiment. As is well 
known, that investigation gives an experimental determination of the 
quantity | d*xe'4*<p| [J*(x), J°(0)] |p>, where |p> is the spin-averaged 
nucleon target state, and J* is the electro-magnetic current. (In the next 
section, we generalize the present results to include spin and internal 
symmetry.) 
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A. Scaling Representation 


One of the remarkable experimental observations is that the invariant 
functions which determine the commutator defined above have a con- 
vergent limit in the “deep inelastic region”, q? and v =p -q both large, 
and w=-—gq?/2v fixed. We may inquire what must be the form of 
<p| [J*(x), J’(0)] [p> in position space, such that its Fourier transform 
exhibit the properly convergent form in the deep inelastic region. The 
desired form, which incorporates this regularity, as well as causality 
and current conservation is the following [20]: 


i<p| [J*(x), J*(0)] |p> 

= [I — 0] e(x°) L(x?) (1/82?) | do w7? -cosæxp Flo) 
+0(x?) f(x, x- p)] (III-1) 
+EP O —p- ep’ + ap") +g""(p - 87] 


o(<°)0(x2)[(1/8n2) | dolwxp)"!sinox-p-F,(w)+ fy(x?,x-p)], 
“1 


x? f,(x?, x: p) zag 0, (III-2a) 
a(x’, x: p) z 0. (III-2b) 


The non-trivial aspect of this representation is that the degree of singu- 
larity at x*=0 is limited, and the most singular terms at that point are 
parametrized by the deep inelastic structure functions of Bjorken [9], 
F,(o) and F,(@). The f;(x?, x- p) determine the subdominant light cone 
singularity, as is seen from (III-2). The degree of singularity at x?=0 
exhibited by (III-1), is no worse than in free-field theory. Also the singu- 
larity is such that the light cone commutator, i.e. the restriction to x* =0, 
exists. 


B. Equal Time Sum Rules 


The representation (III-1) is especially convenient for deriving the usual 
equal-time sum rules appropriate to deep-inelastic processes. Thus by 
going to x°=0 in the 0i components of (TII-1) we obtain the Schwinger- 
term sum rule [20]. 


pl IR), FO] Ip >1x0=0 


+1 (111-3) 
= (i/4n) &5° (x) | dw@-a@~?-F,(@). 
“4 


Similarly by differentiating the ij components with respect to x°, the 
Callan-Gross [21] sum rule emerges, at x°=0. We have after some 
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manipulation 


‚im (p) ~ f dèx <p] êo JO), FO] IPdlo=0 
m l l (ITI-4) 
= | do[F,(@)ö*- F, (œ) (&*- PN]. 


Clearly this procedure may be continued to obtain relations between 
+1 
| do-@”"F,(w)(i=L,2) and 
-1 
Jim (3)"""" fd? x pj LG" **J4(%), JO) IP>Iao-0- 


One disadvantage of these relations is that one needs to compute more 
and more obscure commutators, which depend on dynamics in a compli- 
cated way. Furthermore, it would be preferable to have one relation 
determining F,(@), rather than the infinite number of moment relations 
for f dw w?"F;,(w). The light-cone commutators provide these desirable 
results. 


C. Light Cone Sum-Rules 


Rather than restricting the representation (III-1) to equal times, x° = 0, 
we may just as well restrict it to the light cone, x* =0. It is easy to show 
that 


DIR, J (0)] IP e+ =0 


1 (II-5) 
=6_6_ forsone f do- w7: cosx p* F(a}. 


Comparing this with the model for the light-cone commutator, given 
in (II-31) or (IT-39), we have [22] 


+1 
(47)! | dw-w~? coswx p* F,(@)= <p|S|p>=0.  (II-6) 
=1 


The last equality follows from the fact that the Schwinger term is taken 
to be a c number; hence its connected matrix element vanishes. 
Eq. (ITI-6) contains the Schwinger term sum rule, (III-3); the Callan- 
Gross sum rule for F, (III-4); as well as all the further moment sum rules 
1 


for f dww?"F,(@). It is very important to note that (III-6) may be 
ab 
inverted to obtain the unintegrated result 


F,(@)=0. (III-7) 
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Thus Eq. (III-7) may be arrived at without use of positivity for F; (œ). 
This happens for all the deep inelastic current algebra sum rules: rela- 
tions which in that context involved an integration over œ, are now re- 
derived in an unintegrated form. 

Further results are obtained from the + — components. For simplicity 
we set F, to zero, as indicated by (TII-7). It then follows from (III-1) that 


DLR, J7 (0) [p> |x+ =o 


= —ie(x~) d°(x,)| f (x7 p*)+(p?/8 2) f dwcoswx” p* F (œ) (II1-8) 


—ie(x”) p ô ô (x) (1/87) f dosinwx”p*’F,(w)/o. 


Here f(a) parametrizes corrections to scaling, and will not concern us 
here; but see Chapter IV. 


f(a) =(M?/8n) f dw coswa F,(w)+ d/da[af,(0,«)]. — (III-9) 


‘The expression (III-8) may be compared with the operator formula 
(we ignore c number Schwinger terms) 


BEI Wes) 

=-3 & {e(x” — y~) ô 7x — y) VT (aLy)} (111-10) 
—4 GF {e(x7 - y7) 67(x, -y) LV'(xly) ie Ay) — h.c. 

In a Fermion model, from which (III-10) is abstracted, we have 


J*(x) = B(x) *Oyp(x);  V"(xly)= B(x) YO? yy), 
A*(x|y) = iiP(x) "7507 p(y), 


where Q is the charge matrix. (In the quark model a more specific formula 
may be obtained since Q is given by the SU (3) A matrices.) Consequently 
(III-6) and (III-10) imply [22] 


—(p'/27) j do sinwx” p* F(o)/o=<pliV'(x|0)+h.c. |p>x2-0 
-1 l (11-11) 
= (pl iP (x) Qp 0) + h.c. |p>|,220.- 


As always, the bilocal operator is evaluated on the light cone, x? = 0, 
with x* =0,x, =0. 

This then is the final, remarkably simple result, which summarizes 
all the moment sum rules. It is very gratifying that the bilocal operator, 
whose matrix elements determine F, (œ) has no dependence on the gluon 
field (in the gauge B* =0). As was noted before, the bilocal current is 


Canonical Light-Cone Commutators and Their Applications 17 


a simple generalization of the usual local current. Now we see that just 
as the matrix elements of the local current (the form factors) are relevant 
to elastic scattering, so analogously the matrix elements of the bilocal 
current determine the form factors appropriate to deep inelastic scattering. 
Of course we do not have at the present time an evaluation of 


<pl iP (x) Q?y"p(0)+ hc. |p>, 


just as at the present time one cannot evaluate the elastic form factors 
<p |B) Oy*w(0) |p>. The importance of the result (III-11) lies in the 
fact that it demonstrates the measurability of the matrix elements of 
bilocal operators. Thus these operators should be considered on the 
same footing as local currents in their relevance for practical physics. 

Eq. (III-11) may be considered as a starting point for the calculation 
of the deep inelastic structure functions, Rather than computing the 
cross-sections for all g? and v, and then passing to the deep inelastic 
limit, one needs to compute only the matrix elements in (III-11). 

It is clear that similar investigations can be carried out for the deep 
inelastic neutrino processes. All previous current algebraic sum rules 
which provided relations between moments of the appropriate structure 
functions are now replaced by relations between Fourier transforms of 
these functions; or equivalently between the functions themselves. 

We conclude this section with several further (unrelated) observations 
about the bilocal operators. If ip (x) y*p(0)|,2-, is expanded at small x” 
in powers of x”, one encounters local structures of the form 


py da arp. 


In particular, the Fermion part of the energy-momentum tensor occurs 
in this expansion — though the Boson part is missing [22]. For large 
x”, only small œ, i.e. the Regge region, contributes. Thus the bilocal 
operator at large x” parametrizes Regge poles. Finally we repeat once 
more the caveat that in perturbation theory the commutators possess 
anomalies and are not of the canonical form given here. One consequence 
of this is that F, +0, in perturbation theory, contrary to (III-7) [23]. 
Evidently we are ignoring here these perturbative q number anomalies. 


IV. Fixed-Mass Sum Rules 


As Prof. Furlan has explained in his lectures here [3], the ordinary 
current algebra is insufficient to derive fixed mass sum rules, like the 
one of Dashen, Fubini and Gell-Mann [1,2]. One possible technique 
which does yield fixed mass sum rules, is the p-+oo procedure. Unfortu- 
nately that method frequently fails because of illegitimate interchange 
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of limit and integral and many of the results are invalid even in free-field 
models. On the other hand it has been known for some time that the 
fixed mass sum rules are equivalent to appropriate light cone commuta- 
tors and that the p > technique is an attempt (not always successful) 
to determine the (then unknown) light cone commutator from the (then 
known) equal time commutator [4]. However since we now have a 
model for the light cone commutator we can dispense with the unreliable 
p>x method. We derive in this section corrected versions of various 
fixed mass sum rules [17]. 


A. Kinematic Preliminaries 


We consider the diagonal matrix element of the commutator of two 
vector currents V;", assumed to be conserved, between Fermion states, 
which are not spin averaged. 


Cap, q) = f d*xe'* Cpl CVE), VO] Ip). (IV-1) 


The Fermion state |p) is characterized by a spin (pseudo) vector 
s*=u(p)iy*ysu(p). This vector is orthogonal to p and has the form 
s” = (5°, 8),s°=p-n,S=mnit+pp-n/[E+m]. Here A is an arbitrary unit 
vector specifying the rest frame spin direction, A=<s). Note that 
p=- =m. 

Considerations of parity, time inversion invariance, Lorentz invariance 
and current conservation give the following expression for C#}(p, q). 


CHp, a) =C~9"" + 9"q"/9q71 We 
+ [p*p” —(v/q?) (p*q’ + p’g"(+g"”v?/g?] We" (IV-2) 
+ ie"? s qg WE? + ie p gga: s WR. 

The invariants are functions of q?,v. We shall decompose them into 


parts symmetric in ab, denoted by (ab); and antisymmetric in ab, denoted 
by [ab]: 


We = We iw, i=L,2,3,4. (IV-3) 

Crossing now implies that 
W,2)(q?, v) = — W(q?, —-v). i= L, 2,3, (IV-4a) 
WE?) = WEP, —v) (IV-4b) 


and the opposite symmetry obtains for the invariants which are anti- 
symmetric in a and b. Hermiticity insures that the invariants occurring 
in the right hand side of (IV-3) are real. 
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Since the current has dimension 3, and since our states are covariantly 
normalized, the following functions are dimensionless: W£?, y W2®, y We? 
and y? W$’. Hence these objects approach a limit as v>co with fixed 
—q’/2v=q@*. In the present context, this is not a hypothesis, but a 
consequence of our use of the light cone commutators. In the pre- 
asymptotic region, we shall frequently consider the above functions to 
depend on wand °. In that case W£? will be denoted by (— 1/20) F#°(c, 4”), 
and the others by F¥°(c, q?). In the scaling limit, the F%°(w, q?) become 
F#°(@), i= L, 2, 3,4. The quark model commutator structure which we 
employ implies that F??(w) vanishes, as will be seen below. We shall 
need the correction to scaling for this function. Hence we define 
lim q?F#°(a, q?) =G (œ). One can show that F°? and G“? vanish 


-q2> 


for |o| >1. 


B. Infinite Momentum Derivation of Sum Rules 


We summarize the p > co sum rules which emerge from the 0v components 
of ([V-1). Some of these results are well known, others presumably have 
been recorded in the literature. Eq. ([V-1) and the usual current algebra 


imply 
em) J dq? Calp, a) =i farch T. 


ol Va (0) Ip = p*r; . 


(IV-5) 


We have taken the Schwinger term to be a c number. Throughout 
subsection B we set p-q=0. A change of variable is now performed 
0 0 
q = vip". 


an | (aA CED) Sifa To A-9 


The above is exact, however it is not a fixed mass sum rule since 
q? = v7/p§ — q? varies with v. To obtain fixed mass sum rules, p, is set 
to œ in the usual fashion. Then by equating independent tensors in 
(IV-6b), the following non-trivial relations are obtained. 

From «=0 in (IV-6), the Dashen, Fubini, Gell-Mann [1, 2] result 
emerges. 


fdv WER, v) =n faele g <0. (IV-7a) 
o 


3 This is in complete analogy with the electroproduction results discussed above. 
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From «=i in (1V-6) the “good-bad” sum rules are a consequence 


fv w£1(q2, ») =0, (IV-7b) 
0 
{dvv mel, ») =0, (IV-7c) 
0 
fay Wag, ») =0, (IV-7d) 
0 
Tayweng?, ») =0, (IV-7e) 
o 
far, y=0 gs0. (V-79 


The sum of (IV-7b) and (IV-7c) is the Bég sum rule [24], which when 
evaluated at g?=0 gives a relation for anomalous magnetic moments. 

If the results ({V-7) are compared with the explicit calculations in 
the free quark model, it is then found that (IV-7b) and (IV-7f) are not 
valid, while (IV-7c), (TV-7d) and (IV-7e) are trivially satisfied, in the 
sense that both sides of the equality vanish. Numerical evaluation [25], 
which checks the validity of these in nature is consistent with (IV-7a), 
in that the Cabbibo-Radicatti relation [26] which is a consequence of 
(IV-7a), appears to be satisfied. However the Bég sum rule [24] is not 
verified experimentally. Finally ([V-7f) cannot be true since it can be 

1 


rewritten in terms of the scaling variables as | dw- œ~? . FS (w, q?)=0. 
1 = 
When q? is then set to 00, one gets J dw:+w~?- F¥”(@)=0 which is 


inconsistent with the MIT-SLAC Specht which show that 
FY” (w) + 0 [5]. 

In the free-field model all integrals converge, since the invariants 
are 6 functions. Hence the failures in that context cannot be attributed 
to divergences. In a Regge pole model the relation which diverge are 
(IV-7d) and (IV-7f). Note that the Bég sum rule converges in a Regge 
pole model [27]. 


C. Light Cone Derivation of Sum Rules 


To derive the fixed mass sum rules with light cone techniques, we set 
q* to zero in the +y components of (IV-1), and integrate over q`. 


(27) Jar Car (P; Da =o 
=fd x en A <p| [J dx” V3 (x), V0] IP>\x+=0 2 


(IV-8) 
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With this procedure, we obtain directly a fixed mass sum rule, since 
q’=—q? when q* =0; i.e. q? does not vary in the course of the inte- 
gration. However, the right hand side of (IV-8) involves a light cone 
commutator, rather than an equal time commutator. We assume that 


this object can be computed from our Ansatz, (11-39), (II-40), and (II-41). 
pl Lf dx” V(x), VEO] Ip |x + =o 
= [dx <p (V(x), VO] IPyle+ =0- 


The nature of this assumption shall be discussed further in the Con- 
clusion. Finally we change variables in (IV-8) by setting v=p* q> —p,-q,. 
This yields 


Qn)! | ap) Carp.) 


(IV-9) 


q*=0 
a =(v+ Pp. 'a.)ip* 


= fdèx dx emit t: <p EV (0), KOJ Pero YO 
We now use the model given previously to evaluate the commutator 
in (IV-10) for «= +, —. (No new information emerges from « =i). The 
Schwinger terms are ignored as they are taken to be c numbers. Also 
since the current is conserved, the commutators do not have symmetry 
breaking terms in them. Clearly the right hand side of (IV-10) will 
involve matrix elements of the bilocal operators. These can be given 
by the following form factor decomposition. Define 


Ve (xy) = 2 Vy) + 2 Vals), (IV-11a) 
FFD =a VEC) i VEO), (IV-11b) 
Ha (xly)= 2 Aa(xly)+ 3 AU), (IV-11c) 
AUX y) = 4; Ally) Fr ABIX). (IV-11d) 


The form factors are by definition 

pl Vi (x10) [p> = p“ Vi (x°, x - p) 
+ (x, xX: p), 
<p| F#1) lp) = p* V(x", x- p) 
FREI EKD), 
Cpl A (x]0) [p> = s" A3 (x7, xp) 

+ p"x-s Al(x?,x-p) (IV-12c) 

+x"x +s A3(x?,x-p), 
Cpl La(x|0) [p> = s“ Aa (x°, xp) 

+ p#x-sA2(x?,x-p) (IV-12d) 


+x"x +s A3(x?,x-p). 


(IV-12a) 


(IV-12b) 
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T inversion invariance eliminates a possible structure of the form 
eth x pgs, in (IV-12a) and (IV-12b). It is clear that (IV-12a) and (IV-12c) 
are even in x, while ([V-i2b) and (IV-12d) are odd in x. 

It is now straightforward to extract the sum rules by equating in- 
dependent tensors in ([V-10). The results are [17] 


{av WEG? v) = T fane Tes (IV-13a) 
ja WIR in fi j dat (0,0), (IV-13b) 
Ta WU, v) =} T fare f daA2(0, 0), (IV-130) 
0 
F avwa, ») =0, (IV-134) 
0 
dv W{(q2, ») =0, (IV-13e) 
0 


8 


J dv(v/—q?) W$ (q, v) =} ndase fda Vt, a), q?<0. (IV-13f) 
0 o 


The Dashen, Fubini, Gell-Mann sum rule [1,2] (TV-13a), as well as 
(IV-13d), and (TV-13e) are rederived; but the Bég sum rules [24] (IV-13b) 
and (IV-13c) as well as (TV-13f) are found to have corrections. This is 
extremely gratifying, since it is precisely these which fail both in free- 
field theory and in nature [25]. The corrections are expressed in terms 
of integrals over matrix elements of the bilocal operators. In the next 
subsection we shall show that these matrix elements are measurable, 
even when an integration is not performed. Here we merely demonstrate 
how (IV-13b), (IV-13c) and (IV-13f) can be exploited. 

Observe that the right hand sides of (IV-13b), (TV-13c) and (IV-13f) 
are independent of q?. Let us rewrite the left hand sides in terms of the 
scaling functions. 


I dv WI\(g?, v) = {doo 1 FIV g?), (IV-i4a) 
i dvy Wi*l(q?, v) = [dow Fo, ), (IV-14b) 
jav v/—q°) WiebXg?,v)= fiolo? 'Feb(o,g?), g?<0. (IV-14¢) 
According to the sum rules, the integrals are q? independent, and may 
be evaluated by letting —q*->0o. Thus the right hand sides of (IV-14) 


are expressible in terms of the deep inelastic structure functions F?"(a). 
Performing a similar analysis of the remaining sum rules in (IV-13) 
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gives finally the following results [17]. 


Farmer? „= fdoo™ Fa) = nl; (IV-15a) 

i dv Wi4"l(q?, v) = fawo Fo), (IV-15b) 

i dvv Wi#*l(q?, v) = daw! Flow), (IV-15c) 

j dv Wi2"l(q?, v) = duo" Fo), (IV-15d) 

Í dv Wh (q = |dw Fw), (IV-15e) 

(dvv(-g2)>! ing 2,9) = [doo iro o). (V-15 
0 


The sum rule (IV-15f) has already been derived by an entirely different 
method, by Cornwall, Corrigan and Norton [28]. Note that this relation 
diverges in a Regge pole model. Cornwall, Corrigan and Norton discuss 
a truncation technique, which possibly might give meaning to this 
divergent sum rule. 


D. Measuring Bilocal Operators 


It was seen in our discussion of deep inelastic scattering that matrix 
elements of the bilocal operators, with light-like separated arguments 
are completely measurable in terms of deep inelastic structure functions. 
Clearly the same is true here. This can be established in one of two 
ways. In complete analogy to the discussion of section IV, one may 
exhibit a position space representation for (22) 1 {d*gq e~'"*C#?}(p, q) 
which incorporates scaling as in (III-1). Alternatively one may work 
directly in momentum space and use the light-cone version of the BJL 
limit discussed in Section II. The results are as follows [17]. 


1 
dave V (0, 0) + faye V (0, a) = (27)! | daw el Fb), (IV-16a) 
SUL 
1 
day Ad (0, &) — fanc Ae (0, )=n' | dwe Fs (w), (IV-16b) 
1 
1 
dape% Az (0, &) — fane% A20, 0) = 77t | dwet Fiw), (IV-16c) 


=1 


1 
daped V20, 0) + furetV2(0,2)=¢s in f dom regeln). (IV-16d) 
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This is the generalization to spin and internal symmetry of (III-11). 
Also one finds [17] 

Fr?’(o) = 0 (IV-17) 
which is the generalization of (III-7). (Here F% = F{® + i Fl2?)), 


Additional results follow from (TV-16). Setting « to zero in (IV-16a) 
we have, since V(0, æ) is odd in «, 


1 
fare Vi(0, 0) = 107! (dom! Fl). (IV-18a) 
0 


The right hand side can be evaluated from the Dashen, Fubini, Gell- 
Mann [1,2] sum rule, evaluated in the scaling region, (IV-15a). Thus 
(IV-18a) becomes 


fave Ve, 0) = fare! - (IV-18b) 
Let us recall now the definitions of V!(0,0) and T,. 


P” farce Ve (0, 9) = fabe Cpl VF(0I0) |p> 
= fan.<pl V*(00) |p> ; 


p*r, = <pl VO) |p> - (IV-20) 


(IV-19) 


Thus (TV-18) shows that the proton matrix element of the bilocal operator 
fave V,'(x|y) reduces to that of the vector current as x y. This of course 
is obvious in the field theoretic model considered where V*(x|y) 
= P(x) 5 A,y"yp(y). However, if we generalize and abstract from the model, 
we see that (TV-18) places an important model independent constraint 
on the bilocal operator V*(x|y). 

We now assume that the same is true for the axial vector bilocal 
operator, as is of course the case in our model. We have from [17] 
(IV-16b) 


d25cA. (0, 0) s” = dane < p| A2(0|0) |p> 
= dqy-<p| APO) |P) = dap TAs" (IV-21) 


1 
=2s'n~! | dofo) 
o 
(41 (0,0) vanishes by symmetry.) The sum rule which has emerged, 


relating T^ to the spin odd, isospin even deep inelastic cross section is 
similar to a relation first obtained in the quark model by Bjorken [29, 30]. 
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V. Conclusion 


It is clear that the approach which we have here developed yields very 
informative and elegant results. A question which naturally arises is 
how generally valid are the conclusions of our investigation. It is this 
topic which we shall discuss in this Section. We shall not elaborate again 
on the canonical, formal nature of the present results, which invalidates 
them in perturbation theory. 


The bilocal operators which occur in the light cone commutators 
are certainly model dependent, even when they are interaction inde- 
pendent as in (II-39), (II-40) and (II-41). As we have already mentioned, 
scalar electrodynamics gives entirely different expressions, and the same 
is true of a Yang-Mills theory. These various commutator models may 
be said to differ among themselves in the space-time tensor structure. 
In addition to this there is the internal symmetry model dependence, 
which we have taken to be the SU (3) triplets of the quark model. 


It is plausible to suppose that nature favors the space-time tensor 
structure of the Fermion model, since the other models have non- 
vanishing F, or F, =F, [17, 20, 21]. In this connection it would be 
most interesting to explore what model independent information can be 
obtained about the [V,*, V] light cone commutator, if it is assumed 
that for «= + the form is as in (IJ-39)*. The internal symmetry structure 
which we have used does not, at the present time, possess any convincing 
experimental verification. Our results are true only in the triplet reali- 
zation of SU(3) where 4,2, = (i fanc + dan.) A.. TO generalize beyond this, 
one could for example set A,A, = i fabcAe + av) Where dias) is no longer 
dabee but contains all the symmetric parts of the 8&8 representation. 
Even further, one might generalize by setting 442, = i fan + dian where 
Jian is no longer the pure octet f,,.4, but rather is a general anti-sym- 
metric SU (3) matrix. 


With these internal symmetry generalization, the commutators would 
retain the same form as before, except that the S U (3) content of the bilocal 
operators would be more complicated. One would need to replace 
dayeVe'(xly) and d,,,Ae(xly) by Vew(xly) and Al(x|y), and similarly 
for the antisymmetric combinations. Consequently (TV-16) would possess 
left hand sides with more complicated SU(3) structures. Thus the 
measurement of the internal symmetry of the deep inelastic structure 
functions will provide important information concerning the validity 
of specific models. For example the inequality of the proton and neutron 


* This problem is analogous to the model-independent determination of the [V?, V$] 
equal-time commutator from the [V,°, VP] equal-time commutator [31]. 
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data [5] shows that the bilocal operators are not SU (2) singlets. In this 
context, we see that (IV-18) is an important, model independent constraint. 


The corrected deep inelastic sum rules ([V-15 b), (IV-15c) and (IV-15f), 
as well as the first equality in (IV-15a), (IV-15d) and (IV-15e) are quite 
model independent, since they are consequences merely of (IV-9), 
causality and scaling. The point is that according to (IV-8) and (IV-9) 
a fixed mass sum rule is given by an integral over x” and a Fourier 
transform with respect to q, of a light-cone commutator. This com- 
mutator is local in x, due to causality, i.e. it is composed of 6 functions 
of x, and derivatives of 6 functions. Consequently in Fourier space 
the integral over v of an invariant function W(—q?, v) must be a poly- 
nomial in q?. But the degree of the polynomial is fixed by scaling. Speci- 
fically, for example for (IV-15b) we can conclude quite generally from 
(IV-9) and locality 


ce) 1 
(dv WE- ai, v) = [daw ! Flo, — q?) 
0 0 
(V-1a) 
Fi 2)" clab] 
= q n t 
n=0 z 


However scaling requires that the limit g7->00 of (V-la) exists; hence 
we learn that 


foe) 1 
[dv wie 2, = CH! = [daw Flo) (V-1b) 
0 0 


which is (IV-15b). The second equality in (IV-15a), (IV-15d) and (IV-15e) 
is however a consequence of the specific model considered. 

Finally we come to the question of the validity of (IV-9), which is 
seen to involve an interchange of integration over x” with the limit x* —0. 
A careful investigation of this has been given elsewhere [17] and the 
result is the following. In the terminology of Adler and Dashen [32] 
the class 2 graphs, i.e. the graphs which have fixed mass singularities 
in the current lines, must satisfy the same superconvergence relation 
which is required for the validity of the poo technique. Thus as far as 
the class 2 graphs are concerned, the present considerations have nothing 
new to offer. On the other hand, the Z graphs, which were improperly 
treated by the p->co method, seem to be now accurately included, as 
is seen from the free-field example. 

Put in another way, the poo technique for (IV-7b) to (IV-7f 
requires the absence of all fixed poles; while the present method can 
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account for fixed poles with residues which are polynomials in q?; see 
(IV-15b) to (IV-15f). If there are non-polynomial residues, then the light- 
cone method also fails due to the invalidity of (IV-9)°. 
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Introduction 


Locality being one of the basic concepts of particle physics should not 
only be preserved in approximations but actually be used as a guiding 
principle for finding approximation schemes. It has been the objective 
of a number of investigations to develop local one particle approximations 
[1,2] and local approximations involving the two particle structure of 
Greens functions [2] which go beyond simple perturbation theory. 

It is the purpose of this talk to present an approach to the systematic 
construction of local approximations which can be completely charac- 
terized in terms of intermediate states. It could be considered a local 
form of a Tamm-Dancoff method [3]. The method has been developed 
in anumber of papers in collaboration with Gromes, Rothe, and Stech [4]. 


1. Local Saturation 


We consider a neutral scalar field A(x) of mass m which is relatively 
local to itself, i.e. 
[A(x),A(0)]=0 for x?<0. (1) 


The Haag expansion [5] of this field, i.e. the expansion of A(x) in terms 
of Wick-monomials of “in”-fields has the form 


AR = Au + nd | (i d'k) 3" r- >) 
n= =1 t 


i= i=1 


` h (ky, 2 k,) Anlkı) Bes Ain Ky): > 


30 H. D. Dahmen: 


where the Fourier transform A (k) is defined by 
Akx)=(2n)? f dtk e™™** A (k). (3) 


The properties of the coefficient functions h,(k,,..., k,) are the following: 

1. Lorentz invariance. 

2. Symmetry in the arguments. 

3. Restriction of the arguments to on shell values k? = m?. 

4. h,(k,,...,k,) =h*¥(—k,,..., —k,). 

5. h, are boundary values of analytic functions in each time compo- 
nent kP,(i=1,...,n). 

6. The Fourier transforms of the h, are retarded in all arguments. 

We note that these properties do not insure the locality of A(x) in 
the sense of Eq. (1). 

Because of the properties (1-6) we make the following ansatz 


fallky +k + ie)?) 

m? —(k, +k, + ie)? i 

fslikı +k tie)’, (ky +k3+ie), (ky +k +k; +ie)?) 
m —(k, +k, +k3+ ie}? 


hy(k,, k,) = 
(4) 
h3(k,, kz, k3) = 


etc., where the ie-presciption refers to the time component. The f, are 
real functions of the Lorentz invariant variables 


SalZis ees Zn) = fat (ZÑ, vs Zn) + (5) 


The ansatz (4) assures the weak convergence of the “interpolating” 
field A(x) to the “in”-field A;,(x) and to the “out”-field A,,,(x) defined 
through 


Aou(k)= Ain(K) + 27i elko) ö(m?—k?)j(k), (6) 

where the current j(x) is related to A(x) by 
j(x)=(C +m’) A(x). (7) 
Again the locality of A,.(x) relatively to itself is not guaranteed by 


Eq. (6). 

To insure the locality of A(x) as well as A,.(x) we have to study 
the infinite set of equations which derives from the insertion of the Haag 
expansion (2) into the commutator 


K(qg)=(2n)"* J d*xe'**[A(x), A (0)] . (8) 


This commutator can also be expanded into a sum of Wick-monomials 
of “in” fields 


K(q)=¢o(q) + 5 (n!)7} i I atk, ck... Ky) Ain(ky)»» Aig(Ky): - (9) 


i i= 
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Similar expansions hold for the retarded commutator 


K"*(q) = (27)! | d*x e'4* 0 (x0) [A (x), A(0)] (10) 
and the commutator of the “out” field 
Kr (q) = (27) + f d*x e!* [Aut (x), Aou (0)] (11) 


with the coefficients cH*(q,k,,...,k,) and ce(q,k,, --., Ky) respectively. 
Obviously there are relations connecting the functions f, with the co- 
efficients c*: 


n ’ 


Jar 1 n(m = q’) (m? —(q—-k)’) a (q, ky, tte kn)la-02=m2 (12) 

where B 
q=} k; and kurızg-k. (13) 

i=1 


The following relations can be derived with the help of Eq. (2): 

1. (q) = eldo) 5(m?—q?). 

2. Org, k)=0. 

3. co(q) is the Fourier transform of a local function for any arbitrary 
choice of the fys. 

All other coefficients c have to fulfil an infinite system of nonlinear 
integral equations relating the coefficients c, to the functions f,. This 
set of equations is most directly derived by 

(I) inserting the expansion (2) into (8), 

(II) expanding the product of the two Wick monomials into the 
products of contractions and single Wick monomials according to Wick’s 
theorem, 

(III) comparing the result of (I) and (II) with the expansion (9). 

We propose an approximation scheme for the solution of the condi- 
tion deriving from locality. The approximation consists in a truncation 
of the expansion in Wick’s theorem at a fixed number of contractions. 
The number of contractions taken into account defines then the order 
of our approximation. In the truncated equations every coefficient c, 
gets only a finite number of contributions involving a finite number of 
Ja only. Thus, to second order the coefficients c, are determined by only 
few of the functions: 


Co by h, 
cy by f, and fs, 
C2 by has fs and Sa, etc. 


This approximation scheme has the virtue that it bears a close 
connection to intermediate states and physical thresholds: The contri- 
bution of order | to a coefficient ¢ corresponds to intermediate states 
of particle number / (and a suitably chosen set of partly disconnected 
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intermediate states of higher particle number). It is different from zero 
only above the threshold at (Im)?. 

The problem is then to find functions c, and f,, such that the equations 
obtained in a certain order of approximation are fulfilled together with 
the Eq. (12) relating f,,, and cst, 

As we shall see in the following two sections the solutions to the 
first and second order approximations given there have the following 
properties. 

1. The coefficients c, are Fourier transforms of local functions. 

2. The TC P-conditions hold (up to the same order of approximation), 


1.€. 
(27)? fa= <0| j(O)| pis <--> Pas in) = KOI O Pi, ---, Pps Out>*. (14) 
Here the operator j(0) and the “out” operators, creating the “out” states, 


have to be understood as being defined in terms of the “in” operators 
as in Eqs. (2), (6) and (7). 


3. =0, nèi. (15) 


2. Saturation to First Order 


The first order local saturation is an approximation in terms of one- 
particle intermediate states. The equations for the coefficients c, and c3 


read 
&(g0) ôM? —q? 


q7>p-q 
an ST filp~4+i6))— ie ) (16) 


Pq, p) = 


&(qo) ö(m? —q?) 
m’ —(pyt+pr.—qtie)” 
x f3(*(p, + Po + ie)’, (py —q + iz), (pi +p -g+ ie)”) 
&(q—p,) 6(m? —(q—p,)’) 


cWq, p)= 


2 # )\2 2 + \2 x 
[m —(q+ ie) ] [m - (pı +p2—4+ie)] (17) 
x falla+ie)?) falpıtP2-g+ie)’) 
u ee +P2—-4 
17 Pi, P2>Pa2 
The solutions to these equations are [g+ =(qo+ié, q)] 
polynomial 

cre i = ; (18) 

re GP ama -p 

2 (q, P1P2) 


E polynomial | 1 
(m? —q*.) [m? — (pı +p2-4+)] (m? -— (q+ —Ppı) 


7 Hmm) . (19) 
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Thus, the solutions are uniquely determined up to polynomials. The 
equations for the higher coefficients c, are of similar structure and the 
solutions again are essentially of pole type. The requirements (1-3) 
stated at the end of the last section are fulfilled as far as they apply to 
this order of approximation. 

The first order approximation can easily be generalized to commuta- 
tors of currents involving a whole set of one particle intermediate states 
for their saturation. Also additional constraints like exact or partial 
conservation of currents can easily be taken into account. The equal 
time structure of current commutators can be studied in particular [6]. 

The approximations obtained from tree graphs in the effective 
Lagrangian approach [7] or via the hard pion technique [8] are of the 
same type as the first order approximation of this section. 


3. Saturation to Second Order 


In the second order local approximation (l= 2) the question of unitarity 

enters the game. No general solution has as yet been given, however, the 

solution for c,(q, k) exhibiting locality and elastic unitarity is known. 
The coefficient 


¢,(q, p) = (2) * f d*x e'4* <0] [A(x), A()] |p (20) 
is completely determined by the functions f) and f3:(qi =q+ié) 


cP (q, p) = cq, p) 


Lg.) 5-4), (P-q-)”) a 


+ fla) (m? -q2) [m?-(p-q_”]  \p>p en 


Here c{'(q, p), the first order approximation is given by Eq. (16); the 
s-wave scattering amplitude t is determined by 


të = Im | dk, [O(K9) O(g? —K2) +0(— k2) O(—q° +K) 


22 
Ser Bie ale 


Making use of the fact that c,(q, p) is the qo-discontinuity of c/°*(g, p) 
and of Eq. (12) relating f, and c{' we obtain the unitarity relation for fz 
as a special case of Eq. (21): 


Im f (47) = (qo) f2(44) t(42,m’) (Watson theorem). (23) 


34 H.D. Dahmen: 


The unitarity relation for t, can be obtained by the same manipula- 
tions from the second order approximation to c,(g, p1, P2). These relations 
turn out to be identical to the TCP requirements (14) for n=2, 3. 

The solution to the above equations is then given by (neglecting all 
possible polynomial dependences): 


_, fom’) HRA) + fas —p)*)—2 falm’)} 
(m? —q*.) [m —(q..—p)*] 


ert 


q, P) = (24) 


where f,(q) is determined by the expression [9] 


g-m? © 5(q’*) dq’? 


ae Tera 25) 


fa) = Sm?) exp 


Here ö(s) is the s-wave scattering phase shift as determined by tọ. 

The s-wave projection of the coefficient c3"(q, p1, p2) of the expansion 
of the “out” field vanishes. Thus, the four conditions stated at the end 
of Section 1 are fulfilled also in the approximation of second order. 
The solution obtained is no longer of pole type, both the functions fz 
and og, p) contain two-particle cut contributions with a threshold 
at 4m?. The function f} contributes only through its on shell s-wave 
phase ö(s), which is not restricted in this order of approximation. 

A general expression for c,(q, p) and c{"'(q, p) can be given in terms 
of an arbitrary choice of intermediate states by the following ansatz 


zn 9°) 
a°- -p-a 


eg, p) =r f dq? dQ? (26) 


with 


8 (qo) 8(—Qo) o(q"*, Q?) 


=(2n)""? f d*zd*y ee <0) A(z) jO AQ) lq +ay2=n2- = 
This representation gives a local expression for c{*(q, p) for every choice 
of intermediate states inserted to the left and to the right of the operator 
j(0). Even though, we know of no rigorous proof for the representation 
the Au of it looks very plausible: the pole residue of c¥*(q, p) at 
(p —q)* =m? must reproduce f,(q7) with the correct analytical properties, 
cC,(q, p) must exhibit the &(g,) and e(Po— 4o) structure and possess the 
reality properties contained in Eq. (21). Of course, this representation 
is consistent with our requirements and the approximations of first 
and second order, Eqs. (18) and (24), correspond to the following graphs: 
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q 


| 


a) b) C) 


Fig. 1. Graphical representation of a first order b and c second order contributions to the 
spectral function 


4. Conclusion 


The aim of presenting a local approximation scheme in terms of inter- 
mediate states has been reached in part. The first order approximation 
can be carried out, the solution has the form of pole structures for the 
matrix elements of the retarded commutator. The second order approxi- 
mation has been given for the vacuum one particle matrix element of 
the commutator. The problem of local approximations for matrix ele- 
ments with more particles is under investigation. 
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Introduction 


The concept of duality takes various forms. It has been applied to high- 
energy electroproduction in three ways; in historical order, these are: 

a) To build up at least part of the structure functions from an infinite 
sum of resonances, through a Veneziano-like amplitude. 

b) To relate the behaviour of the structure functions at extremely 
high values of v, q% to what is found at more moderate values 

c) In a parton model, where the partons are identified as quarks. 
Here duality provides a dynamical input, which gives information about 
the relative magnitudes of the contributions from quarks and anti- 
quarks. 

I shall talk about a) and c), leaving b) to be discussed by Dr. Rubin- 
stein. 

Although the first attempt [1] to explain the data for v Wf” through 
a Veneziano-like model preceded the other applications of duality, at 
present a realistic approach by means of Veneziano-like models has 
to be supplemented by the information contained in the parton approach 
[2]. The parton model determines the relative importance of I =4 and 
I =3 baryon trajectories, and also of isoscalar and isovector photons; 
at present we do not know how to do this using the Veneziano approach 
by itself. Also the parton model shows us how to incorporate non- 
resonance, pomeron-exchange contributions. 

One might be worried about identifying the partons as quarks, if 
quarks do not exist as physical particles. The excuse for doing so is 
that, when duality is incorporated, this is the only way at present known 
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in which we obtain results that are in fair agreement with all known 
data on the following experiments: 
(i) deep inelastic electroproduction on a proton [3], 

(ii) deep inelastic electroproduction on a neutron [3], 

(iii) deep inelastic neutrino scattering [4], 

(iv) electron-positron annihilation [5], 

(v) muon-pair production in nucleon-nucleon scattering [6]. 

There is of course a theoretical problem, if quarks do not exist. 
However, the mass of the quarks does not appear in the results of the 
calculations, and when the parton model is formulated [8] through a 
smoothed-out field theory, it does not seem necessary for the quark 
propagator to have poles. That is, it is possible that the quarks appear 
just as fields, and not as physical particles in asymptotic states. In more 
phenomenological language, perhaps the quarks are bound very tightly 
to their parent hadrons, and cannot escape. We certainly do not under- 
stand whether this is theoretically possible, but it is not obviously 
impossible. 

I have said that the results of the theory are in fair agreement with 
all known data. However, in most of the experiments the errors quoted 
are huge. It is worth saying that one should not expect the agreement 
to be perfect when the experimental errors get smaller. This is because 
duality is not an exact notion, being essentially dynamical; its predic- 
tions for purely hadronic reactions are in no more than fair agreement 
with experiment. 


Quark-Parton Model 


In the quark-parton model, the electromagnetic current and the weak- 
interaction current couple directly to quarks and antiquarks. Whether 
one treats the partons phenomenologically [9] or by the smoothed-out 
field theory [8], in the deep inelastic limit the dominant contribution 
to the virtual Compton amplitude arises from diagrams having the 
structure of Fig. 1. There the internal line is either a quark (Fig. la) or 
an antiquark (Fig. 1 b) and there is a sum of six such diagrams, correspon- 
ding to the photons coupling to each type of quark p,n,/ and to each 
type of antiquark p, 7, A. So if I write in explicitly the factors arising from 
the quark charges, then in an obious notation 


Pe = QP (FP + FP) + GP (P+ PP P+ F’). i 


There is a similar expression for v W£". If I use charge symmetry to 
relate the (p, n, A) contributions for the proton to the (n, p, A) contributions 
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for the neutron, this is 
FEN = @ (FH F)+Q) (FP + FP + FH FR). (2) 
In terms of quark fields, the weak interaction current is 


Tres =D y,(1—ys)n (3) 


where, for simplicity, I have set the Cabibbo angle equal to zero. From 
this, \ 

F}? =2(F" + PP) 

F3? = DPF + F"). 


In (4) the quark labels identify the quarks attached to the bubble in 
Fig. 1. If all the quark propagators are equal, the identity of the upper 
quark that propagates between the two current vertices does not affect 


(4) 


Fig. 1 


the numerical value of the contribution. Thus in the first term of F3?, 
an n quark is emitted by the proton, is struck by the positively-charged 
current and changes into a p, which is changed back into an n when 
it emits the second current. In the second term the sequence is (p, n, p). 
The factor 2 arises because each of these two sequences contributes two 
equal terms, one arising from y, coupling at each vertex, and the other 
from y„ys at each vertex. (The cross terms, with y, at one vertex and y,Ys 
at the other, contribute only to the structure functions F}, and not to F,). 

Now make a dynamical assumption, based on the two-component 
theory of duality that is familiar in the study of hadronic reactions [10]. 
Each amplitude is supposed to consist of two parts: 


R+P (5) 


where IP corresponds to pomeron exchange, while contributions to R 
are very small if they correspond to duality diagrams that contain a 
closed loop. Here this means that, so far as their contributions to the 
part R is concerned, the antiquark terms and the A-quark term in (1), 
(2) and (4) are negligible. These terms correspond to duality diagrams 
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of the structure of Fig.2a, which have a closed loop; the duality dia- 
grams (Fig. 2b) for the n and p quark terms do not. Notice that the 
duality diagrams are introduced just for book-keeping, and have no 
dynamical significance beyond this. I write 


FA=FR=FÄ=F4=0. (6) 


The duality diagrams say nothing about the diffractive part PP of 
an amplitude, and in fact quarks and antiquarks contribute equally 
to this. This is because the pomeron has even signature, so that its con- 
tribution to the bubble in Fig. 1a does not change when we apply crossing 
to that bubble to obtain Fig. 1b. Since the pomeron is also isoscalar, 
we thus have 


FR= FR = Fe = F2=D(w) say. 


x (7) 
Fi= Fé =D’(w) say. 
In phenomenological language, when we calculate the part R of an 
amplitude we may imagine that the nucleons consist just of three quarks, 
as in the most naive quark model. But for the diffractive part IP we must 
also take into account a sea of equal numbers of quarks and antiquarks. 
Notice that I am not going to make any detailed assumptions about 
the nature of this sea — it may be that some of the quarks and antiquarks 
in it are bound together to form “gluons”, but this will not matter. 
Since in the naive quark model the proton contains twice as many p 
quarks as n quarks, 


FR=2FR=R(o) say. (8) 
[More strictly, what can be shown is 


| (do/2) (FR-FR) =2,  { (de/20) (FR-FR) =1. (9) 
1 1 


These equations follow because it can be shown that applying the opera- 


tion f dw/2w to the sum of Fig. 1a and Fig. 1b results in the difference 
1 
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Fig. 3 


between the formfactor diagrams of Fig. 3, and the value of this difference 
at t=0 is known for each type of quark, from the isospin, charge and 
baryon number of the proton. We use (6) to remove the antiquark terms 
from (9), and then (8) follows if we make the additional assumption that 
the momentum distributions of the p and n quarks in the nucleons have 
the same shape. Verification of this last assumption comes from the 
scaling property of the elastic form factors of the nucleons, though this 
is not a proof because the momentum distribution of the quarks come 
into the elastic form factors in a slightly different way from that in which 
they contribute to Fig. 1a]. 
So now I have 


F5” (œ) =R (w) +} D(a) +5 D*(w) 

F3"(0)=3 R(w) +? D(w) +3 Do) (10) 
F3’(o)=2 R(w) +4 D(a) 

F3’(w) =4 R(w)+4D(o). 


I can find the function R (œ) in two ways. One is to try to calculate it; 
I will talk about this later. The other is to get it directly from experi- 
mental data: 


Fz? — Fy" =3 R(o). (11) 


The data [3] for the difference (F$? — F$”) is not yet very good (Fig. 4); 
for the present I shall pretend that it corresponds to the curve that I 
have labelled “m = 3” in the figure. Then the part R (œw) for F$” corresponds 
to the curve drawn in Fig. 5. This figure also shows the data [3], and 
I should like to interpret the difference between the data and the curve 
as the diffractive contribution (42 D (œ) + 3 D*(w)). We know quite gener- 
ally that the contribution to F, from the exchange of a reggeon of inter- 
cept &, is proportional to w*~* for large œw, and so for the pomeron 
it should be constant for large w. The diffractive contribution has been 
calculated [11] in a Veneziano-like model in which it is dual to cuts; 
it comes out the right shape, but its magnitude is a free parameter. If 
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the pomeron coupling obeys SU 3 symmetry, D+(œ) =D (w). With this 
extra assumption I can now calculate F3” and F3”. I can also [2] calculate 
F}? and F3” from R(w) and D(). I obtain, for the total cross-section 
for neutrinos scattering on propane, 


0.47 G ME/n per nucleon (12) 
and the experiments [4] at present give 
(0.52 +0.13)G?ME/r per nucleon. (13) 


There is similar agreement [2] with the details of the neutrino experiment, 
so far as they are known. The numerical predictions are not sensitive 
to the assumption of SU 3 symmetry for the pomeron coupling. 
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Muon-Pair Production 


Consider now the experiment [6] 


P+P->u* u +hadrons (14) 


where the total-energy variable s and the muon-pair mass variable |/ q’ 
are both large, such that 


t=q°/s (15) 


is not large. Drell and Yan [6] argued from the parton model that in 
the asymptotic limit the differential cross-section takes the form 


do/dyq? = (8202/3) (q?) >? y(t) (16) 


and they identified one contribution to p(t) as arising from the diagram 
in Fig.6. Here one proton emits a quark and the other an antiquark; 
these annihilate to produce a virtual photon, which in turn produces 


Fig. 6 


the muon pair. Drell and Yan calculated the contribution to y(t) from 
this term, and in the notation that I introduced in (1) obtain from it 


p(t) =) Q | dada, F*w,) F(z) ô(@ 0-1). (17) 


Here the summation is over all the types of quark and antiquark, and 
Q, is the charge of the quark a. The limits of integration over œ, and œz 
have to be chosen to correspond to the kinematical conditions of the 
experiment. From the foregoing analysis of electroproduction I know all 
six functions F“ and can plug them into (17). I obtain a different result 
from Drell and Yan because, as we have seen, the antiquark functions 
are rather smaller than the quark functions F?” and F”, since the former 
contribute just to the diffractive part of electroproduction, while Drell 
and Yan took all six F* to be equal. The result is the curve B drawn in 
Fig. 7; this figure also shows the data [6]. The curve stops at V7 =4 
because after that it is rather sensitive to the details of how the diffrac- 
tive contribution to electroproduction behaves near œ = 1; all I can say 
for sure is that it drops very rapidly indeed between the two vertical 
lines. 
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Fig. 8 


The analysis of the parton model in the smoothed-out field theory 
approach reveals [7] a contribution to p(t) from a second term, Fig. 8, 
in which a pomeron is exchanged after the protons have emitted the 
partons. (That is, it is associated with “wee” partons.) Strictly, this con- 
tribution is probably not a function of t alone, but also of log s. The 
details will depend on the unknown nature of the pomeron. However, 
as logs is a slowly-varying function, we can forget it for practical pur- 
poses. We cannot calculate the contribution from Fig. 8, because we 
do not know enough about the pomeron. But one can show [7] that 
for small t, the contribution takes the form* 

t Ẹla) (17) 

* Note added in proof : This statement is now withdrawn; see Nuclear Physics B36, 
642 (1972). 
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where p(t) is regular at t=0. Between V? =1 and Vr =2, varies 
between about 0.02 and 0.08 in the experiment, so this term probably 
dominates here. If I choose Ẹ(t) to be a constant to agree with the left- 
most experimental point, I obtain the curve A in Fig. 7. I expect the 
curve to continue to fall for larger values of Ver. Remembering the 
logarithmic scale, we see that the sum of the two curves A and B agress 
well with the data. 


Veneziano-like Models 


In the original formulation of the two-component theory of duality, 
it was hoped that the part R would be constructed mainly from reso- 
nances. Although there are indications now from hadron physics that 
cuts are likely to be very important too [10], in electroproduction the 
only known way to calculate R(w) is by means of resonances. 

The general problem of constructing Veneziano-like amplitudes for 
currents is very hard [12]. Particular obstacles are the requirements 
arising from current conservation, and from factorisation of pole residues 
on all trajectories, that is on both the parent and on every daughter. 
Fortunately, there are reasons to suppose that these requirements have 
little connection with the properties of the amplitudes in the scaling 
limit, so that a much more simple-minded construction [1] suffices. 

The reasons for this are, first, that W, arises from the double-helicity- 
flip virtual Compton amplitude and current conservation merely relates 
this amplitude to others, for which we do not need to construct models 
[13]. Second, the factorisation problem imposes conditions on the inte- 
grand of the representation of the amplitude in a particular part of the 
integration space; the scaling property arises from properties of the 
integrand in an altogether different part of the integration space, namely 
that part which gives the Fubini-Gell-Mann-Dashen fixed pole. 

There is, of course, no general argument that such a fixed pole should 
be associated with the Compton amplitude. However, we see from (10) 
that, because of the dynamical assumptions that I made, the function 
R(w) that arises in electroproduction also comes into the neutrino 
amplitudes. In fact, from (10), 


R(o) =3 (F2’(@) — F}”(@)). (18) 
The amplitude A(s, t, q7, q3) whose imaginary part gives (W7? — Wz”): 


Wz? — Wy? =2 M Im A(s, 0, q?, 4°) (19) 
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(the variables are defined in Fig.9) has /=1 in the t-channel, and is 
just the amplitude that has the fixed pole: 


A~ —2F,(t)/xs 
s large, t, q7, q3 fixed (20) 


where Fp is the Dirac isovector elastic form factor of the nucleon. 
Although the scaling limit is a different asymptotic limit, 


S70, q = =’, s/f =1-w (21) 


the presence of this fixed pole, together with the analyticity properties 
associated with duality, plays a crucial role in ensuring that there is 
scaling in the scaling limit. (In the parton model also, there is a close 
connection between the fixed pole and scaling). 

Because of crossing, I can write 


A(s, t)= Als, t)— Au, t) (22) 


and identify A(s, t) as the yP amplitude, or alternatively as its resonance 
part, since the diffractive part cancels between the two terms of (20). 
Then I attribute to A the structure 


Als, )=¥, Vals, t)+ È Vp-(u, t) 
B B’ 
+ ¥ Vep-(s, u). 
B, B‘ 


Here B, B’ denote the various baryon trajectories N,, A etc. that can 
be put in the s,u channels of the yP amplitude. The functions Vz(s, t) 
will be constructed to have the following properties: 

a) Poles at «(q?)=1,2,3,... and a(g3)=1,2,3,... where æ is the 
exchange-degenerate (9, œ, f, A,) trajectory. As the currents carry spin 
one, the poles at « =2,3,... correspond to the spin-one daughters of 
the leading trajectory resonances. 

b) Poles at «(t)=2,3,... (There are no poles at a(t)=0,1 in the 
double-flip amplitude). 


(23) 
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c) Poles in the s-channel corresponding to the trajectory B 


1 


d) Regge asymptotic behaviour t?%"% when tœ with s,q?, q2 
fixed 

e) Regge asymptotic behaviour s*~? when s>oo with t, 47, q3 
fixed, together with a term 


— ba(0)/ns (24) 


Vg (u, t) will have similar properties. The fixed poles of all the terms in 
(23), together with those of A(u, t), must add up so as to give the be- 
haviour (20): 


2 dl) - 2 bp (t) = Fo?) (25) 


(remember that u ~ —s in the limit (20)). 
Consider the function 


(du, du,dz/u,u22) (1+ 1/uy + 1/7 


og 


(1+ 1Y/zP9"?(1+2zX (u ,u,,2)P 9% Y(u,,u2, 2). = 
I require this integral to converge for negative values of the exponents, 
so when either u, or u, is large the function X is bounded and Y tends 
to zero. If also X and Y are neither zero nor infinite when the integration 
variables tend to zero or when zo, the integral has the same pole 
structure in the variables q?, q3, s and t as is desired for the amplitude 
V,(s, t, 7, q3). The s and t channels are dual to each other, their poles 
being associated with different ends of the same integration z, and their 
resonances produce Regge pole terms in the appropriate asymptotic 
limits. For example, when [14] B(s)>—- with t,g?,qg3 fixed, the 
dominant contribution to the integral arises from the region of integra- 
tion where zX is small. The Regge pole term arises from small z; to 


extract it make the change of variables z= —2z’/B(s), take the limit 
under the integral, so that for example 
(1+2zX (uy, ua, z)P 9 >e7?*unu2.0) (27) 


and perform the z’ integration: 


FR) (= Bi)" | (du, Aufn) (14 1u) 
y 28 
(1+ 1/u,) DT 1(X (uy, uz, 0)? Y(u,, uz, 0). (28) 


The functions X and Y must be chosen in such a way that the integral 
(26) has also the term (23) in its s-channel asymptotic behaviour. As 
the coefficient of s~* in this term depends on t but not on g? or q2, 
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the term must arise from the part 


Uj, U2 —> 00 (29) 
z finite 


of the integration region. That is, the function X must tend to zero in 
this region. A choice of the function X that has this property is 


Riis ee a. (30) 
i aa č+u +u) ` 


There will not be an attempt here to establish the essential uniqueness 
of this choice of X, as ultimately the choice depends on what is desired 
from the spectrum of daughter particles; the analysis of this is altogether 
beyond the scope of the present work. However, it will be remarked 
that the choice of X is very severely constrained by several other re- 
quirements, in addition to those already mentioned on its behaviour 
at the ends of the range of integration. Among these is a requirement 
that there be no further structure in the s or t channel asymptotic be- 
haviour (so that, for example, X does not tend to zero when u, alone 
is large) and that the function goes to zero exponentially when s— oo 
at fixed u=(p;— q2). Another natural requirement is that X (0, 0, z) 
= Y(0,0,z)=1, so that the amplitude for the purely hadronic process 


0 +p>g +p 


extracted from the residue of the twofold pole at x(q?)=«(g2)=1 be 
an ordinary Veneziano beta function. In order that the residues of the 
s channel poles have the appropriate factorization properties, it is also 
required that X (u,,u,, ©) and Y(u,, u,, 00) reduce to suitably factorized 
forms in u, and u,. The details of this are unimportant here, but it will 
be remarked that the form factors obtained by factorizing the s channel 
pole residues can be made to vanish as rapidly as one pleases for large 
current “mass”. An explicit example was given in Ref. [1]. 

I still have to make the correct function ¢,(t)/x appear as the coef- 
fizient of (—s)~' in (23). Make use of the known analyticity properties 
of dz(t) to write it as the transform: 


650) = | (dz/2 (1+ 1/2- bal). 31) 


(By making the change of variable 1 + 1/z = e* one can see that this is es- 
sentially a Laplace transform). Choose the function Y such that in the 
region (29) 

Y ~ (a/n) (1+ 2) (u, +42)" dsl). (32) 
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Then on making the changes of variable 
u, = — B(s) ui 


(33) 
u, = ~ B(s) u3, 
taking the limit under the integral and performing the u}, uw, integra- 
tions, we find the result (23). 

Consider now the Bjorken scaling limit (21). In this limit the factors 
in (26) that involve q? and g3 reduce to finite exponentials when one 
makes the changes of variable (33) and takes the limit under the integral 
[14]; that is the dominant contribution again arises from the region of 
integration (29). Using (32), we can perform the u; and u, integrations 
and find the asymptotic form 


— Pplt, w)/s (34) 
where 


a tg, 41/207! ba(2) 
p(t, œo) =n Ad z+(1-o)"! 


(35) 


If #,(t) is real for a(t) <1, its transform $,(z) is real. Hence $,(0, œ) 
is real for œ <1. If it is now continued to >| it acquires an imaginary 
part due to the vanishing of the denominator in the integrand: 


Im $,(0, ©) = — 0(@ - 1) 016 ,((w—1)7). (36) 


The existence of this imaginary part implies that @,(0, œ) has a branch 
point at œw =1, even though (26) has been constructed to have only poles, 
and no branch points. The explanation of this is familiar: the apparent 
cut is synthesised by the infinite numbers of poles in the variables s 
and q’. The absence of these poles in fact prevents the asymptotic be- 
haviour (34) from being strictly valid in a wedge containing the positive 
œ real axis, But it is supposed that it would be valid on the real axis 
in a unitarised theory, where the trajectory functions have cuts and the 
poles are displaced from the real axis. 

Because of (25), when I add together the contributions from all the 
terms in A(s, t) I obtain in the Bjorken limit (21) 


A(s, t, g’, q’) = 66, w)/s + Plt, —w)/u (37) 


where #(t, œ) is given by an integral like (35), but with ¢,(z) replaced 
by Fp(z) defined by 


Fy(t)= (dj) (141/21 F (2). (38) 
ie) 
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I have assumed that the terms V,,-(s, u) do not contribute in this limit; 
certainly it is hard to constructa Veneziano-like model in which they do. 
Because of (36) the first term of (37) has an imaginary part in the region 
&>1 of interest, with the final result 


R(o)= 09-1)! F,(o-1)"'). (39) 


The curves in Fig. 4 and 5 correspond to the choice 


r(i-a(t))  T(1+m-—a0) 


Ft) = 0 
DO= Firma) F(a) (30) 
which, for large t, has the behaviour F,(t)-t"". This gives 
T( — (0 
R(o)= ( +m ot ( )) ag — tt, (41) 


T(1-&(0)) T(m) 
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1. Introduction and Outline 


The large number of speakers who will deal with different aspects of 
highly inelastic lepton scattering at this summer school is evidence for 
the importance currently attributed to this subject. As the first, I will 
begin on an elementary level by discussing the kinematics in some detail 
and making a few preliminary remarks which should provide some back- 
ground for the following lectures for those who are not familiar with 
the subject. 

Next I will turn to the parton model. I will outline the usual “deriva- 
tion” of the model but the philosophy adopted subsequently will be 
that the parton model provides a useful heuristic device from which 
results can be abstracted which might be true more generally. In fact I 


* Work supported by the U.S. Atomic Energy Commission. 


52 C. H. L. Smith: 


will then show that these results can be rederived formally using infinite 
momentum commutators (although superficially different this part of 
my lectures will overlap with those of Gross and Jackiw). However, it 
should be noted that this formal derivation of the parton model is not 
necessarily superior to the traditional one — the difficulties are simply 
buried in formal (possibly meaningless) expressions. 

In the last section I give references to papers about processes which I 
shall not discuss and make a remark about the use of local duality in 
inelastic lepton scattering. 


2. Preliminary Remarks 
a) Electroproduction 


Kinematics 


We shall consider the process in which an electron scatters from a nucleon 
and only the final electron is observed. To lowest order in «, this is 
represented by the Feynman diagram: 


where k,k’,q,p and pp are four momenta, E and E’ are lab. energies 
and @ is the lab. scattering angle. 
With states normalized covariantly, 


PIP’) =2E(2r)’ö°(p— p') 
the differential cross section is given by 
do =(2M2E) +M (2r) > (2E) +d?’ k (27) ôt(k+p—-k' —~pr) (1) 


according to the usual rules. 
The invariant matrix element is 


M = e?ü(k‘, s’) y„u(k, s) (i/q*) <prl J,(0) p> 


where zu=2m, and we use the metric q? =q% — q°. On squaring and 
averaging (summing) over initial (final) spins we obtain: 


IM? (Qn)*ö*(k+p-K-pr)=(4re*/g*)m,,W** 2) 


Parton Models of Inelastic Lepton Scattering 53 
where q=k-k', 


m,.= >. Vik, s) y,u(k, STk, s’) pulk, s) 
=} Trk+ m) y,(k +m.) Y, 
= 2(k,k,+k,k,—g,,k- k’) + O(m?) 
W, (F)=4¥° <p| J; (0) F> <F| J,(0) Ip) (22)°04(q-+ p — pr) 


and È denotes a spin average over the initial states, which will be under- 
stood henceforth. 
If we sum over all the final states F, we may write: 


aM (F) 


(gu, 40/4?) Wi + (Py 4.4: Pla?) (p,— 4,4: P/Q”) W (4) 
W=Wiv,g?), v=q:p=M(E-E). 


(3) 


This is the most general tensor which we can construct out of the available 
vectors which respects current conservation (q’W,,=W,,q’=0) and 
parity conservation. Combining Eqs. (1)-(4) and ‘putting m,=0 we 
obtain: 


d?a/d|q?| dv= (n/M EE’) d*o/dQdE’ 


= (4n0?/q* M?) (E' JE) (W, cos? 0/2 + 2 W, sin? 0/2). ©) 


We shall collect some useful properties of the W which describe 
electroproduction below when we consider the structure functions which 
describe neutrino production. 


Scale Invariance 


We have defined W,,, (Eqs. (3) and (4)) so that it is dimensionless. Consider 
a theory in which all masses and coupling constants with dimensions 
are zero, so that no natural scale is defined and the theory is invariant 
under scale transformations x— Àx and p>p/A. If we consider the form 


of W,,, we see that in such a theory we would have to make the replace- 
ments: 
W, >F, (œ) 
W/M? > F,()/v (6) 
(w= 1/x=2v/((—q’)) 
in Eq. (4). 


Bjorken proposed [1] that Eq. (6) might hold in the real world in 
the limit v> oo with œw fixed, in which case masses might be irrelevant. 
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In fact it seems that it holds to quite a good approximation at rather 
modest values of v and q?. The relevant experimental results are discussed 
by Drees [2]. Here we only remark that in perturbation theory 


W/M’ sce Falo) + f(o)log(g*/M?). (7) 


The second term is not excluded by the SLAC data, although it turns 
out that f/F <1/10, if we fit the SLAC-MIT data [3] in the range 
1 Gev? <|q?| <10 Gev? with this form for œ = 4. 

Henceforth we assume that Bjorken scaling holds, i.e. that f (œ) =0. 
Note that 

a) If this is not the case, the rest of these lectures are empty. 

b) If it is the case, arguments based on perturbation theory may be 
irrelevant in the scaling region. 


b) Neutrino Production 


The differential cross section for neutrino production can be calculated 


as above with 
el + G//2 
Was = LI <P Je OMY FIIO IP Cta) (8) 
miy = Trey, (1 Fyk +m) F Ys). 
Quite generally W,, may be written: 

Ws" =g, Wir? + M? p,p Wz"? — (2M?) Nie, yagp*q? W3"* 
+M ?q,q, WE" +2 MP)" q,p,+q,Pp) WS" (9) 
+(2M*)~*i(q.Py—4yP,) We. 


Whatever the form of the lepton current, m,, may be expanded in a 
similar way with p,—k,. Hence m,,W*” is a second order polynomial 
in k- p= ME} and: 
E3(d?a/dq? dv) ~ m, W® ~ A(v, q?) + By, q°) E, + C(v, 4°) Ey. (10) 
This result [4] (sometimes referred to as a consequence of “locality”) 
depends only on the assumption of a first order vector interaction: 
u 


y 
—— Spin] + SpinO 
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It therefore allows a test of this assumption. From our point of view it 
is important because it implies that only three combinations of structure 
functions can be obtained without polarization measurements. 

In the conventional case W, , 3, are effectively determined, since the 
other structure functions are multiplied by q, and q,l*~O(m,). Ex- 
plicitly, Eqs. (8) and (10) give: 

d?a”*/d|q?| dv=(G?/2 M n) (E'/E) [Wz cos?0/2+2 W;*sin?@/2 (11) 
+W3M~ "(E+ E’) sin? 6/2] + O (m2) W,, 5. 

In principle W, 5, can be determined by measuring the muon’s 
polarization (the explicit expressions are given, e.g., in Ref. [5]) but this 
is hard, since it is necessarily 100% longitudinally polarized in the ap- 
proximation m, =0 in the V — A theory. 

For some purposes it is convenient to define the total cross sections 


for the fictious process of current — nucleon scattering. First we define 
the normalized, conserved polarization vectors: 


= (— q’) "7 (q5, 0, 0, qo) 
e® = 2-00, 1,1, 0) 
ek = 2-176, i, 1, 0) 
[4 = (qo, 0,0, 93)]. 


M 


(12) 


Contracting efe, with W”” and putting in the appropriate factors we 
obtain: 


Er 11/,27a44. 27472 
ORL v+rgnR (Witz /M*— q°/M“ W: u 


= (W,(1— v’/M?q’) — W,) 


su 


where, according to an arbitrary convention, we have used the flux 
factor for a zero mass (q7=0) current with the same s= (p + q)’. 

According to Regge theory, or our knowledge of real cross sections, 
we might expect that 


v> © a—1 


O Pixar Y > 
where a(t) is the leading Regge trajectory which can be exchanged in 
elastic current-nucleon scattering. Hence 

Ww, pied re 

W, > ven (14) 


03-1 
W; > yv ; 
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Presumably «, =a, =1 since the Pomeron can be exchanged. How- 
ever, since W, is the V— A interference term (it multiplies the only 
pseudotensor in Eq. (9)) it requires odd G parity exchange in the t channel, 
if only first class currents are present. Therefore we expect a3 ~ 1/2. 

Since Cg z, s Z0, according to their definition, we see that: 


W,(1—v?/M*q’) 2 WZM Yv -Mg w0. (15) 


In fact, the complete exploitation of the fact that W,, is a positive semi- 
definite Hermitean matrix is somewhat more complicated. The necessary 
and sufficient condition that the diagonal elements remain positive 
semidefinite under all unitary transformations is that all the subdeter- 
minants of W,, be positive semidefinite. This yields two more nonlinear 
inequalities involving W,, W; and W, [6]. 

Some other useful properties which follow immediately from the 
definition of W,, are: 

1. W,,=W- W, real 

2. If T is conserved, W,,(Pa; qp) = W*"(p%, q’) 

Wi... 57 real 
W; — imaginary . `. W, =0. 

(The only effect of W; is to give the final muon a polarization out of the 
reaction plane. It is easy to see directly that this transverse polarization 
is forbidden by T invariance in the approximation «= e?/4r =0 - see, 
e.g. Ref. [5]) _ 

3. The W7 have no kinematical singularities or zeros. However, 
since W,, is finite as q’>0 the W, satisfy 


Wy =" 0(q’) 
(W? + v? W3/M*q’)~ O(q’). (16) 
4. We can define different structure functions by summing separately 
over strange and non-strange final states and put 
W7 > W cos? Oo + o}? sin? be 
REN ie (17) 
FP — FY’ cos* Oo + f” sin" Oe 
where 6. is the Cabbibo angle. The assumption of the charge symmetry 
condition for the As =0 weak current: 
Jr = verte Ih Je eri 
gives: 
FL B) Jf p> =<FU, =B) Jin). 


Parton Models of Inelastic Lepton Scattering 57 


On summing over all states F the isospin reflection is irrelevant and 


we obtain si = er ig 
wer Wi", WP=W, 


a (18) 
FPP — F}”, FPP = FE" 7 


war) f d*x <p JF (x) JEO) DIET (a) 
=(4n)~* J d4x <p! [JF (x), JO) p> e* (b). 
It is easy to check that the additional term is the second line makes 


no contribution in the physical region v>0. Eq. (19b) defines the W, 
for v<0 according to the crossing relations: 


(19) 


W (v, g) TS We(- V, q’) 


; 20) 
W? (v, q) ~ Ww’(-», q ) > 


which are derived by examining W5. 
The advantage of this definition of W,, for v<0 is that the optical 
theorem holds for all v: 


W,,=(22)"'Im Ty, 


K . 21 
Tp =i f d* 0 (xo) <p LIF (x), OJ] p) =+ Polynomial inv. 9) 


T,,y is, of course, the amplitude for forward current scattering so that 
Eq. (21) is the usual relation between forward scattering and the total 


cross section: ; 
feos Nea ee 
F 


Thy = — Guy +M Pa ht (22) 


We may write 


in correspondence with Eq. (9). The T;(v,q?) are analytic functions in 
the cut v plane which satisfy the dispersion relations: 


+2 Wa, g?) dy’ 
ET J 

2 te Was, sV, Oy 4 = 
73,50, 9)=2 | Tr, dv 


assuming Regge asymptotic behaviour (Eq. (14)). (The equations for T, 4 
appear superficially divergent. However, using the crossing properties 
(Eqs. (20)) and the fact that W” — W7 does not receive contributions from 
the Pomeron, we see that in fact they are convergent.) 
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Scale Invariance 


In neutrino production, Bjorken scaling obviously takes the form 


wW, nr Fi (o) 


vW/M? => Fœ) (i1). 


Present data do not allow these predictions to be tested directly. 
However, there are two direct consequences of scaling, which can be 
confronted with the data already. The first is the behaviour of the total 
neutrino cross section. In a theory without any parameters with dimension 
the only possibility is (s = (p + 4)}°) 


o~G’s 


since G? ~ M~‘*. The second is for the average Q? at fixed s which must 
satisfy: 
<Q?) ~s. 


Both these predictions are compatible with the data [6, 7], although 
the values of s involved are rather small. 


3. Parton Model 


a) Intuitive Derivation 


The parton model consists essentially of applying the impulse approxi- 
mation to inelastic lepton scattering, the partons being Feynman’s name 
for the constituents of the nucleon [8, 9]. The impulse approximation 
is supposed to apply when 


A E interaction >A E binding 
or, in coordinate space, 
Tinteraction < TL 


where 1, is the lifetime of the virtual states of the nucleon. In this case 
it is argued that the constituents may be treated as free particles during 
the interaction. The necessary condition is achieved by considering a 
frame in which the proton’s momentum is very large (|pl>o0) so that 
T; is lengthened by time dilation while 1,,,~ 1/90 & 1z. 

A further condition is that the scattering from the constituents is 
incoherent. This is achieved by assuming that the partons momenta 
transverse to the nucleon’s momentum is sharply cut off (which is not 
unreasonable since the debris produced when the nucleon is broken 
up by collision with another hadron satisfies this condition) and con- 
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sidering conditions in which {gq | is very large compared to the cut off. 
In this case diagrams in which the current is absorbed by different partons 
cannot interfere (final state interactions being ignored since they turn 
one set of complete states into another — as we shall see in detail below). 

To see that these conditions may be satisfied in Bjorken’s limit, and 
to put the discussion above on a more quantitative basis, we consider 
briefly the arguments of Drell, Levy and Yan [10], who considered a 
field theory of pions and nucleons endowed ab initio with a cut off in 
transverse momentum (we refer to their papers for complete details). 

They use the interaction representation in which the complete Heisen- 
berg current operator J,(x) is related to the free field current operator 


ee: 3,09 =U") (2) UO 
where : 
U(t)=Texp—i | H,()dr'. 
Therefore: 2; 
w.=+), L <Upl j,(0) U(0) |F> <F| U” \(0) j,(0) [U p) (27)? (24) 
- “(q+ p— pr) 
|Up> = U(0) |p» 


H 
=VZ, {le +X om zn 


‚ m> <m| Hıln) <n| Hılp> 
nn") 


where & indicates summation over all states except |p and we recall 
that in this “old fashioned perturbation theory” the virtual particles are 
on mass shell and momentum is conserved at each vertex but not energy. 

If |g | 00, then in any finite order of perturbation theory the only 
diagrams which contribute have the form: 


= q 
P 


B 


since the transverse momentum cut off requires an infinite number of 
exchanges to connect the groups A and B. We may therefore write 
|F =|A> |B> and U|F> = U|A)> U |B}. 

Consider now a state |n) in U|p) which has n constituents with 


momenta: p=np+k, 
k;- p=0, yn=1. 
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Then 
‚lim E, -E,=(1/2p)(M? -k+ +MP)/n;)+O(L/P?), it O<n,<1 
= O (jpl) otherwise, (25) 


provided |k;| is cut off so that this expansion makes sense. The energy 
denominators in Eq. (24) therefore favour the intermediate states with 
0<n;<1 by a factor p? provided this is not compensated by factors 
of |p| from the numerator. It is possible to show [10] that such a compen- 
sation does not occur in the diagrams which contribute to W,, provided 
we restrict our attention to jọ and the component of j which is parallel 
to p. 

Similarly if we consider U |F = U|A> U|B) when |q |, we find 
that 


Eur ~ Er ~ O(4/\p|) + O(1/\n p + 41) 


in the leading diagrams. 

Hence, as |pl>c0 we may replace qg+E,—Ep by qo + Ey,—Eyr 
in the 6-function in Eq. (24). In this case, in frames in which |g ‚| and 
Ip: 
aoe Wyy = (4n)”! f d*x<Uplj,(x) j,(0) [U p> e'* 

pi? 
and we have “derived” the parton model: the bare current scatters 
incoherently (because of the transverse momentum cut off) from free 
on-mass-shell constituents (which have positive momenta along p) and 
energy and momentum are conserved at the elementary vertex. 

Drell and Yan [10] have analysed the set of frames in which this 
picture holds. The traditional choice is the ep centre of mass in which: 


do =(2v + q7)/4 |p| 
Igıl=V -g 


|p| and v= with œ =2v/(—q°) and qo fixed, but this choice is without 
deep significance, since the electron is merely a device to provide a virtual 
photon. 


b) Parton Model Calculations 


According to the discussion above the structure functions are given by 
adding the contribution of each type of constituent in the parton model 
and 


WO, q’) = f 2 U;(x) Wilx, v, q°) dx (26) 
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where U,(x) is the probability of finding a type i parton with a fraction x 
of the proton’s momentum p and Vj is its contribution to W, represented 
by the diagram 


xP+q 


We have put p',~ xp, since |p| oo and |k;| is cut off. The free particle 
kinematics yield the on-mass-shell 6-function ô(2 xq - p + q°) which gives 
scaling. 

More explicitly, for a free particle 


<pi| J2 (0) |p>> = Q;(p + p’), + spin dependent terms 
2 IF <FI= (On)? f |pr> d*prölpr — mF) <prl 


W}, =} Q? (2p + q).(2p+g),ö(g”+2q: p) + spin dependent terms. 


The spin dependent terms vanish as g—0, so that putting p,>xp, we 
find: 
W: = (1/x)[2.x° Q? òl +2xq- p)] 


where 1/x is put in to compensate for the different normalizations per 
unit volume of the proton (~2£) and the parton (~2E;=2x E) states. 


Hence 
vW? =} xu,(x) Q? 


(x=1/@= —q?/2v). 
In order to proceed we must decide: 
1. What are the partons? 
2. What are the u;(x)? 
Experiment already provides a clue about the spin of the partons. 
Consider the parton photon interaction in the Breit frame 


(27) 


-q7 2 


—— um 


9% 


If the parton has spin zero it carries no angular momentum in or out 
along q. Hence it cannot absorb a transverse photon and (er = or + cr): 


-— 7 


R=ds/67=0 (spin 0). 


If the parton has spin 4 its helicity is unchanged by the electromagnetic 
interaction in the limit |q| > oo and therefore, since its direction is reversed, 
it must absorb one unit of spin from the photon. Hence o,=0 in this 
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case and 
R=0s,/0r=0 


(28) 
or 2F,=wF, (spin 3). 


Experimentally it is known [3] that R is small and can be fitted by 
R=--q?/v? which is zero in Bjorken’s limit. Therefore the majority, 
and perhaps all, of the charged partons must have spin 5. 

Since we wish to consider models in which Gell-Mann’s current 
algebra hypthesis holds, it is rather natural to assume that the charged 
partons are quarks. We consider the results of this model in the next 
section. First we collect some results which follow from assuming the 
conventional weak quark current 


J =py,(1—ys) (ncosOc+Asin 8o). 
For AS=0 processes on free quarks this implies that: 
Wyn’? = 4 Tr(p+M)y,(1 +75) (@p+4+M) yl FY) 6(p +4)? —M?). 
For large v this gives 
Wy"? = ô(x + 2v/g?) 

vW"’P/M?=2xö(x +2v/g?) (29) 

v Wy"? /M? = F2ö(x +2v/g?) 
for the W; is Eq. (26). 


4. The Quark Parton Model 
a) A Difficulty 


An immediate difficulty with the quark parton model is that, according 
to the discussion above, the final state should consist of two separate 
groups of particles (A and B) each of which has non-integral baryon 
number and charge. Since quarks are (presumably) not being produced 
at SLAC, we must avoid this impasse or reject the model as irrelevant 
at present energies. A possible escape route is provided by the “wee 
partons” which carry a vanishing fraction of the proton’s momentum 
as |pl> 0. These “wees”, which we ignored above, have little sense of 
direction and therefore one which originates in group A may be assi- 
milated by group B, and vice - versa. Drell and Yan [10] argue that, 
although the wees play no role in finite order perturbation theory, they 
may be important if all orders are summed. Their suggestion derives 
some support from the work of Chang and Yan [11] who found that 
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in a 6° theory (albeit in a different situation) the “wees” contributed 
— ¥ —— =s%"'>0 as soo if G>1. 


For the moment we take refuge in a different philosophy — namely 
that we are using the quark parton model as the simplest model available 
in which scaling and current algebra obtain, from which we can abstract 
results which might be true more generally even if the naive discussion 
here is invalid. However, Dashen has advanced arguments which suggest 
that this difficulty persists in the formal approach, which may imply 
the existence of real quarks, but not necessarily that they are being pro- 
duced at SLAC; we refer to Ref. [12] for a discussion of this point. 


b) General Results of the Quark Parton Model 


In this section we will subscribe to the philosophy outlined immediately 

above and see what results can be obtained in the quark parton model, 

if we keep the distributions U;(x), which characterize the model as 

general as possible. The results will be rederived formally in Section 5. 
We already found (Eq. (28)) that 


2F,=0F), 


since quarks have spin 3. Another immediate result is that, since scale 
invariance and chiral invariance are both broken by the quark mass, 
chiral invariance will obtain in the scaling region. The tensor W,, will 
therefore be conserved, in which case: 


Fy’=0 2FP =F’. (30) 
(The first result follows from the positivity conditions for W,, together 
with Eq. (28). It is interesting to observe that these conditions give 
v Ws =0 if either Eq. (28) or Eq. (30) obtains.) 

We shall only consider the 4s = 0 structure functions here (the result 
for the As = 1 structure functions are given in Ref. [5]). Using Eqs. (28), 
(29) and (18) we see that we are left with 6 structure functions, which we 
may take to be: 

PL PNY BN Kar, 


which can be calculated in terms of the six distributions for quarks and 
antiquarks in the proton: up, Up U; Up, Up and uz (the distributions in 
the neutron are directly related by an isospin rotation). The A quarks 
do not contribute to 4s=0 weak interactions and in electromagnetic 
interactions they clearly contribute in the combination u, + uz. Effective- 
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ly, therefore, we have six structure functions determined by five distri- 
butions. We can evidently obtain one relation which is [13] +: 


12(F}? — FIN) = F3?— Fa". 


(31) 
(Sakata, Fermi, Yang: 4(F}?— F}") = Fy? — F3") 


as is easily checked using Eqs. (26)-(29). Here and below we give the result 
of the Sakata and Fermi Yang models to indicate the model dependence 
of our results. 

We already noted that u, + uz contributes to electromagnetic scat- 
tering only. Since u;20, being a probability, we can therefore obtain 
the useful inequality [13]: 


F}? + FY" S18(F]? + F}")/5 


32 
(Sakata-Fermi Yang: Fj]? + Fj"=2(F]?+ F?"). oe 


In the Sakata (Fermi-Yang) model the À is neutral (does not exist) 
and there is therefore an equality since the six structure functions are 
determined by 4 distributions. 

This exhausts the local relations. We now exploit the fact that the 
u, must yield the correct values for the proton’s conserved quantum 
numbers so that: 

Strangeness = 0 = f dx(u; — u3), 


Charge = f dx[$(u, a u;) —3(u, Ty u;) | ’ 
Baryon number = 4 f dx(u, + u, — up — uz). 


The first relation is useless in 4S =O reactions, since u, — uz is not 
measured. Using Eqs. (29) and (27), the second relation gives: 


[(E"— F3?) dx/x=2. (33) 


This is the scaling limit of the Adler sum rule [14], which depends only 
on the isovector nature of the weak current (and the validity of the deri- 
vation) but not on the nature of the constituents. Using the spin 4 relation 
2F, =a F, it may be rewritten: 


{ (F]"— Fy?) dx=1. (34) 


This is the Bjorken “backward” sum rule, which depends on the spin 
of the constituents. (The right hand side is zero in the spinO parton 
model and in the algebra of fields). 

1 In this paper the lowest moments | xdx of Eqs. (31) and (32) were also derived in the 


gluon model. At that time I was not brave enough to consider [3"J/2t", J] which gives 
the moments | x"dx for all n. 
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Using Eqs. (26) and (29) the baryon number condition gives 
(RP + 3") dx= -6 


. (35) 
(Sakata - Fermi Yang: 6>2) 


— a sum rule discovered by David Gross and myself [16]. The right 
hand side depends essentially on the non-integral baryon number at- 
tributed to the quarks; if it is correct it would verify the quarks algebra. 

So far we have not exploited momentum conservation which reads: 


1=fxdx} u(x) (36) 


according to the parton assumption that the constituents all move in 
the same direction in the |p|—>œ frame. The difficulty in using Eq. (36) 
is that the right hand side contains contributions from neutral gluons 
which may be present (and serve to bind the quarks together), whose 
distributions cannot be measured, since they do not contribute to 
lepton scattering. However, we may separate their contribution and 
write: 
6 


rar) (37) 
OSeS1 


where e is the fraction of the proton’s momentum carried by the gluons. 
Eq. (37) may be used to obtain the interesting upper bound: 


1 
| (F3? + F3") dx <3 (38) 
0 


which is satisfied by the data (unless big surprises await us at unexplored x) 
the left hand side being ~0.32. 

A question which arises is whether e is zero. A possible test involving 
electromagnetic data alone is provided by the relation [5, 17] 


FR? dx2(1-— e)/9 (39) 


— the left hand is certainly >$ for the proton and very likely also 
for the neutron (in the model of Bjorken and Paschos [9] the integrals 
are 24 if no gluons are present, but this requires more specific assump- 
tions about the u,(x)). 

Finally we note that we can derive the relation [4, 17, 18] 


e=1+ [dx(3 F3r+"/4— 9 F3? +"/2). (40) 
The result of the CERN experiment gives [6, 7] 
o’P+’"2 = G? M E (0.52 + 0.13)/n 
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and hence 
{ Fr +" dx 21.08 +0.27. (41) 
Inserting this in Eq. (40) and using [3] 
{ F3? +?" dx = 0.28 + 0.04 (42) 
we find 
620.52+0.38, (43) 


i.e. the present data seem to require gluons. 

Eqs. (41) and (42) are compatible with (32). In fact the model is quite 
compatible with all present data [5] as is illustrated by the more specific 
model discussed by Landshoff in his lectures at this summer school. 

A comparison of Eqs. (41) and (42) with Eq. (32) shows that the Sakata 
and Fermi-Yang models are excluded by the data. In fact any model 
with integrally charged constituents will tend to. give too small a value 
for F3/F}. F} depends on Q? while Fy depends on the “weak charges” 
which are zero or one, depending on the isospin but not on Q;. Therefore 
we expect that in most models with integral charges F>/F} will be less 
than in the quark model. Since the data (Eqs. (41) and (42) is already 
approximately equal to the maximum value permitted by the quark 
model (Eq. (32)) we see that models with integral charges are likely to 
be excluded. This has been shown explicitly by Nachtmann in several 
cases [30]. 

Nachtmann has also derived several new inequalities in the quark 
parton model [30]. The most interesting is 

ı- FE") _ Su, tu) +5, +u +u, +u) 
4 


NA 


<4 


Fo) (Up + Up) + S(Un + Uz +u, + u) 
which follows from the fact that u; 2 0. The latest data [31] give a value 
of the order of 4 for this ratio (albeit with large errors) for x ~1. More 
accurate measurements of this quantity will be very interesting. Nacht- 
mann has derived further inequalities which depend on the fact that 
the nucleon belongs to an isodoublet or that it belongs to an octet 
(assuming exact SU(3) symmetry). All these results can be rederived 
using the formal methods outlined in the next section or (equivalently) 
by means of light-cone commutators. 


5. Formal Derivations 


The formal derivation which we shall present starts from the Bjorken- 
Johnson-Low (BJL) expansion [19,20] of T,,(Eq.(21)). In the limit 
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go>i%, the exponential provides a very rapidly convergent factor 
exp —|go| fin the t-integration. It might, therefore, be sensible to expand 
the commutator in a Taylor series about t=0 which gives 


T, 222 Ý (i/go)"*! f dxe'4* Cpl [J (x, 0/80", J,(0)] [p> 
n=0 


(44) 
+ polynomial in go 


provided the commutators exist. Eq. (44) is the BJL expansion. It is at 
this point that we part company with perturbation theory in which [21] 

1. The equal time commutators do not all exist. 

2. Even when they do exist, the values obtained using Eq. (44) do 
not coincide with the values given by naive canonical commutators. 

However, as explained above, we assume that scaling occurs and 
this gives us courage to proceed blindly using Eq. (44) and naively mani- 
pulating operators according to the canonical rules. 

In order to avoid inessential complications we consider the case of 
scattering a scalar current in which 


T =i | O(xo) e'*’* <p] [J (x), J (0)] [p> d*x + polynomial in go 


2 Wi’, q*) dv’ 
= T(0, q?)+2v f AI my 
2 wo Wt+aWw 


To, Orga oot a ea 
(w, q?)= T(0, q?) + al SATa do’ + (45) 


wrrriw*r w2"t2Ww 


10,9 +22 ay 2n+2 =, POETE Jao + 


W*=Wo, q’) W(-w,q°), lol<1. 


Next we can take the limit gy >i% in a frame with q = 0 and compare 
the coefficients of (1/qo)" in Eqs. (45) and (44) (assuming scaling), following 
which we divide by (po)" and let pọ ©. More directly; we can let gy ico 
with w= 2po/(— qo) fixed (w| <1) and compare coefficients of œ which 
gives: 


2 | F*daw/o"*? = ‚im, i(—i/2 po)" * f dx <p] [3 J(x,0)/8 t", J(0)] |p> (46) 
~— ifn odd, + if n even, F(w)= lim W(w,q?). 
© fixed 


In the case of vector currents the moments of different structure functions 
are given in essentially the same way by the commutators of different 


68 C. H. L. Smith: 


components of the currents e.g. 


2) FË do/a" +? = „im i(i/2po)"* + f dx<pal [FJ (x, 0/2 t", JZ (0)] Ip2> 
o > 00 
F : > 47 
— if n even, + ifn odd, FE =F +F}. er 
The complete set of these relations is given in Refs. [4, 17]. 
In order to calculate time derivatives of the currents we must specify 
the interaction Hamiltonian, which we take to be the sum of the re- 
normalizable interactions: 


Hr=g PP +g PPYP + GBP. (48) 
Formally we can construct the commutators 
f dx[0"J (x, 0)/62", J(0)] (49) 


completely using H;. However, we need only pick out those pieces 
whose diagonal matrix can grow like p&* ' (or faster) when py— œ, since 
only these pieces contribute in Eq. (47) and the analogous equations 
for the other F,. The pieces must obviously be parts of Lorentz tensors 
of rank n+1 or higher. 

For simplicity we consider first the case g,=0. We note that the 
equal time commutators never introduce inverse powers of masses or 
fields. Therefore, in Eq. (49) we have an operator of dimension M"*? 
constructed out of positive powers of masses and fields, from which we 
must pick pieces of tensors of rank 2n+1. Noting that the operator 
must contain at least two quark fields, it is easy to see that the only such 
operators are [5, 17]: 

PO) Vas Baa «++ Gags s PO) 


BO) Lays Yard Pag +++ Fans 2'P (0) 


Therefore we may set g, = g, = M, = M,= Mọ =0 in calculating the 
operators. Hence we have derived the parton model in this case, since 
we have shown that all moments of the structure functions (and hence 
the structure functions themselves, barring pathologies) are given by 
the matrix elements of the same one-body operators as in a free field 
theory of massless quarks. 

We have treated the case g, +0 by this method elsewhere [17]. 
Again it turns out that we can put M,=0. Because of the gauge in- 
variance it is then hardly surprising that we need only change the free 
field theory results by putting iö,—iö,+g,B, everywhere, which does 
not change any of the parton model results which depend only on the 
SU (3) and tensor properties of the operators. We do not give details 
here since the vector case will be elegantly treated by David Gross in 
his lectures at this school. 


(50) 
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The discussion above reestablishes all the results of Section 4 except 
the sum rules and inequalities based on momentum conservation. To 
derive them we look at the explicit form of Eq. (47) in the case n=1: 


| Fldx= lim (1/273)<p.16O) ((7.0.+9472B.)0?vOp.> (51) 


[E+ yndx = „lim (1/2 p%) <p. p0) (iyz0,+ 9,7zB:) (4B+2 Y) y(0) IPz> 
where Q, B and Y are the usual SU(3) x SU(3) matrices and 
Q? =2B/3+ Y/6+1;/3. 

An important point in the following is that the matrices B, B— Y 
and 2B+ Y+2I, make positive semidefinite contributions whenever 
they appear on the right hand side of Eq. (51) [13] (this is needed in the 
formal derivation of Eq. (32)). If we call one of these matrices A, then 
we can consider a function F? defined in terms of py,Ay just as F} is 
defined in terms of #y,Qy. In the analogue of Eq. (51) for F}, Q? will 
be replaced by A?«A. The left hand side is positive semidefinite and there- 
fore the right hand side must be so also (in parton language this corre- 
sponds to the assumption that 4,20). 

Next we note that 


‚lim (1/2 pt) <p- PO (i7.8.+967-B.) WO) Ip.) 
= lim (1/2 p3) <p- 9- — 68, — 1.) 
=1-:20 
e= lim <p.| 6%, 1p:2/2 Po 


where 0,,=6£/ög,, is the symmetric energy momentum tensor and 
6%, is the part of 8,, due to the free gluon fields. Clearly e <1, since the 
left hand side is positive semidefinite according to the argument above. 
If we could establish that &> 0 this would give all the results at end of 
Section 4. Naively this seems to be the case since 6%, is a sum of squares. 
However, this is not a sufficient condition since the subtraction of the 
vacuum expectation value is understood here (and throughout) [23]. 

For (pseudo-)scalar gluons we can consider scattering the current 
0, and show that their contribution to ¢ is positive semidefinite using 
an argument similar to that in the paragraph following Eq. (51). For 
vector gluons we have not yet been able to establish positive semi- 
definiteness, although the parton analogue (that all constituents move 
the same way in the infinite momentum frame) seems to be true in per- 
turbation theory in a cut-off field theory [24]. 

To summarize: The “momentum sum rules” at the end of Section 4 
can be derived in models with scalar and pseudoscalar gluons. It is 
possible that they can also be derived with vector gluons. 
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6. Miscellaneous Remarks 
a) Other Processes 


We refer to Refs. [10, 25—27] for a discussion of the application of parton 
models to some other processes. 


b) Parton Models and Local Duality 


It has been suggested [28, 29] that many of the properties of the structure 
functions near w =1 can be calculated in terms of the behaviour of the 
Born term as Q?->00. I am sometimes asked how this could be under- 
stood in parton models, since they deal with incoherent scattering, while 
the Born term involves the coherent process of elastic scattering. The 
answer is that it cannot be understood and in fact the parton model 
(or the formal calculations using [J, J] commutators) is incompatible 
with these calculations, as can be seen by considering spin $ parton 
models. 

Consider deep inelastic en-scattering. The arguments which are 
appartly successful [29] in the ep-case (see the accompanying lectures 
` by Rubinstein) can be applied giving: 


o Born term 
1. R= pen ~ = 
la | Q’-00 > 
since or=0 for elastic scattering off a spin 0 particle. 
1 ge vw u noe term) _ 0. 
vW?” loxi Q*> x 


In contrast, R=0 in spin 4 models and S is presumably not zero (it is 
one in the quark model and in all models in which there are no local 
isotensor operators with “twist” — equals dimension minus spin — less 
than or equal to two). Therefore the extreme use of local duality in 
inelastic lepton scattering is incompatible with constituent models or 
formal calculations involving the commutator [J, J]. 
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Introduction 


The purpose of these lectures is to review and discuss some known and 
new results concerning duality and scaling. Particular emphasis is given 
to questions concerning the possibility that electroproduction and photo- 
production are more closely connected than expected, or more precisely, 
that the limit g?0 is very smooth. 

Most of the considerations described here are contained in papers 
written in collaboration with Rittenberg. 

In the first part we review shortly the kinematics and notation. We 
discuss a few general concepts, like the relations between Regge and 
Bjorken limits. 

In the second part we introduce the concept of duality and discuss it 
in the context of fixed q? sum rules as done originally by Bloom and 
Gilman. 

We emphasize the difference between strong and variable q? processes. 
In the latter, duality has a meaning irrespective of the model one uses 
for the asymptotic form of the amplitude. 

In the third part we discuss analyticity and local duality in the whole 
v—q’ plane. We arrive to these ideas, or some of them by different 
roads (Regge, light cone and finite energy sum rules). 

In the fourth section we discuss the extent to which duality is locally 
obeyed and introducing an appropriate scaling variable we find a 
remarkable connection with photoproduction. Properties of the neutron 
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proton structure functions difference are analyzed. Questions concerning 
the time like region are briefly mentioned, and predictions that will 
test the model in depth are discussed. 

We conclude by discussing the relevance of these results under the 
light of current available models. 


Kinematics and Limiting Properties of the Structure Functions 


The double differential cross-section for inelastic electron scattering for 
an inclusive process in which only the electron coordinates are observed 
is given by [1]: 
do _ 40°? E? 
dQ’dE a 


[2 W,(v,q7) sin? 0/2 + W,(v, q°) cos? 0/2] (1) 


where E’ and dQ caracterize the energy and solid angle of the outgoing 
electron, 0 is the angle measured with respect to the beam direction, 
q’ is the mass of the photon, that in these experiments is kept spacelike. 


v=E-E. (2) 


E being the initial electron energy. 
The aforementioned quantities are related by the equation 


q’=4EE sin? 0/2. (3) 

We also compute the missing mass to be 
s=W?=2mv+m?- gQ?, (4) 
u=—2mv+m’—q’, (5) 


where m is the proton mass. 

The main difference with hadronic reactions is that the external 
mass can be varied. The v— q’ plane can then be studied by properly 
defining lines as seen in Fig. 1. Constant W? lines have positive universal 
slope and negative intercept, while fixed angle lines are given by the 
equation 


q =4E? — 4 Ev sin? 0/2. (6) 


We also introduced constant wy lines defined by 


2mv +m? 
2_ Be) 
q OE a (7) 


where a? is a constant. 
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Fig. 1. The v—g? plane. The quasi-horizontal lines correspond to #=6°, E=7, 10, 13, 5 
and 16 GeV discussed in Ref. [2]. Dotted lines correspond to fixed w’, the other to fixed 
s-values 


The structure functions W, and W, are the absorptive parts of the 
forward Compton scattering amplitude (spin averaged). The formulae 
reads 


W, =[(v—q’/2m)/4n?«] or, (8) 


2 
W, = [0 q?/2my4na] i (er +s) 0) 


where Gg is the scalar and cy the transverse photon cross-section. Notice 
that W, is kinematically constrained to vanish at q?=0. This simply 
reflects the transversality of physical photons. Since oy is finite for 
q? =0 we have 


wW, 1 crl =0) 


limg’—0 = 0 10 
==; q? v Ana = (10) 
an expression we will use later. 
Another quantity of importance is 
2/2 2 W, 
patea ta Ema |. (11) 
Or W, qty 


The Regge and Bjorken limits are different asymptotic domains of the 
structure functions. It has been conjectured that they may be related [2] 
though there is at present little or no evidence that this is the case. 
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We first keep q? fixed and consider 
lim v 00 W, ~ B; (0, 2) vH O(v* 2), (12) 
lim v> W, ~ B,(0, q?) v9"? + 0 (v2 9), (13) 


is the t=0 intercept of the leading trajectory that couples to forward 
compton scattering like the Pomeron or A,, for example. f’s are residue 
functions and we have not put, notice, a Regge scale factor. The ex- 
ponents are different because W, and W, have different number of 
helicity flips. 

Since these functions are absorptive parts, contact terms and fixed 
poles make feel their presence only indirectly and through dispersion 
relations [3]. 

The Bjorken limit obtains when 


q? and v>o = = = Wy = finite. (14) 


Scaling is the statement that these limit exist 


Wi <=, Fi), (15) 
VW, ssw? F2(o). (16) 


In perturbation theory this conjecture is not verified [4]. One obtains 
scaling up to logarithms. This problem is to be kept in mind. A stronger 
connection has been proposed by some authors [2]. The same term 
given by (12) and (13) that is leading in the Regge limit is leading in the 
Bjorken limit. This implies 
Cg lea Mae (17) 
Bla) (a). (18) 


This behaviour obtains in simple ladder theories [5] (neglecting loga- 
rithms!) though counter examples exist as well. We will not assume 
these connections in our analysis but it will be helpful to use the language 
here and there. Using light cone analysis it is also clear that the Bjorken 
and Regge limits obtain support from different regions in configuration 
space [6]. Not surprisingly, the Bjorken limit tests short distances while 
Regge distances are related to the masses of the hadrons involved, or 
sometimes to the universal slope 1 (GeV)?. 


Fixed q° Duality 


We will assume in what follows that Bjorken scaling holds for the proton 
data to a reasonable degree of accuracy (see Fig. 2). Hence, we assume 
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Fig. 2. Values of W, and v W, plotted against œ for W = 7, 2 GeV and varying q? for R =0.18 


(12) and (13) to be valid. For large q% and v, in the neutron data this is 
satisfied provided some additional assumptions are made [7]. 

Duality in electroproduction carries slightly different connotations 
as compared to purely hadronic processes. In the latter case one makes 
a statement about analyticity of the amplitude under discussion and 
then extrapolates at low energy the asymptotic form of the amplitude. 
This last step demands some model, usually a Regge parametrization, 
and a saturation hypothesis, allowing to cut the integrals at low values. 

In photon initiated processes one does not need to do that. No 
model is needed to determine the asymptotic form of the amplitude 
and to extrapolate it. At fixed œ a point can be reached either by high 
v and q? or both being proportionally smaller. Hence independent meas- 
urements test the equations directly. 

So, no hypothesis about the nature of the J plane singularities 
controlling the Bjorken limit is really needed. A fixed q? sum rule can 
immediately be written using the hadron techniques to read 
1+W?/q? 


2m e 
LM f ymo, qg) dv= | do'F,(o’ (19) 
1 


0 
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where we have introduced œ' as defined in the coming section. Here 
we do face an ambiguity. Though the asymptotic form œ of the variable 
is fixed, the variable that averages the low v, q? data is not prescribed a 
priori. The possibility of using expressions other than œw as a scaling 
variable was pointed out at SLAC [8] and since then and even before 
other improvements have been suggested. 

If we look at Figs. 3-7 we see the typical behaviour of some data for 
q? = 1.0 up 2.5 from Gilman’s Tel-Aviv talk. 

The resonances very much follow the scaling curve and, what it is 
striking, is that each peak is changing and adjusting itself to the new 
scaling curve. The scaling variable œ is given by (2mv + m?)/q?. 

These figures could be drawn for W, as well. Traditionally vW, has 
been preferred because at small angles W, could not be measured accu- 
rately because of the sin? 0/2 depressing factor. 
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Fig. 3-7. v Wz(v,q°) plotted against a’ for fixed q? as stated. The solid line is a smooth fit 
for q? > |W) 2 GeV. The arrow indicated the position of the elastic peak (see Ref. [is}) 
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The fact that the scaling curve averages the resonances to 10% from 
the very beginning is quite remarkable. This is precisely what local 
duality means: locally the asymptotic form is equivalent to the low 
energy amplitude. 
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Historically, early theoretical models [2] conjectured that the re- 
sonances could disappear fast when q? grows and that the scaling proper- 
ties witnessed at low q? reflected properties of the background or, if 
we use Regge language, the diffractive Pomeron. We have just seen that 
this model is ruled out by the data. The agreement between both sides 
breaks down in the absence of resonances. We will come to that later 
on again. So the next theoretical question is: can one produce models 
in which direct channel resonances build a substantial part of the vW, 
function. This is indeed the case and some of these dynamical models 
[9] are discussed by Landshoff in his lectures. 

We want now to discuss the predictions of the local duality model. 
We warn the reader that strict local duality is of course untenable as 
resonances peaks do indeed exist. Nevertheless the narrow resonance 
model is used to show that its predictions are semi-quantitatively correct. 
The sum rules can be written down to read 


1+W?/q? 


2m * 
ae Savy Waly) = | davw). (19) 
0 1 


Strict local duality demands that the integrands are roughly equal point 
by point, or say resonance by resonance. We parametrize the excitation 
form factors of the system for large values of q? as a power law that we 
assume is universal, then 


Ga) VI”. (20) 
Near threshold we have that 
yW, > (w — 1} (21) 


since the structure function vanishes when w—1, dynamically. By 
evaluating the integral at a resonance locally, and remembering that 
the contribution to W, ~ô(W?— m?) G?(q’) we obtain the relation, 
keeping q? large: 


n=p+l. (22) 


This relation holds quite well experimentally [10]. For the nucleon form 
factor n=4 and the fall off to threshold of vW, seems to be consistent 
with p= 3. Data however does not determine the slope very near thres- 
hold, only up to œ = 1.5. The argument implies that the relation holds 
universally for all resonances. This is also compatible with the data [8]. 
The corrections one would expect due to finite widths are not easy to 
estimate but the results seem quite encouraging. 

Similar arguments can now be applied to obtain an estimate of the 
ratio of the neutron and proton structure functions. Since the normali- 
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zation of these form factors is known one integrates up to about the 
nucleon pion threshold (Wnresn). Taking the derivative of both sides 
of (19) and keeping always q? large one obtains 


Wo =1+ Wiiresn/G) = (0 — 1) (—g? dG? (q?)/dq?) (23) 
where G(q’) is essentially at large q? the magnetic form factor. 
C@) aap GP V7 + Gila’) (24) 


It immediately follows that [8] 


2 
p =r (=) =0.47 (25) 


P 


where one has assumed the same functional form of the neutron and 
proton form factors for large q?. This result is compatible with the present 
date but one only has values that are not that near threshold. It is inter- 
esting that recent parton models [11] have predicted the ratio to be 


(| =0.66. (26) 
P 
This value is not ruled out but seems not be preferred by experiment. 

To conclude this section we should like to emphasize some crucial 
points that seen by now definitely established. 

a) Resonances survive in the large q” regime. 

b) Though they may constitute the whole y W, at small w’, back- 
ground cannot be ruled out, specially for high v. 

c) At least for fixed, relative large q? duality holds locally in the 
sense that after every resonance the sum rules holds to better than 10% 
in accuracy and that their presence is essential to get agreement in 
Eq. (19). 

d) If one assumes a Regge picture and this is something unrelated 
to the previous points, the presence of resonances demands the existence 
of “normal trajectories” and at most some diffraction. The study of 
Regge-like models and how they fit is beyond the scope of my lectures. 
However, I feel that there is much to be desired about their ability to 
explain the data. We will come back to this point in our last section. 

Llewellyn Smith has also pointed out that if local duality holds for 
pion targets then the Born term predicts a very large (really infinite 
05/07 ratio). This of course contradicts the quark model algebra. I think 
that if local duality holds there too, the conclusion is inescapable and 
would be a point where the marriage of partons and duality may irrevo- 
cably dissolve. 
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Duality in the Whole v—q? Plane 


As stated in the introduction varying q? one could learn much more 
than in the purely hadronic case where the external masses are fixed. 
Since duality is not rigorously local and form factors are rapidly varying, 
new information can be gained by writing sum rules along lines of the 
form av +q?=b where a and b are constants. Two problems must then 
be solved. One, is to establish enough analyticity so as to write the re- 
levant sum rules. Second to make some dynamical assumption so that 
the sum rules gain content. We first discuss the problem of existence. 

The forward scattering amplitude for massive Compton scattering 
(neglecting spin complications) can be written in three different ways 


T,= Set *dxé,<p| Da), 3,0] Ip> (27) 


where &,=0(x,) for the retarded function, 
€;=6(—x,) for the advanced function, 
é,=T, time ordered product for the scattering function. 

All of these agree in the physical region since they differ by terms 
that vanish when the four momenta are physical. This is so due to the 
stability of the one nucleon state. However, whenever an analytic 
continuation for complex values of v and eventually q? is performed the 
results will be, of course, different. 

To study analyticity one should use some representation that ex- 
plicitates singularities. We use for this purpose the Deser-Gilbert- 
Sudarshan representation [12] (DGS). This representation is based on 
microcausality and some technical assumptions about the particular 
model under discussion. Nakanishi has proven that every connected 
Feynman diagram satisfies a DGS representation, and these Feynman 
diagrams are in turn connected to the expression (27). Hence we are able 
to establish results to all orders in perturbation theory. This is perhaps 
not enough because of what we said before about the Bjorken limit [4] 
but it is the best we can do. A derivation of the equations can be found in 
the Boulder Lectures by Brown. We simply quote the results and outline 
the procedure to study the analytic properties. The retarded or advanced 
functions read: 


h(a, p) 
(-q°+B—otie(v/2+ p) 


Tp = | ap fac (28) 


=j 

where we have set M =1 and h(o, ß) is a weight function whose support 

is established from the spectral conditions from the commutator. 
A slightly different representation holds for T, 

En: doh(o, ß) 


= J ae (-Q?+Bß-o+ie) 


(29) 
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so that 
1 œ% 
—nI,T = J dB) dod(B-—o—q’)h(B, o). (30) 
ai 


Hence the problem is to see when the line 
-q +B-c=0 (31) 


intersects the support of h(o, ß). For v and or g? fixed the analyticity 
is easily established. If the mass g? is timelike then after the first threshold 
the v axis becomes fully cut. Real analyticity is lost and some procedure 
of analytical continuation must be proposed. The new function will 
not agree either with T, or Tp. If the vectors q and p are real the support 
in the physical region of h(o, ß) is related to a parabola v?=4q’. Any 
line av+b=g? not entering the region bounded by this parabola has 
a gap in the v axis and hence a standard finite energy sum rule can be 
written. Otherwise the problem of continuation must be faced. Of course 
these are the most pessimistic possibilities. Models need not have these 
singularities. The generalization to lines that avoid the parabola for 
lines of the form av+b=g? was first proposed by Leutwyler and 
Stern [13]. 

However, in general the sum rules can be written and they will read 
ive) Vb 
{[vW,(av+b=q’)— f(@)]dv+ J VW,—fl@))dv 
Va = 0 


(32) 
+ f DW-flodv=0 


complex cuts 
where the last term reflects the possibility of having to integrate over 
complex cuts. 

Fortunately we need not worry because we add now the hypothesis 
of semi-local duality [14]. This dynamical assumption states that these 
integrals must vanish by themselves since the amplitude equals locally 
its asymptotic part. Hence, we have established that, if semilocal duality 
holds, finite energy sum rules are expected to exist in any direction in 
the v—q? plane. 

Now we want to check whether these new equations check experi- 
mentally. The sum rules originally tested by Bloom and Gilman were 
the fixed angle sum rules. These read 


E (v (om) — Flow) (q? a. a?)dv . (32a) 
jp (2m+4Esin”0/20,) 

Hence, a logical starting point, partly determined by the available data 
is to choose a direction perpendicular to a small angle. This way a fixed 


« line of about 2-3 is the best choice. 
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We combined the data of a large number of experiments since these 
are performed at either fixed missing mass or fixed angle. After processing 
the data to reconvert it for our purposes we plotted v M, = vW,/f(w') 
as the ordinate against a variable roughly proportional to v. It reads 


¥=(v+m/2) (w — 1) = W?/2m. (33) 


This variable has the advantage that whatever œ the resonances appear 
at the same value of y. This is needed to test properly local duality. 


2.0 Cpt yee pe org 


0.5 1.0 1.5 2.0 2.5 
v(Gev) 


~ vW, A ; 
Fig. 8. Fixed œw sum rule. y W, = Tor is defined by Eq. (3). The origin of the experimental 
w 
points is given in the upper corner of the figure. The first column corresponds to 2 <a’ <2.5 
the second one 2.5<w'<3. The curve is traced by eye. Representative errors are shown 
for two points. The position of the resonances is shown; for the first three of them their 
width is given 


Again to compare points with different & we divided v W, by the 
value of its asymptotic form. For this region of œ the variations of the 
asymptotic form are still considerable. This way all points are plotted 
against a universal function. Notice that we choose w’ the same variable 
as in the previous section. We will show later on that there is a better 
one Wy but they only differ for small q?. Hence, to avoid ambiguities 
we restricted ourselves here for large values of q?. The acceptance Aw’ 
was rather large to be able to gather enough data. We did separate the 
points between 2 and 2.5 and 2.5 and 3 but no systematic variation was 
observed. We finally plotted all points together. By looking at Fig. 8 
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we first observe that the local average obtains to a very good approxi- 
mation when we plot the resonances against œ. œ does not average the 
resonances in this direction either. 

If we cut the sum rule at the bottom of the valley after every resonance 
we find that both sides are equal to about 10%. 

The contribution of the nucleon pole was computed explicitly. 
However, its local contribution is undefined and we arbitrarily spread 
it over a typical hadron width. The curve was traced by eye aided by the 
known resonances widths and positions. We have no way to estimate 
errors in a reliable fashion and hence we took seriously the ones quoted 
in each experiment. Under these circumstances the agreement seems 
beyond expectations. There are several qualitative aspects of the picture 
that deserve attention. First we notice that in Fig. 8 and 9 the relative 
contributions of the A (1236) is changing relative to the next resonances. 
This indicates that the form factors of higher mass resonances have a 
higher tail than A but at the same time there is a crossover point. 

In Fig. 8 all resonances contribute with about equal strength. In 
fact for Ÿ up to 1.5 the prominent resonances build most of vW,. Small 
contributions from daughters may make the resonance contribution 
almost 100%. This possibility has been also noticed by Gilman [15] 
in the fixed q? case as well. When Ÿ grows the interpretation becomes 
less compelling and nothing can be said. Hence one must conclude here 
again that resonances scale and they are a big piece of vW.. 

To further exploit those sum rules one should continue and use the 
left over equations. Going to the nucleon pole in the u channel and 
fixing m= —gq? we compared this contribution to the scaling function 
and found it far off. 

This result is perhaps not surprising since there is no reason to expect 
the mass singularities to be dual to the intermediate hadronic states. 
This may be one reason why models for currents including duality are 
still facing difficulties. One way to see how individual resonances may 
scale in theories (dual models) where there is a growing caracteristic 
mass is the following. One worries about having large masses in the 
theory since the trajectories raise continuously. However, if the current 
couples through vector mesons also through members of these trajectories 
the universal slope is the same in the current and t channels and the 
scaling obtains as the slope cancels [16]. 


Extension of Local Duality as a Function of q? 


We want now to extend duality and scaling up to q?=0. If we look at 
Fig. 9 we see that the averaging performed by w’ is not very good for 
low values of q? for v W, or the experimentally measured cross-section. 
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Fig. 9. a, b plot of v W, versus both against œ (dashed dotted line) and wy for a?=0.15 
(dashed line). The smooth averaging curve corresponds to high q° and v data assuming R =0. 
d?o 
dQ'dk’ 
tively, going upwards, to a? = 0.25; 0.15 and zero 


c, d plot of against W. Data from Ref. [13]. The averaging corresponds respec- 


Because of kinematical factors and the existence of both longitudinal 
and transverse contributions the two tests are not equivalent. So any 
averaging variable should affect only low q? values. We propose 

_ 2mv+m? s+q° 


= = : 34 
Ow gta gta (34) 
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Fig. 10. 2m W, plotted against &,. Data from Ref. [13]. Circles are points calculated from 
data from Ref. [5]. Rhomboids correspond to the 1520 MeV resonance, Ref. [7] 


We first assume R =o;/o7=0 (we will remove this assumption later) 
and check local duality with a? = 0.15 —0.25 and notice that local duality 
is much better verified. We cannot adjust better a? with the data we 
had but an effort to test the formula with the world compilation of data 
is already under way [17]. 

Now we can combine local duality much as in the previous section 
with our variable and make further tests. We consider electroproduction 
in the region of the A (1236). Data exists for q? = 0.1 [18, 19]. We computed 
W, from the experimental transverse cross-section, which is known 
separately by averaging the values at the top of the resonance and at 
W = 1390, corresponding to the center of the subsequent dip. This is 
of course a rather brutal, first order approximation to averaging, but 
if duality holds to about one resonance width it should give a reasonable 
answer. The data are very scarce (6 points!) but we can use both fixed q? 
and fixed wy sum rules. The results are shown together with the ex- 
perimental data [7] for R=0.18 in Fig. 10. The agreement is really ex- 
cellent. A few points coming from Ref. [20] are plotted too and they also 
fit well. 

We next compute R=<(os)/Xor). We assume that the ratio itself 
is smooth. Values of the order of 0.2 are obtained from the 4(1236) 
including points with q? = 0.1. Unfortunately errors are too big to allow 
for a meaningful determination of R as a function of wy. It is again no 
small triumph of duality that the resonance value of R is the same as 
obtained at high energies and high g?. 

The next excercise is to discuss v W,. Here we have as shown in Eq. (10) 
a kinematical zero. Hence we must modify the structure function to 
eliminate this effect. The simplest way to eliminate this “threshold effect” 
is to consider [14] 


(v W,) i, = @/@y(v W2). (35) 
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We can test our variable against the SLAC data and see whether the 
threshold change of v W, as a function of œ as explained. The agreement 
is excellent. In fact we predict that for 


w =m ja? (36) 


the slope will change sign. This seems to be indeed true. The results 
are shown in Fig. 11. Unfortunately, there is some averaging over œ 
that precludes an even stronger test, but the œ dependence of the slope 
is beautifully explained. We have found that properly interpreted the 


j 16<w<24 

re el 

o | 23 45 0 I| 2 3 4 8 
q2(Gev/c )2 a2 (Gev/c)® 


Fig. 11. v W, behaviour as a function of q? for fixed averaged values of œ. From Ref. [7] 


whole electroproduction amplitude is a function of one variable through- 
out the whole v and q? range. The A resonance contains all the informa- 
tion about electroproduction. This is not completely true since there 
is a bound on wy of about 6. This makes the connection with diffraction 
satisfactory since at the 4 the background is small and the high v limit 
seems diffractive in nature. 

Our connection between duality and scaling seems to hold. It is a 
very strong restriction on theories. Available dual resonance models do 
not seem to have the local averaging variable as demanded by the data. 
By far the most intriguing aspect of the problem is the ability of the A 
resonance alone to reproduce the asymptotic scaling function. Narrow 
resonance models in which other trajectories are to be included as 
additive terms have no contribution to the imaginary part at the A 
position. Hence, it might indicate that one needs more complicated 
than additive models to explain the data, since the A by itself by varying 
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q? reproduces the whole structure function to which these other trajec- 
tories contribute as well. Hence duality seems here even stronger than 
in the strongly interacting case. 


Connection with Photoproduction 


In the previous section we have extended duality to q?=0.1. One can 
immediately ask if one can take the limit g7=0 and make a fruitful 
connection with photoproduction. 
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Fig. 12. 2m W, plotted against wy. Data from Ref. [13]. The curves are calculated from 
photoproduction data, Ref. [14], for a7=0.15, 0.20 and 0.25 respectively starting below 


To that effect we can simply compute the structure function W, 
by taking the total cross section data [21] and going by the same technique 
of averaging from top to bottom of each resonance. The results are 
satisfactory but very sensitive to the value chosen for a? (see Fig. 12). 

The possibility of using real photon data allows us to make some 
striking predictions concerning the neutron proton structure function 
difference. 

The argument is very simple. The total inelastic cross-section for y 
rays on protons or neutrons near threshold is the 1 pion production one. 
At threshold one has a low energy theorem, the Kroll-Ruderman theorem, 
that states that 
(37) 


oO (o 


yp>pr? = Oynonn0=0> 


Oy ponnt — Fyn pre + (38) 
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Hence the total cross-sections are equal. Using formula (8) or (35) if 
one wants vW, we predict that around that value the difference must 
vanish. Of course, duality is not strictly local but a detailed study of the 
multipole models that fit the data very well [22] demand in fact that 
throughout the first few hundred MeV o,,,20,,. Hence we predict that 
the difference between neutron and proton must vanish or become 
slightly negative at the point wy ~ 6-7, where ¥ is some value around 
the first resonance. 

Though the analysis is slightly model dependent if one accepts the 
predictions, at the second resonance the proton has taken over again. 
See the accompanying Fig. 13 [23]. 


E (GeV) 
02 04 06 08 10 
00 


o(ab) 


ce) 
13 I7 21 25 


E"Ticevy?] 
Fig. 13. Total cross section for the reaction y+ n>p+n? as a function of incoming energy 
given by ++++,y+p>n+n* given by ——, y+p>p+n® given by —-—-— (see 
Ref. [23]) 


The value here for wy is about 8. For large values of œ there are 
data [21] and the difference between the proton and neutron does tend 
to zero. 

If these oscillations really take place, then the usefulness of the Regge 
parametrization at moderate values of œ can be put in doubt. 

Though it is too early to decide because of problems with the data, 
more tolerable errors in this difference will be of importance. 

If one goes to electroproduction with small q? at the 4(1236) re- 
sonance, one moves to smaller wy and hence one expects, if the theory 
works, that the proton cross sections will grow, relative to the neutron 
one, with q?. There may be some evidence already in favour of this pre- 
diction. 

The most puzzling aspect of this work is the origin of a”. A detailed 
study of the data may allow for a? slightly bigger, if one allows the m? 
parameter in the numerator to change as well. It is perhaps no accident 
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that both constants are roughly the two numbers that seem to charac- 
terize strong interactions: the transverse momentum length and the 
universal slope of Regge trajectories. 

From the theoretical point of view the breaking of scale invariance 
seems to be characterized by a small number a?. There have been some 
attempts [24] to relate by use of Regge theory photoproduction and 
inelastic electron scattering. This model does not predict the proper 


Fig. 14. vW,,,—v W,, as a function of œw. From Ref. [7]. Typical Regge fits from Ref. [9] 


scaling variable but introduces the good a? dependence of the residue 
function ß>(g?+a?) ' that gives reasonably well the threshold behaviour. 
However local duality is missing in this model. 

It will be fascinating to see, if the marriage of duality and scaling will 
produce a theory with some success. 

Both concepts came from high v and q* speculations, but seem to 
apply to embarrassingly low values of these variables. We hope that this 
uncertain situation will stimulate further theoretical work. 


Postscriptum, November 1971 


Since the lectures were delivered, the work reported in Ref. [17] has been finished. It 
lends strong support to wy as a scaling variable. The parameters are slightly changed 
and the solutions though still not conclusive seem to favor slightly larger a and m. 
Also, at Cornell, the ratio of structure functions for protons and neutrons is reported 
to be even smaller than the one predicted by (25). This is no blow to duality since the 
additional assumption, as explained in the text, of form factor scaling was used. In fact 
this assumption is also known to fail. The situation with the difference in structure functions 
is unfortunately still confused after the Cornell Conference and awaits further experiments. 
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1. Introduction 


Operator product expansions in the vicinity of the light cone [1] have 
been proven to be a useful tool to study processes involving high virtual 
mass photons. However, with the exception of the well known case of 
deep inelastic electron photon scattering, supplementary assumptions 
are necessary. 

In the present lecture we will show that further tests of these ideas 
can be made for the electroproduction of a final state as well. Here too 
some supplementary assumptions will be used. We consider that the 
general structure obtained in the Björken limit for these processes will 
provide a stringent test of the validity of the operator expansions. 

Our result is that for the one particle electroproduction 


“”+A>B+C (1) 


the invariant amplitudes A; can be expressed in the Björken limit 
(v-+ 00, q?-+ — 00; w and t fixed) as: 


Ay, g’, t) as v“ f(a, t) (2) 


where t=(gq- p}, o=2M,v/—q*, v=p4q/M,, q is the photon four 
momentum. The fact that the d; are independent of œw and t is a test of 
the conjecture that states that the singularities of the operator expansions 
near the light cone are c-numbers. As will be discussed below the power 
of our method may be increased if we conjecture that the numbers d; 


* Lecture delivered at the International Summer Institute in Theoretical Physics, 
DESY, July 12-24, 1971. 
** Work supported in part by Deutscher Akademischer Austausch-Dienst grant. 
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may be obtained assuming naive dimensionality for the operator 
products. The same type of analysis can be applied for many particles 
electroproduction: 

“yw” + A>B+C,+C,+—-4+C, (3) 


if the particles C, are fragments of A and B is a fragment of “y” and 
henceforth to the inclusive process 


“y” + A> B+ anything 


if B is a fragment of “y” and we are in the fixed missing mass Björken 
limit. 


2. Application of Light Cone Expansion for one Particle 
Electroproduction 


We illustrate the method by studying the scattering of scalar photons 
on scalar equal mass isospinless particles 


“y” (a) + A(p*) > B”) + C°). (4) 


Straightford application of the reduction formalism allows us to write 
the amplitude as 


T ~ J dxe'**0(— xo) ALL), ja(0)]|C> (5) 


where j(x) is the “electromagnetic current” and jg is the source of the 
scalar field. We have omitted possible terms containing equal time 
commutators, a procedure that we justify later. The approximation 
scheme of the expansion demand that one finds a frame in which light 
cone dominance can be proven. We choose the rest frame of particle A. 
In this frame: 
4=(4o,0,0,43); p” = (po; 0, p2, p3) 
p“=(M,0,0,0);  p°=(p5;0, p$, pS)- 


We consider the limit v>oo and œ, pS and p§ fixed (particle C is a 
fragment of A). This corresponds to fixed t limit since 


|p°|=(1/2M) [tt -4M?)]*, 
cos 0c = (2M? — wt)/2Mo|p°| (p$ = |p| cos 8o) . 


(6) 


(7) 


On invariance grounds, the matrix element <AILi(x), jp(0)][C> can 
depend on x?, xp“, xp© and t which are all independent of v. Hence all 
the v dependence is isolated is the exponential e'4* since qo =v, 
g32v+M/o. 
Now 
expiqx = expiv(xy) —X3)+iMx,/o. 
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The integral picks up contributions from |xo—x3| + 1/v and 
[xo + x3|<@/M. Hence |x5 — x3|< w/Mv. From causality we know that 
the matrix element is non zero only for (x3 — x3) — (x? + x3) 20 so that 
we have finally x? @/Mv-—0 and hence light cone dominance. We 
now apply the light cone expansion of operator products: 


8(— Xo) Li), jpO)] = C(x) F(x) ++ (8) 
where 
C(x) = 0(— xo) |(— x? + iexg)* — (— x? — iexg)4 (9) 


and F(x) is an operator which is nonsingular for x? =0. Inserting (8) 
in (5) we get: f 

T= | dxe'*C(x) F(xp4, xp“, t) (10) 
where F(xp4, xp“, t) is the matrix element of the nonsingular operator F. 
In order to do the integrations in x we consider the double Fourier 
transform of the matrix element: 


Fxpxp = | f dads g(a, dexpilap*+pp)x. (11) 


The limits of integration in (11) are obtained from the spectral conditions. 
(In writing (11) we have assumed that there are no end point singularities.) 
We use now the identity 


f dxe” C(x) = C(d) h(k?, ko) (12) 
where 
C(d)=n2”"*’T(d+1)T(d+2) 
h(a, b) =e" *** f(a +i tT? —(—a— ie) 4-7] (13) 
+ (a+ ie)~4~? — (a — ie) t? 
and get 
1 © 
T~C(d) J $ dadß g(a, B, t)h(E?, čo) (14) 
where 


Šo = (qt ap* + BP) 
é =(qt+ap* + Bp)’. 
This is all what light cone expansions give us. In order to prove scaling 
we shall now assume that only finite « and ß contribute to the integrals in 
(14). In this case 
fo =v 
éE ~2Mv[a+ B(L — 1/w)—1/o]. 
A similar assumption was not necessary in proving scaling from light 
cone dominance for the deep inelastic total electroproduction cross 


(15) 
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section. However, it can be shown that at least in some models our 
assumptions are valid [4]. In fact the same type of assumption is needed 
in order to prove light cone dominance in mass dispersion relations for 
form factors [5]. 

From (14) and (15) we get finally 


T~? flow, t) (16) 


where d is independent of œw and t and depends only on the asymptotic 
dimension of the commutator of the two considered currents (one 
electromagnetic and the other one hadronic). The expression (16) gives 
the scaling behaviour for one particle electroproduction in the fixed t, 
Björken limit. 

The fact that in the Björken limit the real and the imaginary part of 
the amplitude behave in the same way proves that there is no need for 
contact terms and this justify the omission of terms containing equal 
time commutators in the reduction formula (5). 

An application of these ideas to pion-electroproduction may be 
found in Ref. [2]. 

It is very important to mention that the scaling behaviour (16) 
was also obtained in the ladder approximation [6] of the ° theory 
using the method of Altarelli and Rubinstein [7]. Since this model admits 
a light cone expansion [4] this leads support to our form (16). 

A different behaviour of the amplitude in the Björken limit is obtained 
in a dual resonance model for currents [8]: 


T ~ v* f (ca, t) (17) 


where č(t) (independent of œ!) is related to the Regge trajectories included 
in the model. May be that there is some incompatibility between light 
cone expansions and usual dual models but at the present stage it is 
very difficult to decide since on one side supplementary assumptions 
were necessary in our derivation and on the other side dual models for 
currents have still their own problems. It is worthwhile to mention 
Müller and Rühl’s conjecture which also leads to the Regge type 
behaviour (17). 


3. Generalization for many Particles Production 


The same procedure may now be applied for the process 
“P+ A>B+C,+--4+C, 


for fixed t=(q—p®)* and the particles C,,...,C, being fragments of 
particle A. Now 


Thai {dx e'4*6(— xo) <A), ia(O)]|Cy, wig are (18) 
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Considering again the rest frame of particle A. the particles C,,..., C, 
having finite momenta, the matrix element will not depend on v and 
again the whole v dependence is isolated in the exponential. We get 


Tari ~v? fla, t, PP) (19) 


where d is the same like in the one particle electroproduction due to 
the fact that d is related only to the light cone properties of the commutator 
of the same currents. 

The fact that the exponent d is independent of how many particles 
are electroproduced, enables us to consider the inclusive process 


“y”+A>B + anything 
The differential cross section of this process is a function of 4 invariants: 


poda/d p” = f (v, °, M’, t) (20) 
where 
M’ =(q+p*—p®)?; v=qpM. 
We now consider the Björken limit in the fixed missing mass limit. 
We shall prove that in the limit v> 00,9? > — œ, M’, œ, t fixed we have 
SO, g, M’, =v" fo, M’, 1). (21) 
In fact 


d Ci d Cn T 2 n i 
phdojdp?~ Dalaran 5%) e 
n Po Po ~ i=1 


where 


A= s? + M*+q*+2M7q? +2sq*—-2M’s 


but in the fixed missing mass limit we have (we are again in the rest 
frame of particle A): 


n 2 n 2 
(È p°) =M? [pi~ x p°) =t (23) 
i=1 i=1 
and hence 
u : 1 
Y Po = Mr) (24) 
i=1 


which means that all the 3 moments of the particles C, remain finite 
in the Björken limit. Thus we can apply our standard technique to get 
the asymptotic behaviour (19) and from (22) we finally obtain the scaling 
rule (21) with d = —2d — 5. 

All these derivations have included a first step in which the right 
reference frame and the proper reduction formula has been used in 
order to get the light cone dominance. However light cone expansions 
are just in one variable (x?) and can be used just to do one integration. 
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This was enough for the forward Compton amplitude in order to prove 
scaling but this is not the case in general. Thus the next step consists in 
either making models for the remaining unknown function or, as we 
did, to make some general assumptions which allowed us to prove 
scaling. As concerning the first step let us remark that instead of using 
the reduction formula (5) we could use 


T~ [ax e'"*0(xo) CAIDO, j0] C (25) 


But in this case pz = M4 is finite so that the proof does not hold. In 
fact the matrix element is v dependent and we are unable in this way to 
“see” the light cone. 

We have been able to prove scaling like behaviour only for the differ- 
ential cross section of exclusive processes in certain kinematical regions. 
This does not mean that the integrated cross sections should scale and 
a priori we have room essentially for any kind of behaviour for the 
integrated cross sections. Although from some models scaling has been 
conjectured also for these quantities [11], such models seem to have no 
obvious connection with light cone dominance. 

During the preparation of this talk I have enjoyed discussions with 
R. A. Brandt, G. Preparata, B. Renner and H. R. Rubinstein. 
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Duality and the Pion Electromagnetic Form Factor 


KERSON HUANG 


From general principle the pion electromagnetic vertex, as illustrated 
in Fig. 1, should have the form e(p,—p,),F,(t), where t=(p, — p}. 
The factor (p2 — pı), insures gauge invariance when the two pions are 
the mass shell. The invariant function F,(¢) is called the pion electro- 
magnetic form factor. 


P2 


t = ( P2 =, Pı ) 
T 
Pi 
Fig. 1. Pion electromagnetic vertex 


The normalization condition F,(0)=1 defines the pion charge 
to be e. 

We have only limited experimental information about the pion 
form factor. The pion electromagnetic radius has been cited to be [1] 


r, = {6 [d F, (9/dt],-0}? =0:86 +0: 14x 10713 cm. (1) 


Colliding e”e” beam experiments [2,3] have measured F,(t) in the 
neighborhood of the g resonance t=m2. When fitted with a Breit- 
Wigner formular with constant width, these measurements yield a ọ 
width of about 100 MeV, which is noticeably different from the width 
of about 130 MeV that has been obtained from two-pion production in 
purely hadronic reactions. These measurements, however, are not 
sufficiently accurate to determine the detailed shape of the ọ resonance. 
Indeed they can also be fitted by a squared Breit-Wigner formula [4]. 
Data on the large |t| behavior, either for t>0 or for t <0, are lacking. 
On the theoretical side the situation is equally unsatisfactory. The 
simplest phenomenological model based on the assumption of vector 
dominance would give 
F,(t)=C(t—mo+im,I,)*. (2) 


One might add to this contributions from other vector mesons. However, 
the t dependence of C and T,, which became crucial as one goes away 
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from the resonance position, are unknown. Within the assumption of 
vector dominance, a theory of F,(¢) simply becomes a theory of vector- 
meson propagators. 

If,one adopts the view of Lee and Wick [5], namely that propagators 
are regularized by unstable particles of negative metric, then one could 
arrive at 

F,,(t) = | dm? - A(m?) (t — m°)! 


+ B[(t- M? + iMy)"! +(t- M?- iMy)"'], 8) 


a form explicitly proposed by Schwinger [6] in his theory of sources, 
in which B is such as to make F,(t) vanish faster than t~! as t>w. 
As noted by Schwinger, the unconventional last term in (3) is permissible 
because his theory “does not have to contend with the restrictions 
imposed by an assumed underlying operator field structure, nor with 
preconceptions about analyticity”. The Lee-Wick theory represents a 
definite way in which the usual “preconceptions” are dispensed with. 
In the theory of Kroll, Lee and Zumino [7] vector dominance is given 
a more dynamical role via the concept of “field-current identity”. They 


find F,(t) = F,()/(t—m,)’, (4) 


but the “g form factor” F, (t)is not calculable within that theory. Although 
there is no definitive answer about the form of the g pole from these 
theories, we see that, in the form factor, the ọ pole may have a more 
complicated structure than a simple Breit-Wigner formula. 

None of the above theories take into account hadron dynamics, 
or what little we know about it. On the other hand, there have been 
attempts to construct hadronic currents by emphasizing the duality 
aspect of hadron dynamics, and then try to satisfy gauge invariance and 
current algebra requirements [8]. These have not been successful so 
far, and I shall not go into them. 

I would like to explore here a very crude model [9] that combines 
duality and current algebra requirements in a simple way. Our object 
lies not so much in producing a formula that can be compared directly 
with experiment, but rather in exploring the type of qualitative behaviors 
that could possibly emerge. In this sense it is an experiment in a “theoreti- 
cal laboratory”. 

If the pion had no strong interaction, its electromagnetic vertex 
would be e(p, — pı),, as represented by the first Feynman graph in Fig. 2, 
and the form factor would be unity. The strong interactions give rise 
to final state interactions, represented by the second Feynman diagram 
of Fig. 2, in which A is the pion-pion elastic scattering amplitude. In 
principle one may have other graphs, in which the two — pion inter- 
mediate state is replaced by other hadronic states, namely states con- 
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taining any number of bosons and baryon-antibaryon pairs of arbitrary 
spin. However, it does not violate any known principle to make the 
following model: Only pions pairs are coupled to a photon through the 
elementary vertex e(p,—,),- All other charged particles are coupled 


P2 


P, 


Fig. 2. Our model for the pion form factor 


to the photon in a non-elementary way, namely by first becoming a 
pion pair via the strong interactions. This is our model. If philosophy is 
important, we might adopt two different interpretations of the model: 

(1) The model is a mathematical exercise. We really believe that 
other hadrons also couple to the photon in an elementary way, but we 
consider only the pion contribution for the moment. 

(2) We believe that the electromagnetic field single out the pion as 
an “elementary” particle. That is, pions are partons. 

I tend to favor the second, if only for the reason that it is less apolo- 
getic; but it must be admitted that, if taken literally, this view seems to 
have been ruled out by experiments. The SLAC-MIT experiments on 
deep inelastic e — p scattering showed that the virtual photon in that 
process is predominantly transverse, and this would require the partons 
to have spin 4 instead of 0. Thus we have to retreat a little and say that 
we would perhaps like to use spin 4 partons, but for technical reasons 
(namely, that a dual resonance model for boson-fermion scattering is 
unknown) we pretend that they have spin 0. 

The advantage of introducing electromagnetism in this manner is 
that gauge invariance and current algebra requirements can be easily 
satisfied. Gauge invariance is automatic, being a consequence of the 
gauge invariance of the elementary vertex e(p, — pı), Current algebra 
sum rules require [10] that be the residue of a fixed angular momentum 
pole at J =1 in the crossed channel of the Compton amplitude for the 
scattering of a charged photon by a pion. We satisfy this by defining the 
compton amplitude by the Feynman graphs of Fig. 3. It can be shown 
easily that the residue of the fixed pole in question is obtainable from 
the graphs of Fig.3 by “pinching” together points 1 and 2, thereby 
yielding the form factor. 

Our plan is then to take the Feynman graphs of Fig. 2 and substitute 
for A a suitable dual amplitude. It should be noted that our approach 
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is different from one in which one first tries to calculate approximately 
ImF,(t) and then use a dispersion relation to obtain F,(t). The latter 
approach would require a knowledge of the amplitude for 27 (arbitrary 
hadrons), and we want to avoid having to make models for that. In our 
approach the effect of the many hadron intermediate states are swept 
under A. 


Fig. 3. Our model for pion Compton scattering 


We decompose the pion amplitude A into the contributions: 
A = Apom + A, > (5) 


where Ap,,, contains the Pomeron trajectory and has no resonance 
poles, while A, is approximated by the Veneziano amplitude [11] 
which contains only resonance poles, the lowest one being the g. Since 
the t-channel has isospin 1, we have [12] 


A,(t,s)= — PV(t, s)— V(t, u)], (6) 
V(t, =r =x) P(l—a,)/P(l—a,—4,), (7) 
a =a+bt. (8) 


The trajectory «, is completely determined by the conditions 


a(ms) =1, 


(9) 


alm) =3. 


The first states that the g has spin 1, and the second is required to satisfy 
the Adler condition [12]. They lead to 


T. 


b =4(m} -m})"' = 1 (GeV/c)? . 


BEE: GER ER Big? RE a OE I 
a=7 zm, (m, m Xa, 


(10) 


The constant £ in (6) can be related to g meson parameters by calculating 
the decay rate g > 2 within the Veneziano model, leading to the relation 


B = 16n(6T,/m,), (11) 
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where I, is the ọ decay width. Noting the experimental relation 
6I,/m, © 1, we obtain 
Bx 16x. (12) 


Thus there are no unknown parameters in A. 

An important assumption is made here, namely that the amplitude 
A, is the same as that on the mass shell. The two intermediate pions, 
of course, go off the mass shell in the Feynman graph of Fig. 2. How the 
amplitude should depend on their masses is an open question. In field 
theories, such dependences are determined by the way current operators 


P2 
(1) (p2-Py Sp = —~<OL 
Py 
o OY DIE 


aM 


Fig. 4. Calculation of the form factor 


are constructed out of field operators, in accordance with current 
conservation and assumptions of locality. We do not have such princi- 
ples here. Therefore we make the simplest assumption, albeit a strong 
one, that there is no mass dependence. The consequence of this assump- 
tion will be clear, as will be the consequence of modifying it in other 
simple ways. 
Corresponding to the decomposition (5), the form factor can be 
written as 
F(t) = 1+ Gpom(t) + Go (t), (13) 


and we shall calculate G,(t). To sketch the calculation we refer to the 
graphs of Fig. 4, where on the first line G,(¢) is represented as an Feyn- 
man integral involving a Veneziano amplitude V. On the second line 
we note that the Veneziano amplitude is dual, that is, it can be represented 
as a sum of poles, either in the s or in the t channel. Thus, we see on line 
three that G,(¢) can be expressed as a sum of triangular graphs. Each 
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triangular graph can be evaluated using standard Feynman graph 
techniques, and the answer turns out to be an integral whose integrand 
is of the form (t + x,)"!, where x, is some complicated function indepen- 
dent of t. The important point is that it is again of the form of a pole 
in t. Summing over these pole terms again gives back a Veneziano 
amplitude. In this manner we obtain 


G, (t) = (8/47? t) ie V(t, — x) g(t, x), (14) 
ò 


where the function g(t, x) is explicit given in Ref. [9]. We give here only 
the simple form it approaches when the pion mass is neglected: 
g(t, x) — (4 -x/)In(l-1/x)-1. (15) 


m0 


The logarithmic factor gives rise to a branch cut of G,(t) from t=0 
to œ (4m2 to œ, if m, had been kept finite). 


(Pd Gplt) = -o< 


Fig. 5. Why second order poles occur 


This elastic cut, of course, arises just from the Feynman propagator 
used for the two intermediate pions. Any inelastic cut that G, (t) might 
have can arise only from the cuts of V(t, —x), which in the Veneziano 
model are approximated by an infinite distribution of poles. 

From (14) and (7), we see that G,(t) is proportional to F (1 —«,), 
which has simple poles at ,=1,2,3,... corresponding to an infinite 
number of vector mesons 0,0',0",..., with equally spaced squared 
masses. When one attempts to calculate the residues of these poles from 
(14) by setting a, = integer, one finds that they diverge. If one first 
calculates the residue function for «, +n, and then approaches «, =n, 
one finds that the residue function itself has a simple pole at œ, =n. 
Therefore near «, =n, G,(t) has a second order pole, plus a simple pole. 
This can be seen quite directly by noting that by duality we could have 
formally replaced line 3 of Fig. 4 by the graphs shown in Fig. 5. We see 
that each graph separately diverges because of the self-energy loop, 
and the divergence is directly traced to the assumption that A has no 
mass dependence when extrapolated off the mass shell. 

It is not difficult to show that if the amplitude is damped by a factor 
(m?) "" for each external mass m off the mass shell, then the first n vector 
mesons g,0',0”,... will appear as simple poles, but the (n+1) st and 
higher one will continue to be second order poles. An exponential mass 
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damping will make them all simple poles. Since there is no a priori 
reason why vector mesons should appear as simple poles in the form 
factor, we stick to the original version of the model, which is free of 
arbitrary parameters, and certainly more intriguing. 

Let us examine G,(f) in more detail near the vector meson poles, 
for in these neighborhoods F,(t) = G,(t). Thus, at the g and g’ poles we 
have 

F, (0) © (6/4827) (1 — a) `+ b7? (t — m)? 


16 
+i0(t-m})!, (tmz), 19) 


F, (t) © (8/4827) (2 — a) Ba — 1)! b7? (t-m)? 


17 
+i0(t—mZ)*, (tam), en 


where bm? ~}, bm’? 3. Note that the second order poles are real, 
while the simple poles are pure imaginary. For a = 4 the residue of the 
g’ second order pole is $ that of the g. 

These poles, however, are on the real t-axis, because the Veneziano 
model isa zero-width approximation. How might these results be changed 
if widths are taken into account? We can, of course, only offer a guess, 
for we do not know the correct way to unitarize the Veneziano model. 
Since we know experimentally that 92x is the overwhelmingly pre- 
dominant decay mode of the g, we may assume that near the g pole it 
is a good approximation to assume elastic unitarity, which requires 


Im F,() = Fr() TO), (18) 
where T.(t) is the n-z p-wave, I =1 scattering amplitude at squared 
energy t. 

Elastic unitarity again requires 


Im T, = T Ty. (19) 
Therefore the general solution to (18) is 
FQ=ROT., (20) 


where R(t) is a real analytic function. The entire content of (20) is that 
the phase of F(t) is the same as the phase of T(t), i.e. the p-wave, I=1 
phase shift for z-z scattering. Let us assume that 7,(t) contains the g as 
a simple pole. There are then two possibilities: (a) The second order pole 
in F,(t) arises from a confluence of two poles of 7;(t) in the limit of 
zero width; (b) It arises from a confluence of the ọ pole in T; (t) with 
a pole in R(t), in the limit of zero width. The first possibility can be 
ruled out as follows. 

It can be shown that if 7, contains two simple poles close together, 
then (19) requires that their contribution to the phase shifts add, and this 


Duality and the Pion Electromagnetic Form Factor 105 


leads to 
PAREPA D; (21) 


T(t) = Alm til), (n= 1,2). 
In the limit m, >m, >m, and A, >4,-0, we would have, by (20), 
F,(t)® RL (A, + Az) — m?) + 21d, Aa(t—m’)’], (22) 


which disagrees with (16), because the second order pole is pure imagi- 
nary instead of real. Therefore we can only have the second possibility, 
and since R(t) is real analytic, its poles must occur in complex conjugate 
pairs. This leads to 


F(t) —A,(t—mz +id,)-'[C+D((t—M? +o)! +(t-M?-ie)"')] (23) 


which approaches the form (16) when m>M and 4,>¢-0. It is inter- 
esting that (23) is of similar form to (3), if m, and M are further apart 
than the widths. 

To make some rough comparisons with experiments without having 
to adjust the numerous constants in (23), we may just take (16) and 
replace m, by m,—iI,,/2, where 6I,/m,~1. We then obtain, near 


t=m?, = = 
me |F. (|? = (6r b?) [t — m)? + (mT) 1-7. 24) 
In particular, at the g peak 
|F,(m2)|? = (62) *(bm, 1) * = 30, (25) 


which lies within experimental accuracy [2, 3]. If we simulate (24) with 
an equivalent simple Breit-Wigner formula of the same height and width, 
we obtain 
[Fase Ol? © (67)? (m, Tare) (b mT) it- m3) +(m, Dar) ] ~! 
Tas = (2 — 1P T, = 0.657, 


which could explain the apparently smaller width of the ọ in colliding 
beam experiments. 

Other properties of G, (t) would be less directly related to observations, 
because of the presence of the term 1+ Gp,,,(t). The contribution of 
G,(t) to r„, ignoring Gp, (t), is divergent as m, +0. The leading divergent 
term gives 

[6 G,(0)]? =p [12b In (1/bm})]? = 1.4 x 1072: em, (27) 


(26) 


which is of the same order of magnitude as the experimental value (1). 
At t=0 we find that G,(¢) vanishes in the limit m, —0: 


G,(0) -p — Bbm?/8r. (28) 


This fact is a direct consequence of the Adler condition that A, satisfies. 


Hence we expect Gp,,,(0) to behave similarly. To order bm2, therefore, 
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the normalization condition F,(0)=1 is automatically satisfied. The 
cancellation G,(0) + Gpom(0) = 0 need be enforced only when we “turn 
on” the pion mass. 

The asymptotic form of F,(t) is of particular interest, but it can only 
be a subject of conjecture for us. We do know by direct calculation 
from (14) that 


G,(t) — (B/8n2) F(1—a)(—bt) + 0(1Anbt)]. 29) 


Assuming that Gpom(t) behaves similarly except we set a=1, ignoring 
the divergence of F(1— a), we would conclude that 


— 
|t] +0 


Grom (t) I isa? constant, (30) 


and it is very templing to conjecture that this constant should be exactly 
— 1, thus cancelling the “contact term”. This would lead to 


FD gpp t+. (31) 


It is difficult to construct a model similar to ours for the nucleon 
form factor, for the Veneziano amplitude involving fermion is unknown. 
However, if one writes the usual Veneziano amplitudes for the invariant 
amplitudes (thus ignoring difficulties with the occurence of parity 
doublet), one obtains in place of (31) 


Fyucteon Okres lt] >a nr? (32) 


when «y(0) is the intercept of the nucleon trajectory, taken to be —4. 
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Deep Inelastic Hadronic Scattering 
in Dual-Resonance Model 


KERSON HUANG 


There has been much interest recently in the so-called inclusive 
hadronic reactions, namely a reaction of the type 


A+B-X + anything (1) 


where A, B, and X are hadrons. Most of the works in the literature (1) have 
dealt with special regions near the phase space boundaries of the reaction, 
in which the general behavior of the cross section can be predicted by a 
new type of Regge phenomenology, without having to assume a detailed 
dynamical model. In this respect they are similar to the Regge phe- 
nomenology of two-body near-forward scattering. 


backward fixed angle forward 
scattering scattering scattering 
(fixed u) N (fixed t) 
See A 
A ai $< 
1 0 +4 Z 


Fig. 1. Phase space for two-body reaction at given s 


I would like to explore the behavior of such inclusive reactions in the 
phase space far from the boundaries — a region that I shall refer to as 
the “deep inelastic region” in analogy to the well-known case of electron 
scattering. Here one can make predictions only if one makes detailed 
dynamicalassumptions. In this respect it is similar to wide-angle scattering 
in two-body processes. Before we go into the subject, it is helpful, for 
orientation, to review those aspects of two-body processes that are 
relevant to our developments. 

For definiteness let us consider p— p elastic scattering at high 
energies. The kinematic is completely specified by two numbers: the 
squared CM energy s, and the squared 4-momentum transfer 
t= —2p*cos@, where p is the magnitude of the 4-momentum, and 9 
the scattering angle, all in the CM frame. For fixed s, therefore, the phase 
space @ may be represented by the line segment —1 <z <+ 1, where 
z= cos, as shown in Fig. 1. The shaded positions near the phase space 
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boundaries are those amenable to conventional phenomenological Regge 
analysis, with the differential cross section do/dQ behaving as f(t) s’®, 
or f(u) s*™, respectively. The extensions of these shaded portions shrink 
to zero as soo, In the wide-angle region, which is almost all of phase 
space as s—> 0, we need dynamical models, of which there are few, and 
none so far is completely satisfactory. The experimental data in the 
10 — 30 GeV/c range of incident laboratory momentum seem to show a 
remarkably simple scaling behavior. Orear [2] suggested the empirical 
formula i 
cals, 2)=do/dR = Const. e275"? , (2) 
where a is approximately independent of the energy, and has a value of 
about (160 MeV/c)"!. Krisch [3] suggested that a more accurate fit 
can be obtained by replacing psin@ by (psin@)*. On the basis of a high 
statistics experiment, the CERN group [4] found that a still better 
empirical formula is 
oals, z) X Const. e7 4558 , (3) 


None of these really fit the data in all the fine structure, but they reproduce 
the gross features rather well. 

On the theoretical side, there is only one dynamical model that is 
reasonably precise, based on premises subject to other independent 
tests such as crossing symmetry and duality, and from which predictions 
can be easily obtained, and that is the Veneziano model [5]. It predicts 


a(S, ) 5 CO, 
z fixed (4) 
g(z)=(1 +z) In5(1+2)+(1—z) Ind(i—z), 


where the quantity C generally depends on s and z, and also depends on 
specific choices of terms in the Veneziano model (e.g. st, su, or tu terms); 
but for large s, it varies at most like a power of s, and can be regarded as 
a constant compared to the exponential factor. A very good numerical 
approximation to the function g(z), for the whole range of z=cosß, is 


g(z) x —(2In 2) (sin 0)? . (5) 


Thus by comparing (4) and (5) to (2) and (3), we see that the Veneziano 
model gives a prediction that fits well the gross features of experimental 
data, provided we regard b as an adjustable parameter. 

Now within the Veneziano model b is fixed, being the universal slope 
of nonvacuum Regge trajectories, and should havea value b x 1 (GeV/c). 
On the other hand to fit the data one needs b= 4(GeV/c) 7. One might 
“explain” this by saying that the Veneziano model leaves out the 
Pomeranchuk trajectory, which is expected to have a smaller slope, and 
that if one known how to incorporate it into the Veneziano model, 
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one might get (4) with a smaller effective value of b. We therefore do not 
really have a prediction from a theory, but merely a suggestive result 
from an incomplete theory, the improvement of which, of course, has 
been the major concern of many theorists. 

In the same spirit as the application of the Veneziano model to 
elastic widge-angle scattering, let us apply the dual-resonance model — 
the generalized N-point Veneziano amplitude — to the inclusive hadronic 
reaction (1). The kinematic is now specified by two numbers, apart 
from the squared CM energy s. Many experimentalists use q, and qyp 

i U fixed U fixed 
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Fig. 2. Phase space for inclusive reaction at given s. Deep inelastic region is entire rectangle 

except for the shaded strip, whose width decreases as s increases. Lines of constant qy/p=r2, 

constant q,/p=r(1—2z’)? are also shown 
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the transverse and longitudinal momentum of the produced particle X 
in the CM frame. As we shall see, however, it is more convenient to use 
the CM momentum q of the produced particle, and its CM angle 6 
with respect to the incident direction. We shall introduce 


r=q/p, (OSrs1) 


z=cos0, (-1sz<SI1), 


(6) 


where p is the incident momentum in the CM frame. Thus, for given s, 
the two-dimensional phase space is represented by the r—z plane in 
Fig. 2, on which is also shown lines of constant q,/p and q,/p. The 
shaded strip near the phase space boundaries lzi ={1 and r =Q are regions 
analogous to the forward and backward scattering regions in two-body 
processes, in that the behavior of the cross section can be deduced from 
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the new Regge analysis of A. H. Mueller (Ref. [1]). Other conventionally 
used kinematical variables are illustrated in Fig. 3, with definitions and 
asymptotic forms in terms of r and z given by 
s=(p,+ Pa 4p’, 
t=(pa— px)’ > -3sr(l-2), 
u= (pg — px) > -3sr(l+2), 
M?’ =(pa+pg— px =s+t+u-m2 mi —m} 


(7) 


>s(i—r), 


where “>” denotes the limit s— oo with r and z held fixed. 


X INVARIANT 
MASS=M 
t u 
s 
A B 


Fig. 3. Kinematical variables for inclusive reaction 


The region near z = 1 is called the fragmentation region of the incident 
particle A, because as s— o0, the produced particle X carries off a finite 
fraction r of the incident momentum, with however a finite momentum 
transfer t from the incident particle. Thus it came off as a “fragment” of 
A. The differential cross section in this region is given by 


d’o/dgd2 z[(1+r)/1-r)®, (fragmentation region) (8) 


where a(t) is an appropriate Regge trajectory. Similarly the region near 
z=-—1 is the fragmentation region of the target B. The region near 
r=0 is called the pionization region, from a terminology originating 
in cosmic-ray experiments, because even as s> 0, the produced particle 
X has finite momentum q. Mueller’s analysis predicts the cross section 
to be a function of q, =qsin@ alone, but is unable to determine the 
functional form. DeTar et al. (Ref. [1]) calculated it explicitly from the 
dual-resonance model: 


d?c/dqdQ x~exp(—4bq’sin?@), (pionization region) (9) 
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where b is the universal slope of linear Regge trajectories in that model. 
This yields a transverse momentum cut off of 500 MeV/c, which is in 
qualitative agreement with the experimental value of 300 — 500 MeV/c. 

We are interested in the region complement to the regions above 
(i.e. the region outside of the shaded strip in Fig. 2), and we shall call it the 
deep inelastic region, in analogy with ep inelastic scattering. For 


-1 0 +1 Z 
Fig. 4. Lines of constant œ = — 2myv/t for inclusive reaction 


comparison we give the conventionally used auxiliary variables in the 
latter case: 
q =(Px— pa} =t, 
V= Pg’ (Pa — Px)/Mmg = (t — u + mg — m3)/2mp , (10) 
= — 2myv/g?. 


The deep inelastic region in ep scattering is defined as that in which 
vom at fixed w: 


owxl—M2/t~1+21—n/r(l—2). (11) 


On our r—z plot, therefore, the lines of constant œ are as indicated in 
Fig. 4. 

The deep inelastic region here is mathematically defined by the rule 
that we take the limit s> first, with r+0,1 and |z|+1 fixed. Only 
after this is done shall we approach the phase space boundaries. In par- 
ticular the quasielastic region r=1 means s—oco,r—1, in that order, 
so that even there the final state of the reaction containts many particles. 
As s— oo, the inelastic regions cover almost all of phase space; but of 
course the cross section is vanishingly small, most of it being concen- 
trated in the fragmentation and pionization regions. 

To calculate the inclusive cross section d?o/dqd@ from the dual- 
resonance model, we proceed as did Mueller and DeTar et al. by first 
making use of a generalized optical theorem, which relates d?o/dgqd@ 
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to the absorptive part of a 6-point function in the variable M?: 
d?a/dqdQ=ImB,(M’), (12) 
which may be illustrated by the graphs in Fig. 5. 


UE M2) 
A x X A 
- 2 M2 
u 5 (R2- M?) mem 
X X 
_A A 


M? 
Fig. 5. Illustration of the generalized optical theorem [Eq. (12)] 


is. 


The dual-resonance model is then used to calculate Im B,. Thus (12) 
is an exact statement implied by unitarity, but we approximate the right 
hand side by using the dual-resonance model, which does not satisfy 
unitarity. The dual-resonance model for Bg is of the form 


B,(M’)=n"! 2 R(n)/(n— M°), (13) 


where dependence on variables other than M? has been suppressed. 
Therefore, strictly speaking 


ImB,(M2)= Ș R(n) d(n—M?). (14) 
n=1 


For M? =œ, however, we may perform a local average over M?, and 
obtain 
Im B,(M?) = R(M’). (15) 


This residue function is easily calculated, and has been given by DeTar 
et al. (Ref. [1]). We only have to take their result and evaluate it explicitly 
in the deep inelastic region. The result [6] turns out to be extremely 
simple and suggestive: 

d?a/dqdQ a> a0", (16) 

r.z fixed 
G(r, z)=(1+ rz) Ing(it+rz)+(1—rz) In$(1—rz) 
-(+r)In41+r)-(1-r)In$(1-r), (17) 
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where b is the universal slope of linear Regge trajectories in the model: 
a(t)=a+bt, (18) 


and o, is a function of s,r, z, that varies at most like a power of s as 
soo. Its explicit form depends on the details of the reactions (i.e. 
quantum numbers of the participants). The function G(r, z), however, 
is universal. For sufficient large s, ¢9 becomes unimportant, for according 
to (16) 

F=s ‘d*o/dqgdQ—bG(r, z) + O(Ins/s). (19) 


For comparison with existing data, however, we take o, to bea constant: 
Co =1 mb/(Sr — GeV/c). (20) 


a value consistent with the observed p — p total cross section. By com- 
paring (17) with (4) we obtain the interesting formula 


god? a/dqgdQ a> Ca lS, 72)/oals, r). (21) 
r,z fixed 
That is, in the dual-resonance model, the inclusive cross section in the 
deep inelastic region is determined by the cross section for wide-angle 
elastic scattering. 
The quasi-elastic limit r—>1 is particularly interesting, for (21) gives 


d?a/dqdQ = Const. oalr, 2). (22) 
r>1 


We note from (7) that in this limit t> M?, although M?— oo. Thus even 
though many particles are produced, the cross section is like that for 
elastic scattering, as long as the momentum transfer is much larger than 
M?. 

We emphasize that (16) and (17) apply only to the deep inelastic region. 
They do, however, approach limits that are qualitatively correct when one 
approaches the phase space boundaries |z|=1 or r=0. For example, 
as r—0, one recovers the pionization limit: 


In(d?o/dgdQ)> —4bq?sin?6. (23) 
As |z|>1, one recovers the fragmentation limit: 
In(d2s/dgd@Q)>btln[(i+r)/A nl. (24) 


which differs from the correct expression (8) only in that bt should have 
been a+ bt, where a is the Regge intercept. 
We now compare our model with data by examining the ratio 


b=s7'In(d?a/dq dQ)cy/G(r, 2), (25) 


114 Kerson Huang: 


when (d*a/dqdQ) y is the experimental cross section in the CM frame, 
in mb/(Sr-GeV/c). Scaling is tested by seeing whether b is independent 
of s. The dual-resonance form (17) of G(r, z) is tested by the additional 
requirement that b be a constant, whose value would then be the slope 
parameter of some effective Regge trajectory. 

We use the data at 30 and 20 GeV/c of Anderson et al. [7] and those at 
12.5 GeV/c of Akerlof et al. [8]. To exclude the “strip”, we keep only 
values of |z| <0.7, in all the data. The results are shown in Figs. 6-9, 
which are self-explanatory. The extensive data at 19.2 GeV/c of Allaby 
et al. [9], unfortunately fall entirely in the fragmentation and pionization 
regions, and hence cannot be used for our purpose. 

The validity of scaling, and the specific form of G(r, z), seem to be 
borne out. The general feature is that s~'In(d?a/dqdQ) approaches the 
universal function G(r, z) as s>%,2z—0,r—1. For fixed r, the rate of 
approach is faster the larger the s, and the smaller the z. This is illustrated 
most clearly in the n* data in Fig. 7. As an example of the quality of the 
fit, we compare two data points in the 30 GeV/c proton spectra: 
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Figs. 6-9. Plots of experimental s"!In(d?o/dgdQ)m divided by known function G(r, z), 
versus r for different values of z. The experimental differential cross section (d?o/dgdQ)cm 
is in units of mb/(Sr — GeV/c), and s is in (GeV)?. Data at 30 and 20 GeV/c (BNL 1967) are 
taken from Ref. [7], and those at 12.5 GeV/c (ANL 1970) are taken from Ref. [8]. 
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The effective slope parameter b seems to be different for different 
outgoing particles: 


0.3 for X=p 
0.4 for X=x* 

b= 3 (26) 
0.5 for X=n 


0.5 for X=K*,K”. 


The relevance of these results to the electromagnetic interactions 
of hadronics lies in a comparison with some recent speculations by Chen 
and Harte [12] and by Berman and Jacob [13]. They extended the 
Wu-Yang conjecture [14] 

(27) 


Gep->ep I (Ca)? 
to the inclusive reaction pp—> p + anything, and suggested that in the 
region of phase space for which M?/t remains fixed as r—>1, as indicated 
by the shaded area in Fig. 4, 


d?a/dqdQ~ (oq)? R. (28) 


where R is a slowly varying function proportional to the form factor 
vW, from e+p>e-+ anything. This conjecture is not borne out by our 
result, which resembles more closely (27) with (c,,)? replaced by a... 
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I. Introduction 


These lectures are devoted to a simple current algebra description of low 
energy pion physics. The motivation for a renewed interest in the subject 
is represented by the fact that there is some feeling, at present, to consider 
the small pion mass, rather than a dynamical accident, as the strong hint 
for an approximate SU(2) x SU(2) symmetry of the hadron world. In 
order to establish a common framework let us quickly remind a few 
known facts about currents and current algebra. 

Roughly speaking, one can say that current algebra can be con- 
sidered as the successful attempt to generalize some outstanding results 
of electrodynamics, like the non-renormalization of the electric charge 
by strong interactions and the low energy photon theorems, to the other 
currents used to describe some important properties of hadrons. This is 
achieved by applying to jẹ” a S U(3) rotation and a parity transformation, 
i.e. a SU(3) x SU(3) transformation. As a consequence we have to deal 
with two octets of currents, vector and axial vector, V7 and Aj,a=1... 8. 
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Some of them are used to describe the electromagnetic interactions 
(V3 V) and the weak interactions (V,’”*°, Ab >45) between hadrons 
and leptons. 

The following properties of these currents represent important 


achievements in our understanding of elementary particle physics. 


a) Algebraic Properties 


Consider the weak charges associated with the currents 
Q*= {dx v(x), Q* =f dx Ale). (1.1) 


These objects are assumed to be the generators of the SU(3) x SU(3) 
algebra and to obey therefore the set of equal time commutation relations 


[9° Q*) = ifagyQ’ 
[Q*, OP]=if,,,0° (1.2) 
[0°, OP] = ifapy 0’ š 
The proposed identification represents a non trivial bridge between 
the symmetry properties of the hadron world (i.e. behaviour under 
SU (3) x SU(3)) and the quantities used to describe the lowest order weak 
and electromagnetic interactions of hadrons. Furthermore the validity 


of the above commutation relations is absolute, independent on whether 
the charges are conserved or not. 


b) Conservation Properties 


The Vector Case 


We have 
fon ie =0, 


and 
oH yE? =0 


only in the SU (2) limit, i.e. in the absence of electromagnetism (and weak 
interactions). Moreover 


ge yon =0 


only in the S U (3) limit, i.e. in the absence of “medium strong interactions”. 

These conservation properties reflect well known algebraic sym- 
metries (or normal symmetries). One then obtains, in the symmetry 
limit, familiar results like the classification of hadrons in degenerate 
multiplets corresponding to SU(3) irreducible representations, non- 
renormalization theorems for all charges, relations among vertex func- 
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tions and so on. This realization of SU (3) seems to be the most adherent to 
the real state of affairs. Thus, even if some SU(3) charges are no longer 
conserved in nature and corrections to the symmetry results appear, one 
can still say that there are evident traces, in the physical world, of the 
underlying symmetry. Quantitatively this means that if one writes the 
hadron Hamiltonian as 


H=H,+eH,, (1.3) 


where Hy is SU (3) invariant while H, is not, the breaking parameter e is 
small enough to justify attention to the zeroth and first order predictions. 


The Axial Vector Case 


Let us concentrate on the isospin axial currents and charges A}, 0% 
&=1,2,3. The first hint about the peculiar nature of the conservation 
limit (“chiral” limit) comes from the existence of a longitudinal com- 
ponent of the current, associated with a 0” particle, the pion: 


COL ARINE = ib .p Tufa (1.4) 
<O| HA% my = me fy Sap >» 


showing that the symmetry limit is realized either by m,=0 or f, =0, 
namely either the particle is massless or it is completely decoupled from 
the current. While the second alternative is easily recognized to cor- 
respond to normal symmetries, the peculiar realization by massless and 
spinless bosons characterizes the so-called “spontaneously broken sym- 
metries”, a la Goldstone. A further argument comes from considering the 
nucleon matrix elements 


pa A; Ip) =u,(t,/2) [Ys Ya G4(4°) +4,75 Gp(4°)] u, 
A=p,—py 


and 


(1.6) 
and 
<pal 6" A | py > = iU Ys (1"/2) u [—2 M G4(4°) + 4°G,(4’)]. (N) 
The conservation requirement taken at 4? =0 shows that if 
A? Gp(A*)|42=0 =0 
(as the case when m, +0) then 
MG,(0)=0 (1.8) 


and one is left with consequences of the symmetry, M =0 or G,(0)=0, 
which are both badly violated in Nature. A more appealing result comes 
from the assumption that the axial current is still coupled to a massless 
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pion. In this case Gp(4°) 35 — 2f,9,n/A? and one finds the Goldberger- 


Treiman relation 
Ian” —MG,(0)/f,. (1.9) 


We can consider it as the first example of low energy theorem, ex- 
pressing the nucleon-massless pion vertex in terms of weak interaction 
quantities. The agreement between the above theoretical prediction and 
the experimental value is quite reasonable? (within a 10%). This and other 
successful predictions we shall discuss later, support the idea that there 
is an underlying SU (2) x SU(2) symmetry in Nature, corresponding to 
the conservation of isospin vector and axial vector currents and realized 
by SU(2) multiplets and massless pions”. The immediate consequences of 
chiral symmetry are then of a dynamical nature, since one can derive 
exact low energy theorems for massless pion amplitudes rather than 
algebraic results. 

Actually the pion has a small, but non vanishing, mass (on the hadron 
scale) and chiral symmetry is broken. However, we can say that axial 
currents still show a sort of “partial conservation”, which is concretely 
realized if the divergences D* = 0“ A% (and their matrix elements) are pro- 
portional to some positive power of the pion mass. This is the basic idea 
of the so-called PCAC hypothesis, which, more concretely, reads 
<b| D*|a> ~ O(m?2). Again a unified way of expressing all this is to split 
the total Hamiltonian 

H=H)+eH, (1.10) 


in a part, H, SU(2)x SU(2) invariant plus a part, eH,, responsible of 
chiral breaking: 


O* = | D*dx =i(H, O"] =ielH,, 04]. (1.11) 


Then e is consistently determined to be proportional to m2. 

This scheme supplies the most suited frame to fit all soft pion theorems, 
which can thus be logically understood asa product of the current algebra 
commutators and of PCAC (i.e. approximate SU (2) x SU(2) symmetry). 
Since the experimental verification of these theorems is at hand, the 
problem of evaluating the “corrections” due to the finite pion mass is of 
some importance, both from the practical point of view of having predic- 
tions for physical pions, and with the aim of a better insight into the 
nature of the breaking. 

We shall illustrate in these lectures a simple formalism, based on 
current commutators, to obtain a representation of physical pion ampli- 

1 With frlep = 0.68 m,, G4()lexp = — 1.23, one gets gpwle.7,%12.3 while g.nlexp = 13.5. 

? It has to be noticed that one could look for normal realizations of chiral symmetry, 


resulting in degenerate parity multiplets or vanishing mass nucleons, but the Goldstone 
mechanism seems to be the preferred one. 
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tudes, where both the soft pion limit and the corrections are explicitly 
exhibited. Owing to the chosen realization of the symmetry, the latter 
will vanish with the pion mass and will thus depend both on matrix 
elements and on “higher” commutators involving the axial current 
divergence. The application of the approach will be limited to pion- 
nucleon scattering and pion-electroproduction (with special emphasis on 
the last reaction), namely to “old fashioned” pion physics, where the 
present low-energy description can be considered complementary to dis- 
persion theory or to the static model. Although other phenomena, like 
K,, decay, can be treated with the same techniques, we prefer, for time 
reasons, to concentrate on reactions where there is a sharp distinction 
between the pion and the target and no SU (3) considerations are required. 

Thus we shall only occasionally mention the recent and interesting 
ideas about SU(3) x SU(3) as an approximate symmetry of strong inter- 
actions. This represents a natural generalization of the previous con- 
siderations, which arises combining SU(3) and SU(2)x SU(2). The inter- 
play between these two partial symmetries, the role of kaon PCAC, the 
SU(3) x SU(3) transformation properties of H, in Eq. (1.10), represent 
open and fascinating problems and we hope that the techniques we 
present here can represent a tool for a better understanding of those 
points. 

In Section II we shall review the main features of the formalism and 
treat in some detail the pion-nucleon scattering. Sections IH, IV will be 
devoted to photo- and electroproduction of one pion and toa comparison 
with experiment. 


II. Low Energy Theorems 
a) A Simple Formalism for Low Energy Pion Physics 


One of the remarkable features of current algebra is that several tech- 
niques are available, rather than a single one, to exploit the fundamental 
commutation relations. The various approaches (chiral lagrangians, 
Ward identities, sum rules, etc.) have their own peculiar aspects but, 
after careful translation, the results they give are seen to be more or less 
equivalent. We choose to describe a simple technique of deriving low 
energy pion theorems, based upon the saturation of equal-time com- 
mutators. Of course, the commutators must involve the axial charges 
O* (x= 1, 2, 3) since the symmetry pattern has to be explicit. 
Let us define the operators [1] 


Oi, n=O" + (i/m,) O . 23) 
Os =O" + (1m) O* = Oi, r F (i/m,) Vi. r 
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such that 
<0] OF In’) = — <n? ORIO) = 2 if, m,(2 n)? 6,,6(r) 
<r 0710) = <01 0% In) = <|Q5|0> = <0|O5|2> =0. 
Making use of the relations 

Ot = | dx D(x), D(x) = 6" A(x) 
O5 =m;? [dx(D] + m2) D(x) 
we easily find the matrix elements: 
(27)? f,5(P, — Pa) b| xla) 

i(E, — E,) 


(2 2)? m, f25(Py — Pa) blz la) 
i(E, E E,) (E, + mM, — E,) 


(2.2) 


(2.3) 


<bIO5 la) = 
(2.4) 


<bIÖ5, rla) = 


where we have introduced the virtual pion source y*: 
(O + mz) D*O) = mz fax), 
om) = Kal) =0. 

These simple relations enable us to write down representations for 
pion amplitudes. As an example let us concentrate on the process 
a+ V—n*+b where V can be taken to be an external source, such that 
<0[|V|2z> =<x|V {05 =0, and ja>, |b are the initial and final “target” 
states. The starting point is represented by the introduction of the equal 
time commutators: 


<b| (Qt, V]|a> = <b| E,| a> + (i/m,) <b|X.la) , (2.6) 
E,=[0%,V]; X,=[0% V]=O(m}). 
We then insert in Eq. (2.6) a complete set of physical states, make use 


of the cluster decomposition (see Appendix) and, selecting the pion state, 
find the sum rule?: 


<bIE,|a) + (i/m,) <b| X,|a> 
=) (0105 |n?) <7, b| Vapor 


+ 2 [<b] OF |.) <N IV Ja) — ct] 


EF (- im, fa) (27)? 6(p,— Pr) <blx*In) <n|V lay 
n#N,n (Es — En) (E + m, — En) 


3 From the saturation of axial charge commutators we get, in general, sum rules in 
which combinations of different pion amplitudes appear. It is easily realized that the 
“charges” defined in Eq. (2.1), besides having the symmetry breaking already incorporated, 
enable us to obtain, after use of completeness, a single pion amplitude. 


(2.5) 


where 


(2.7) 


ct.|, 
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), representing the contributions of the intermediate states which are 
N 


degenerate with the external ones (if they are present). In practice the 
E, term will consist of a well-known current algebra commutator (finite 
in the limit m,—-0), plus eventually higher ones, contributing an O(m,) 
to Eq. (2.7); X, will be again a higher commutator, affecting the O(m,) 
part of Eq. (2.7). Of course, all higher commutators are outside the current 
algebra frame, and further, more specific hypotheses are needed in order 
to have information about them. 

Once we assume the l.h.s. vertex functions to be known, we are 
then led to the following representation for the physical pion ampli- 
tude T,,.y,24,,» With the pion at rest*: 


Ty + vone+5 = Xn", b|V ay =(1/f,) b| Ela) 


+ (i/m,f,) <b|X,14> — (A/F) LEDIEN Y <A [Vay - ct] 
N 
(2.8) 
A Mm, f EACH = m,)] Z dg 0*(go; q= 0) , 
o%(q)=i(2n)> % ôP +a- ph) <blx’In><nlVla)- ct. 
nFN,n 

Eq. (2.8) is the general formula we shall use in the following, when 
treating pions at rest; so, let us briefly illustrate its most interesting 
features, referring to Refs. [1,2] for more details. 

A very important point to be stressed is the analogy between Eq. (2.8) 
and a dispersion relation, with both energy and “mass” varying, sub- 
tracted at the point gy =0, the subtraction constant being specified by 
the equal-time term and the continuum integral contributing on O(m,)°. 

As it is easily ascertained, energy and “mass” vary along a certain 
line, connecting the “soft pion point” (q = 0, m, = 0) to a “physical point” 
(q =0, m,). Our Eq. (2.8) is thus to be regarded as a possible prescription 
to extrapolate soft-pion theorems, rigorously valid in the limit of 
vanishing pion four-momentum, where the amplitude is simply expressed 
by a current algebra commutator plus the nucleon term, into the physi- 
cal region. The outcome, of course, can be written in the form T(m,) 
= T(soft)+O(m,), and the nature of the corrective terms is explicitly 
displayed. We have, as already anticipated, higher commutators, in- 

* This is due to the fact that a charge, being an integrated operator, can be connected 
only with pions of vanishing three-momentum. 

5 This is not the only analogy suggested by Eq. (2.8). Remembering Low’s treatment 
of pion physics, still based on equal-time commutators and completeness [3], we can see 
that Eq. (2.8) is substantially a Low representation for the amplitude T,,)W.,«+,, where 
the identification &,.«>0” has been made, and the subtraction constant specified. Another 
remarkable resemblance can be established with low-energy nuclear physics, where one 


has the potential part plus the resonance scattering: it can be shown indeed that the soft- 
pion amplitude is the analogue of the Born approximation in potential scattering. 
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volving time derivatives of the axial charge, which appear since the cor- 
rections to the symmetry limit m, =0 depend on the way the symmetry 
is broken, i.e. the pion acquires a finite mass. A further O(m,) is given by 
the continuum integral, which represents the higher states contributions. 
Actually, it must be emphasized that such an O(m,) might be, as a matter 
of fact, an O(m,(In m,)”), due to possible infrared divergences of the 
integral (see also Ref. [4]). In other words, analytical dependence of 
Eq. (2.8) on m, is not granted, and, consequently, simple power-counting 
in m, would not be meaningful: the only general statement is that the 
“corrections” to the soft-pion theorems must vanish as m,—0. Further- 
more, the evaluation of the continuum is in general a difficult task since 
the vertices <b|y*|n> =F*((E,— E,)?) which appear in g% are energy- 
dependent and, moreover, at timelike momentum transfer. In practice 
however, as a first approximation, this dependence can be neglected: the 
feeling is that such an approximation affects only the less important part 
of the amplitude. 

The next fundamental point concerns the choice of a reference frame. 
We are working out a configuration where the pion is at rest, but p,, Pe 
are still arbitrary and, in principle, any point of the physical region can 
be reached. In other words, the sum rule Eq. (2.8) is frame-dependent 
and we can take advantage of this to enhance (or to depress) certain 
contributions with respect to others. We will choose nearly rest frame 
saturation. The reason is that in the configuration p, = p,=0 a selection 
rule holds, due to parity and angular momentum conservation, which 
strongly reduces the spectrum of intermediate states, allowing only 
s-waves at the vertices <b|Q%|n> and <n|Q%|a). 

In general higher waves will be depressed if both |p,| and |p,| are 
small. We obtain thus a very neat separation between the s-wave and the 
p-wave (and higher), whose presence is eventually related to small 
(nonvanishing) values of p,, p,. While keeping this general limitation in 
mind, we can choose the reference frame at our convenience, depending 
on the process under investigation. In the pion-nucleon scattering the 
choice is immediate, since the configuration p, = p, = 0 is allowed in this 
case, and we reach, of course, the physical threshold. In the electroproduc- 
tion process, owing to the more complicated kinematics, many pos- 
sibilities are available, among which the Breit frame p,= — p,=p will 
be chosen as the most convenient one, for symmetry reasons. 

Our attitude, in conclusion, is to look for a configuration of the 
external particles, such as to have the simplest “correction” terms and 
a good control on them; in so doing the “best” extrapolation will be 
obtained, in the sense that the final result is believed to be as close as 
possible to the soft pion limit m,=0. From the point of view of the 
SU(2)x SU(2) symmetry this evidently corresponds to enhance the 
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Goldstone mechanism (massless pions, isospin multiplets) with respect 
to the other possible ones. 

The best illustration of the general ideas discussed so far is perhaps 
their application to specific pion processes. Let us devote then the second 
part of this Section to pion-nucleon scattering at threshold, reserving 
Sections III, IV to the photo- and electroproduction processes. 


b) Pion-Nucleon Scattering 


We have to take in Eq. (2.8) |a), |b) nucleon states with momenta p,, pz 
and V=Öh. 
The pion term is then: 


5 (2m)? f da. E2' <0| Oz 1m) <x, p2|O6|P1> 
= (fz /m,) (27)? ö(p2 — Pi) Tex 


where TZE = <n", p|x’|p> is the forward pion nucleon amplitude with 
the pion at rest. We specialize moreover to the nucleons rest frame p = 0, 
reaching thus the physical threshold configuration, and introduce the 
s-wave scattering lengths a’’: 


(2.9) 


T2f (th) =82(M+m,) a. (2.10) 


The current algebra information (or, if one prefers, the soft-pion result) 
is contained in the well-known SU(2) x SU(2) commutator 


P| [0*, OP] Pp >= 146) P2 |Q"| p> 
= i&,g,(t’/2) (27)? 2 E5(p2 — P1) 


(2.11) 


and, in addition, the higher commutators ID, Ollpy®, <p\LO, Öllp> 


and <p|[O,O]|p) are required. 
It is convenient to introduce the decomposition into even and odd 


part (under «> f): 
T =ô p 17044150170". (2.12) 
Collecting then Eqs. (2.8)—(2.12) we find very easily 
(1+m./M)a = L(1 + ôa), (2.13) 
(1+m,/M) a*”’=Löa*, (2.14) 


6 Between states of equal four-momenia, <{Q, OD =- <Ë, öl. 
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where a = $ (a, —a;), a’) = $ (a, +2a,), and L=m,/8n.f2 ~ 0.09 mz! 
is a universal length, independent of the mass of the target. Moreover’ 


dgoo‘” ‘God =0) 


da = (2m? f2/nM) | +m? 0,0), (219) 


ine  g0(95 — Mz — ie) 
5a‘) = (2m2 f2/nM) | dgog' (Go q =9) 
ma 40(4o — m; — ie) (2.16) 
+m,' <L, ÕI) +m <À, ËI), 
efla) = 5) dxe <p E), Olp; p=0. 17) 


We have thus a theoretical prediction for the pion-nucleon scattering 
lengths, in terms of equal-time commutators and a continuum (only 
s-wave excitation !); it reduces of course to the well-known Weinberg- 
Tomozawa [5] result in the limit m,—>0 8°. 

Actually pion-nucleon is a very unlucky case, because, while we can 
be confident in an approximate evaluation of the continua, which can 
be handled in a relatively simple way, the nature of the higher com- 
mutators, on the other hand, is quite complicated. Apart from the 
o-term <N|[Q,Q]|N)>, which is very interesting since it gives informa- 
tion about the transformation properties of the symmetry breaking, we 
need the [D, D] and [D, D] commutators, on which we cannot make any 
general statement at all. We are inclined to use field-algebra like models, 
since the present treatment of the low energy region is roughly similar 
to a Lagrangian one, once the identification ¢,.«+ D* is made. In so 
doing we can rely on canonical quantization rules!°, and hope to get 
rid of both [D, D] and [D, ny (actually, the latter is interaction dependent). 
In this respect the electroproduction process is much simpler, since, as 
we shall see, the problem is particularly insensitive to higher commutators 
(involving the breaking). 

7 Owing to the selection rule introduced in Part a), the nucleon cannot contribute. 

® We remind that it may be misleading to ascribe a simple power dependence on Mm, 
to the corrections. Logarithmic dependences, or even a lowering of the power are not 
excluded. 


? It is perhaps worth to mention the consistency condition we get from comparing 
with dispersion relations: 


VI TU + Om) = gy /(M? — (m4 M) + Q/n) | ag im TS? 
0 Ga 


2_ M?°m 


Performing the limit m,—0, and using the Goldberger-Treiman relation, the Adler-Weis- 
berger sum rule follows: 


r3’ =1+@M?/g,n)|v'dvlota"p-otr*p)]. 
1° This is impossible in other models, e.g. in the quark-model. 
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If we try anyway a comparison with the experimental results'!?, for- 
getting the correction effects, we find a very good agreement: 
(a) = 0.166 mz! (a,)”? = (0.18 + 0.005) m; ' 
(a,)°°" x — 0.083 mz!  (a;)”? = — (0.088 + 0.004) m'. 


11b 


(2.18) 


As far as the continuum corrections’*” are concerned, the complete 
structure of the spectral functions o} (© is easily established, by using 
the cluster decomposition, and the result is da = ôa; + day + day), where 
ôa; n,m represent, respectively, the contributions of connected, semidis- 
connected and Z-graphs. Detailed calculations of these terms are not 
feasible, but, however, indications about their order of magnitude can be 
obtained, taking into account the structure of the integrals in Eqs. (2.15), 
(2.16) and the masses of the lowest isobars [2]. It turns out that, in practice, 
it should be reasonable to keep only the leading contribution to ôa, i.e. 
the x — N s-wave intermediate state (only 3” states are allowed in öa,), 
and assume the integrals to be dominated by the threshold region. The 
outcome is, using elastic unitarity 


ôa? =(2m,/3 nL) (az — a3) 
alt? = (m,/3 L) (a? + 2.43). 


Inserting the experimental numbers, this is seen to affect very slightly 
the soft-pion predictions quoted in (2.18). 

An alternative attitude could be to try to infer, from the knowledge 
of the experimental data and from an approximate evaluation of the 
continua, information about higher commutators; in particular, about 
the o-term. This problem has recently received much attention, and is 
really a controversial question. So, although not willing to enter into 
details, let us end this Section by briefly commenting on it. 

The first indication came from Kim and Von Hippel [6] (see also 
Refs. [7, 8]), who generalized the above approach to meson-baryon 
scattering at threshold, in the framework of the SU(3) x SU(3) descrip- 
tion!?. The role played by the continua appears to be really crucial in 
this case; drastic approximations and uncertain experimental data may 
cast doubts on their calculations. However, if we believe in their results, 
we have indications that the (3, 3) +(3, 3) model of SU(3) x SU(3) sym- 
metry breaking [9] fits, and, furthermore, ø ~ 30 MeV. This is of course 
the size we expect if SU(2) x SU(2) has to be a better symmetry (namely, 
less broken) than SU (3) (remember that o ~ m2). 

lia References and values based on more recent data are given in the compilation of 
Ebel et al. [Nucl. Phys. B22, 317 (1971) (Ed.)]. 

116 Strictly speaking, the term “corrections” should apply to a‘~ only, whereas a‘*? 
is given by a sum of terms which are formally of the same order, as is evident from 


Eqs. (2.13)-(2.16). 
12 We refer to Ref. [9] for what concerns the general scheme of SU (3) x SU (3). 


(2.19) 
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On the other hand, an alternative approach to the o-term, based on 
fixed t dispersion relations for the pion-nucleon isospin-symmetric amplı- 
tude TĘ, seems to indicate a quite larger value, namely « ~ 100 MeV 
[10] (which, by the way, is the order of magnitude of the SU(3) breaking). 
The assumptions and the hypotheses made in this kind of calculations 
are of different nature, but the reliability of the results is certainly not 
improved, since, among other things, reasonable uncertainties of the 
experimental input lead to strong variations of the final outcome’?. 
Anyway, a large value of the o-term would be perhaps welcome by 
theorists (actually, it has been, and papers on that are becoming more 
and more numerous), since it could be a fruitful starting point for further 
speculations concerning chiral symmetry. Among the many promising 
ideas let us mention, for instance, the possible connection with a 
(spontaneously) broken dilatation symmetry'?. This is why the problem 
of the o-term, which, in our opinion, is still an open question, deserves 
further investigations: a better understanding of chiral symmetry and its 
breaking will perhaps result. 

As a final remark, let us emphasize that the approach we have 
presented here to perform the extrapolation from soft to physical pions, 
is, in Some way, the most economical and nearest to the original current 
algebra framework. In particular the minimal number of commutators, 
beyond the ones of the SU(3) x SU(3) algebra, has been introduced: 
a prize to pay for this is the very nature of the continua. It is clear that 
an improved knowledge of higher commutators and a more detailed 
information about them would enable one to have more and more sum 
rules, from which to eliminate the unpleasant contributions. We have in 
mind recent approaches to the extrapolation problem, starting from the 
knowledge of light-cone commutators: such an information (that roughly 
amounts to an infinite number of equal time commutators) seems to 
allow a different representation of the continuum corrections, but their 
estimate is anyway not easy [13]. 


IH. Electroproduction at Threshold 


As a further example, we now consider the electro- (and photo-) produc- 
tion of a single pion. This process has represented a test of very different 
approaches, and to try to review the large number of related papers 
would ask for a separate effort [14]. From our point of view it will be 
enough to recall the relevant “soft pion theorems”. In photoproduction 


13 The dispersive evaluation of Ref. [11], for instance, yields the result ¢~40 MeV. 

Note added in proof. The result of Ref. [11] is essentially confirmed by a calculation, 
in which the new CERN phase shifts (Oct. 71) were used (H.P.Jakob, CERN preprint 
TH 1446) (Ed.). 
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the standard application of a method!* due to Low [15] leads to the 
so-called Kroll-Ruderman [16] theorem for the s-wave multipoles, 
stating that: 


lim ES n= (1/2M) gv, = 0), 
; (3.1) 
Jim ES} =0, 
where the usual isotopic decomposition of the photoproduction ampli- 
tude 

Ty = 6.3 Th + $0, l- Th +2, TO. (3.2) 
has been used. 


When the photon is off mass shell, we have the following formulae, 
derived by Nambu and coworkers [17]: 


4M? —x? g,n(0) | G,(x?) K? , 
(~) = aN A v 2 
ESBE) 4M? 2M | G,0)  4M? 2 Guli ) 
(3.3) 
4M? =r? genl) | — r? GR?) 
(0),(+) as nN M 
ER O(S.P.L) ar 6A 
_ 4M?—x? GE(K?) 
EL | er 99 
4M? — x? — G5’(x?) 
LYM Ko (S.P.L,) = Im Gnn(0) ino, (3.6) 


and for a quick derivation of them, which also summarizes the main 
features of the fundamental work of Nambu, we refer to a previous paper 
[18]. In Eqs. (3.3)-(3.6) and (3.1) Ey, Yo, are the s-wave multipoles, 
defined as in Ref. [18], and a definition of the different form factors 
entering the various expressions will be found in going through this paper. 

A non trivial point, which has been first emphasized by Nambu, is 
the frame-dependence which the soft pion limits may in general exhibit. 
From a point of view of dispersion relations, this happens as a rule for 
all those quantities which receive contribution from the Born approxima- 
tion. Actually, a prescription to reach the situation where the pion has 
both vanishing mass and momentum has to be specified. Usually, one 
means by “soft pion limit” the limit attained by setting first the pion 
three-momentum, and then the pion mass, equal to zero. In this way, 
a frame is automatically chosen; in particular, Eqs. (3.3)—(3.6) are derived 
by this limiting procedure in the final pion-nucleon center-of-mass 


14 A rigori, this method would lead in general to “soft photon” theorems. However, 
if the photon is real, it reduces to the form Eq. (3.1). 
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frame, where the multipoles are defined. (This should be kept in mind 
if one wanted to compare our final formulae with Eqs. (3.3)-(3.6)). 

To begin with, we consider the simple configuration where the final 
pion is at rest. The relevant commutator, to be inserted between two 
one-nucleon states of three-momentum pz, p, respectively, is then: 


[Oi V] = Ez- (3.7) 


Here V, is the electromagnetic current, V, = V{® + V{° and Ez is the sum 
of a standard charge-current commutator 


[0% V,] =ie,3,A?, (3.8) 
and of a term explicitly depending on the breaking. 
[0% V]= Jdx[D%, 0), V] (3.9) 


(for the time component the relation [Q*, Vy] =ie,3,D’ follows immedi- 
ately from Eq. (3.8) by the electromagnetic current conservation). 

We do not make at the moment any assumption about the com- 
mutator [D*, V]. Its role in the determination of the s-wave multipoles 
at threshold is practically irrelevant, whilst it has some importance in the 
case of certain higher multipoles. Therefore, in the first part of this 
Section, let us simply ignore it. 

To proceed on, it is necessary to fix the momenta of the external 
nucleons, i.e. to choose a definite frame. For a series of reasons partially 
explained in the previous Section, it seems to be very convenient to choose 
the Breit frame of the two external nucleons, p, + p, = 0. 

Following the rules of the game illustrated in the previous Section, 
we then select the contribution coming from the one pion state and are 
led in a straightforward way to the expression 


1 iu 
Tg, t K’)sa. = A <p2| Egl Pi? — 2 (p21Qrin) <n|V,|pı> + et.) (3.10) 
where T% is the physical single pion electroproduction amplitude 
Ti = 1 CN (po), 2°(q = 0)| V1 N(Py)> (3.11) 


evaluated in the kinematical configuration where the pion is at rest in 
the Breit frame, which means, in a covariant notation, at the point: 


vemE=v.. 
t=K?—m=t (3.12) 
where (p, +K = p’ +q): n PBL. 


v=P:k=P-4; t=(p -p 


3.13 
P=ż}(p +p); P?=M*-t/4; E=|/P?. 01) 
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Although this is not strictly necessary, it is convenient at this point to 
adopt a fully non-covariant notation and effectively evaluate Eq. (3.10) 
in the frame P=0. To this purpose we need the explicit expressions of 
the nucleon-axial current and nucleon electromagnetic vertices, which 
are respectively 


<p|A*| — p> = (07/2) {2EG, (0) [n(o -n)— 0] + no-n D()}, 


a Bi FERN (3.14) 
plAl-Ppr=0; <pID*|— p> =(1°/2ì) ox D(t), 
where 
D(t)= ~2MG,(t)+tGp(t); n=x/|k] (3.15) 
and 
EIN — py = ż <1°, 1) 2M Gg (0) (3.16) 
<p VO) — py = Ir’, 1> io xk GY (t), (3.17) 
where Gg m are the Sachs form factors: 
Gy=F,+F,;  Ge=F,+(t/4M’)F. (3.18) 


Concerning the amplitude T;, we choose its space components as follows 
T=Ao+B[k(-x)—x’o] (3.19) 
and define its time component through the gauge-invariance condition: 


K: T= Ko To (Ko = My) « (3.20) 


By comparing the time and space components of Eq. (3.10) it is now 
easy to deduce the equations, where a “real part” is always understood: 


A = D(t/2f,+6A™, (3.21) 
Alt) (0) — mM a (0) GES(t) + 6.AT»O (3.22) 
ZEISS Í 
a 1 = G,(t) 
= (A+B = — AZ 
tG 4(0) GA) 1 pe) 
A M = 
= 6A+töB]” 
i 2E OATB, 
Tit =m, M G 4(0) Gi (O/4 EB’ fx 
(3.24) 


k 1 
+im, BOONE + ap 6A + tôB] HO 
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where 
> dqoes b(Go, d) 
AO), BO) =2 OSAR 3.25 
ot Mr J 90(95 — mz — ie) ’ ena 


ALO = 2m | f doga” (aot) ; (3.26) 
o qolqo — mz — ie) 


j Sate Mdo, t) 


(3.27) 
— m? — ie) 


5B — 2m, | 
0 


Q4, g are the components eo to the decomposition Eq. (3.19)) of 


o% = i(2n)° Pee Px<P\X'|n><n|V|—p> (3.28) 
and the quantity m,B“) is the contribution to B+» from the 
commutator i/m,<p|[Q%, V]|— p> (it does not contribute to the real 
part of B owing to the crossing properties of 9%). 

The reason for the choice of the particular combinations (3.21)-(3.24) 
is that, following the notation commonly adopted in the literature [19], 
one can easily verify that (we are now in the c.m.s., where Wand 6, are 
the total energy and the pion angle respectively): 


(ll/ll) (doz/d coSO,)p., =(0E7/2W*)/|Trp.l?, (3.29) 
(x |/Iq)) (doz/d.cosO,)p.. = — (° 0/8 W° mz) An... (3.30) 


According to the remarks of the previous section, it is clear that the still 
not evaluated terms in Eqs. (3.21)-(3.24) will be vanishing together with 
the pion mass. If we neglect for the moment the fact that the actual 
corrections can lead to a pion mass dependence of the corrections 
different from the explicit powers multiplying the continuum integrals, 
then, in the specific case of photoproduction, we can summarize our 
results in the simple form: 


G40) 


Th (xk? =0;B.t)= — Of. + O(m2). (3.31) 
G40) 
O2 0: Bt) aa) 3 3.32 
Tr (k* = 0; B.t.) 4M f, +O(m;). (3.32) 
Introducing the physical amplitudes 

y 

T, =I ro; POR 

yn>n p 
: (3.33) 

YP p 


Tor= TTO; 
02 a ynon°n 
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we find: 
(\«\/lq|) do, /dQ,,. = 12.7 ub/sr +O(m2); Exp.: 15.6 +0.50% 
R=do_/do. = 1.35 + O(m?); Exp.: 1.265 + 0.0752 
~ 1.3522) (3.34) 
(\K|/Iq))doo,/dQ,..= 0.04 ub/sr + O(m2); Exp.: 0.07 + 0.0229 
(lK\/Iq\)doy2/dQ, = O(m?) 


whilst the corresponding predictions of the Kroll-Ruderman theorem are 
respectively!®: 19,5; 1; 0; 0. Since it is immediate to verify that the Kroll- 
Ruderman theorem is actually reproduced by Eqs. (3.31), (3.32), we can 
look at them as a very simple and meaningful generalization of that 
theorem. The terms respectively classified as O(m2) and O(m?) are 
therefore, as we expected, “corrections” (unambiguously defined) to the 
soft pion theorem Eq. (3.1) (although the final dependence on m, can 
turn out to be different, in particular Inm, can arise). 

In a perfectly similar way, the analogous terms in the general equa- 
tions (3.21)-(3.24) evaluated at x? +0 will represent “corrections” to 
some results, which would be exact in a world of massless pions. The 
corresponding “soft pion formulae” are different from Nambu’s Eq. (3.3) 
(this is not surprising if we remember that the pion has now been put at 
rest in the Breit frame, and not in the c.m. frame). 

A pleasant feature we want to mention here is that, if one imposes 
that the same results have to be derived by starting from a dispersion 
theory approach and by going to the limit (q=0, m,=0) in the Breit 
frame, then one obtains “consistency conditions” which are just the 
Fubini, Furlan and Rossetti [24] and Furlan, Jengo and Remiddi [25] 
sum rules?°. This is just the analogue of the result shown in the previous 
Section, the “consistency condition” being in that case the Adler-Weis- 
berger sum rule. 


15 We have used the value G,(0)/2 f, =0,92/m,. 
16 The sum rules are 


G4(0/G 40) = Fit) - (2 Mt/g,un) | dv-v Im M§(v, t, g°=0, x? =2), 0 


cont 


Fr)? M=(2 M/g,ym) | dv- v~t- Im Mr, t, q?=0, K? =t), (II) 


cont 


where M, , are the invariant coefficients of $ ys [y 8, y} k], ys(k-8y x — K?y - e) in the elec- 
troproduction amplitude. Usually the sum rules are exploited to get an estimate for 
F3:5(0) and G; (0), but one can also read them as soft pion theorems, i.e. 


Mo, g? =0, K? = Non = (gen/2? MP), a) 
MS, t, g? =0, K? = Dloone = (g,n/2 Mt) [FF (©) — G4(0/G4(0)] . (Iv) 


We shall comment on them in the following. 
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To give a more concrete meaning to Eqs. (3.21)-{3.24) an evaluation 
of the “corrections” is required. Clearly, if such an evaluation could be 
exactly performed in terms of known parameters, then from the measure- 
ment of electroproduction at the Breit threshold one would be able to 
learn something on the nucleon axial vector form factors, which con- 
tribute to the (—) isotopic amplitudes, that is to the charged pion cross 
sections. 

Since the terms to be evaluated are related to off mass shell quanti- 
ties, their treatment requires unavoidably some “smoothness” assump- 
tions. To minimize the role ofthe extra hypotheses needed, a first prescrip- 
tion is to choose a “small” t configuration (say, |t|< 20 m2). This is 
suggested by the fact that, at t=0, the selection rule illustrated in the 
previous Section would forbid p-waves and higher waves to contribute 
in the direct channel, so that their contribution at small t should be 
particularly depressed. Correspondingly, we shall obtain predictions 
for electroproduction in the range 0 < |x?| S 20m. In this region, an 
evaluation of the “corrections” has been proposed [26], whose main 
features can be summarized as follows: 


a) The contribution of the extra terms is a significant fraction (about 
20-30%) of the “leading” (i.e. equal time and one-nucleon) contribution 
(which would be sufficient to reproduce the “soft pion limit”) in the case 
of the (—), (0) amplitudes. However, it appears still licit in this case to 
treat the extra terms as “small” corrections, to be evaluated by making 
recourse to a number of unavoidable approximations. 


b) The contribution of the extra terms is dramatically larger than 
the “leading” one for the (+) amplitudes. Therefore, all the predictions 
for neutral pions depend in a rather crucial way on the approximations 
used to evaluate the unphysical contribution. Owing to this reason, we 
shall consider the latter predictions as purely indicative, and shall limit 
ourselves from now on to treat the charged pion processes. Before doing 
this in more detail, it is fruitful, however, to realize that these formulae, 
although the simplest, are probably not the most popular ones. First 
of all they hold at the Breit-threshold (corresponding to pions emitted 
backward in the c.m.s.) where a small amount of p — and higher waves is 
present. 

A prediction concerning the actual physical threshold would be 
probably more favourably accepted. Secondly, a generalization from a 
single-point prediction to a prediction covering a given energy-range 
would allow a more complete check of the whole approach. 

If we do not want to loose the pleasant features of working in the 
Breit frame, a simple solution is to have predictions for pions which are not 
at rest in the Breit frame, but have just that amount of three-momentum 
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which corresponds, by a Lorentz transformation, to being at rest in the 
c.m.s. 
This aim can be attained by introducing the operator 


Ar(q) = Abla) + io(q)mz? Da), w= +m} 


E j y _ (3.35) 
Aig) = [dx exp(—ig x) Ai(x) — Ot 
with the property 
0 sa’) = 2i@f,(27)* 545 5(q — 7) (3.36) 
all) =0 
and starting from the commutator: 
Lig), V,(- k)]= E + iom? X" (3.37) 
where 
E” =[ 45, Va] (3.38) 
and = 
X*=[D*,V,]. (3.39) 


As we shall see, we will be able to deduce from these commutators an 
information for s, p and d-waves which covers a kinematical region 
ranging from the physical threshold up to the beginning of the first 
resonance, thus avoiding also the second previously mentioned restric- 
tion. 


IV. Low Energy Electroproduction and Concluding Remarks 


A point to be stressed before we proceed is that if one wants to start 
from a commutator like Eq. (3.37), involving a charge density rather 
than a charge, then of course many features of the approach outlined in 
the previous Sections will change. First of all the chiral symmetry 
pattern will be partially lost. By this we mean that one will have now two 
kinds of predictions. The first ones would give back in a particular limit 
all those results which can be derived by starting from a charge-current 
commutator as Eq. (3.8); therefore they will still contain a term which 
reproduces a “soft pion limit” and some “corrections”. The second ones 
are new predictions it is not possible to obtain starting from a charge, 
and thus they do not correspond to any “soft pion result”. Roughly 
speaking, they involve those invariant amplitudes whose coefficient 
vanishes in the limit q,—0. 

The amount of new information one obtains in this way must corre- 
spond to some new imput, and therefore to a loss of generality of the 
approach. In our particular case the use of the current-current commutator 
implies that possible Schwinger terms appear in [A5, V], which are 
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model dependent. In a simple field algebra model which is suggested by 
the analogy of our approach to a Lagrangian description of low energy 
pion physics, the mixed Schwinger term in [A}, V] is zero. Other models 
more generally [27] predict the existence of a gradient term: 


[AG(x, 0), Vily, 0)] = iex3, Aula) 5(x — y) - iS44(y) gu 6.ölx — y) (4.1) 


provided that: _ u 
[D*(x, 0), Voly, 0)] = i€,3yD"(x) ölx y) ` (4.2) 


In general one can easily verify that, if we saturate Eq. (4.1) in the 
Breit frame P = 0, then this term will contribute to only one appropriately 
chosen invariant amplitude. But gauge invariance will allow us to 
replace this amplitude with a combination of the remaining ones, so that 
within this approach the nature of the Schwinger term in Eq. (4.1) is 
irrelevant. 

It remains to discuss the commutator Eq. (3.39). A pleasant feature 
of the saturation of densities is that it allows us to see to which multi- 
poles this commutator actually contributes. If we limit ourselves to 
charged pion processes, i.e., to the amplitudes (—), (0), then it can be 
shown [18] that this commutator is practically irrelevant to the determina- 
tion of the s-wave multipoles'’, whilst it will contribute in general to 
higher waves. Therefore, if we want to push our predictions beyond 
threshold, we must make some definite hypothesis concerning it. A 
simple solution, suggested by field theoretical models, is that this 
commutator is vanishing: 


[Dz V]=0. (4.3) 


We shall accept Eq. (4.3). From the point of view of dispersion theory, 
it can be seen [18] that this assumption has at least the same validity 
as the statement that the Born approximation gives a “reasonable” 
evaluation of the (—), (0) invariant amplitudes at threshold. 

The technical details are now not very important and we shall limit 
ourselves to a sketch of the procedure. The starting point is represented 
by the Breit frame saturation of the commutator Eq. (3.37). The procedure 
should be by now familiar, and the result can be written as [28] 


Sa TAP, 9) = KLAG, VD 
+ (io/mz) <LD*%, V,]> + (one nucleon) (4.4) 
+oRi+q- Si 
where R}, Si, are directly connected to completeness sums. The factors 
œ, q clearly indicate that we are extrapolating a result, originally obtained 


17 More precisely, it contributes to a kinematically depressed part of EP, which, 
in turn, is a small quantity compared with EṢ}. 
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in the limit q, m,—0, to a whole range of the physical region. Eq. (4.4) 
allows us to derive, after appropriate projections, the complete set of six 
invariant amplitudes at a certain physical point, specified by q. In 
particular, we shall choose q so as to obtain that kinematical con- 
figuration which corresponds to the physical threshold. The next step 
is then to transform our language and to introduce the experimentally 
meaningful multipoles. In the low energy region we want to consider, 
namely from threshold up to |g°™'| sm, we are allowed to use, for these 
quantities, power series in |g¢™ | truncated at the lowest power. Moreover, 
even if, in principle, Eq. (4.4) could be exploited in a general configura- 
tion, we prefer to fix 0 = 90° (in the c.m.s.), for simplicity reasons. We are 
able thus to describe electroproduction at 0 = 90°, in an energy region 
from threshold up to the starting of the first resonance, and for not too 
large momentum transfer between the electrons (say, |k”| < 20 m2). 

To be more definite, a number of final formulae is unavoidable. We 
first write down the expression of the differential cross section, including 


terms of order |go™|?: 


do a EĘ kW M 1 gd 


{[A + B cos + C cos? 0] + (e(— x?)/x2)[D+ E cos0 + F cos? ô] 


+ecos2 sin? 0G + 2e(l + £) (— K°)/K cos sin OCH + I cos0]} 
(4.5) 


where EF, E}" are the lepton energies in the laboratory system, Q, = (0, $) 
is the x — N solid angle and ¢ is the “polarization” defined as: 


e= [1 —(2|«"|?/x?) tg? 0/2] 
A, B, C are given by the following equations: 
A=F} + $2? +y?+3yz+3xEo- 
B=2xE + (4.6) 
C=—A4+F?+x?+15Ey. yp 


where 

Eo+ =Eo+(agl=9); Fy, =|Eo+ +3M2- +E,- —3M2+ —2E2+|. 
z=Re(E,+—M,;). y= Re(M,+~E,+—M,-—-E,-). 
y=Re(2M,.+M,-). y= Re(2M2+ + E24). (4.7) 


x=Re(3E, + +M,+—M,-) 
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and for the remaining coefficients we have: 
D=F2 +0; E=2u2,-; F=uW—-12+45 2: LZ); 
G=3z?+3yz+34Ey:; H=% +32)+ Eo-b; 
Er ae ae as 
=| Po +22- -32l Fo = Zo- (igal =0); 
pares +H,-; v0=/,--22+; 
=62,--94,:; C=12E,.—3E,-—3M,-—3M,:. 
The quantities which are needed to describe the process in our con- 
figuration 0 = 90° can be derived by a lengthy but straightforward com- 
putation. Since their explicit expression is rather complicated, we limit 


ourselves here to give the s-wave multipoles at threshold, and refer to 
Ref. [18] for a complete list of the results. 


Bg? (th) = —|/P?/M?(G4(0)/2f,Jly— VP?/M? [4P?f,/1 — x qP 


` [(x%,/2) G ACH) Gig (IR) z (M/P?) v5, Gp(qx) Gul KN Im + JEG? > (4.9) 


E (th) = Pn Cala Files) + |/P?/M? [4P?f,//1—x-q'/P?]~' 


2M q'/P? 
[Gar F5 cr) (1? PP AM?) + m, vy Fi (ic) + Falh) t/4M°)] + 5 BY? , 
(4.10) 
HP froda = Venen //P?/M?[4P?f,/1—K-q/P?]! 
MD G 
a +5 LD Keo , (4.11) 


Prola =V MIP GIO , D VPP ETP 


r 


M vr q Men 
2M+m, |\ P?” M 


GB) u) (4.12) 


mV A 
~ TEP GR PD +82 


M(x? —m2) 1 
in “Mem Ky t; qn = O(m;); vy = O(mz);K-q' = om) 


The “corrections” are still vanishing with the pion mass, and one 
can easily check that Eqs. (4.9}{4.12) reproduce the “soft pion limits” 
Eqs. (3.3){3.6), of which they are therefore a meaningful and unambiguous 
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extension. The higher waves represent, of course, a new result, in the 
sense that, in this case, there is no low-energy theorem to extrapolate. 
We have now to evaluate the “corrections”. This evaluation has been 
performed in previous references [28], and we do not want to quote 
their complicated expressions here. From a general point of view we can 
say that there are some “corrections” which can be very reasonably 
evaluated without making recourse to particular assumptions, i.e. the 
s-wave continuum and the N* contribution in the s-channel (the latter 
being always extremely small). In addition we have to evaluate more 
complicated terms. Our empirical assumption is that the main contribu- 
tion comes from the vector meson exchange in the «?-channel and that 
other graphs, e.g. the 3x contributions, can be ignored. 

In this way we are led to a representation of the various multipoles 
containing a “leading” current algebra and one nucleon term and 
“corrections” whose typical order of magnitude is a 10-20% of the 
overall quantity. The next important point is then the comparison with 
experimental data. This may be done only if we choose a definite expression 
of the axial vector form factors (we suppose that a dipole fit with My =6m,, 
well accounts for the electromagnetic ones). In particular, if we choose 
the following parametrization (we work in the small t region) 


G (t) =(1—1/M3)-?G,(0), (4.13) 
D(t) = 2m? f,gqy/(m2 — 1) (4.14) 


then, fixing G ,(0) = — 1.23, we find, corresponding to the values M , = 6,7, 
8, 9m,, the results for the slope of the photoproduction differential cross 
section at threshold shown in Table 1, where the experimental results for 
(\x\/\q\) do ,./dQ [20] and o_/o, [21] are also quoted. Concerning the 
latter, however, we have new indications [22] which seem to point 
towards the value R = 1.35. 

In the electroproduction case, a recent threshold measurement [29] 
corresponding to the values 

k?’ = —10.4m2; EF=800 MeV; 


4.15 
E; =542 MeV; ¢=0.74 en 


Table 1. Slope of the differential cross section at threshold (ub/sr) for x*-photoproduction 
and n*/n~-ratio, assuming G,(0) = — 1.23 


Ma 6m, 7m, 8m, 9m, exp. 

oe 15.2 15.6 15.8 16.0 15.6405 

gl dQ 

R=o_je, 1.35 1.35 1.34 1.34 1.265 + 0.075 


1.35 
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Table 2. 2*-electroproduction at threshold in units of 107°! cm? (GeV/c)~*sr7+ and 
6,/67-ratio corresponding to the value G,(0)= — 1.23. k? = —10.4p?, EF = 800 MeV, 
Ei = 542 MeV, ¢=0.74. 


Ma 6m, 7m, 8m, 9m, exp. [29] 
2 
grape 4.03 5.15 6.00 6.72 4940.7 
q T 1 
0/67 0.53 0.40 0.32 0.28 


Table 3. da(90°)/dQ (unpol.), da (90°)/d@, both in ub/sr, asymmetry ratio Z and 04/0, 
at 90° as functions of photon lab. energy. G,(0) = — 1.23, M; = 7m, 


do ,/dQ = (x M?|qj/4xW? |x|) (A, +B, cos + C, cos?6) + O(g?) 
A,=F?+y*, B,=B=2xE,,, C,=x?-y +1554 


El(MeV) da/dQ do,/dQ X do,/do, do/dQ 
unpol. n* unpol. z7 
160 4.96 5.26 0.059 0.89 6.11 
170 6.38 7.09 0.119 0.79 7.80 
180 7.17 8.34 0.163 0.70 8.70 
190 7.85 9.89 0.192 0.70 9.45 
200 8.68 10.33 0.190 0.68 10.36 
210 9.47 11.11 0.172 0.70 11.24 
220 10.41 11.70 0.134 0.76 12.30 


gives the result: 
lim lg." d’a(r*)/dQ "dE; = (4.9 + 14%) 10” *!cm?(GeV/c)"?sr”! 
ale 

(4.16) 


In the same kinematical configuration, our results are summarized in 
Table 2, where also the o,/o7 ratio is given: From Tables 1, 2 we conclude 
that the agreement with the available threshold data is excellent if we 
choose M,=7m, (and G,(0)= — 1.23). Therefore we fix from now on 
M,=7m, and try a comparison with data beyond threshold. In the 
photoproduction case we find, at 0 = 90°, the results quoted in Table 3 
and one can see in Fig. 1 the way our points fit the experimental [22] 
curve (something like + 5-10% has to be attached to every point) for x*. 
In Fig. 2 the asymmetry ratio [30] is finally given. 

In electroproduction, no data are available till now. We give in any 
case in Table 4 our results for the ratio o,/or at threshold at various x’, 
for different values of M,. 

A general (and obvious) comment has to be added to our numerical 
predictions. They arise, as we already said, from an expression containing 
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160 200 240 
E y (MeV) 


Fig. 1. n* photoproduction differential cross section at 90°, as a function of the Lab. 
photon energy. The dots (@) represent our numerical predictions 


150 200 250 
Ey (MeV) 


Fig. 2. Asymmetry ratio X for x* photoproduction at 90°, as a function of the Lab. photon 


energy. The curves represent dispersive calculations, by Schmidt ( )and by Donnachie 
and Shaw (— — — ). The dots (@) represent the theoretical predictions. For the experi- 
menatl points see Ref. [30] 


a leading term, with a simple dependence on G,(t) and D(t), plus a 
“continuum” part. Once a theoretical estimate of this last contribution 
is available, one can test from a comparison with experiment the validity 
of fits like Eqs. (4.13), (4.14) for the axial vector form factors, in particular 
for D(t). Conversely the excellent agreement obtained in photoproduc- 
tion, where these fits should be rather good due to the small t value 
involved, represents a good check of our theoretical computation of the 
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Table 4. o,/o;-ratio at threshold, corresponding to the choice G,(0) = — 1.23 


K? M,=6m, M,=7m, M,=8m, M,=9m, M,=6m, M,=7m, M,=8m, M,=9m, 


oz/o7(n") o1/o7 (x) 


— 3m? 0.33 0.30 0.28 0.27 0.53 0.49 0.46 0.44 
— 6m? 0.41 0.34 0.30 0.28 0.78 0.67 0.61 0.56 
— 9m? 0.47 0.37 0.31 0.27 1.04 0.84 0.72 0.64 
— 12m 0.59 0.42 0.33 0.28 1.30 0.97 0.79 0.69 


continuum. However for larger values of t the accuracy of those simple 
pole expressions could fail, in particular D(t) could be not well accounted 
for by (4.14). If one believes in the extended validity of our computations, 
then measurements of or and o, at threshold could provide a useful 
information about the effective t-dependence of these form factors. 

The importance ofa careful estimate of the continua can be appreciated 
by comparing our indication M,~7m, with a recent fit for G,(t) by 
Nambu and Yoshimura, obtained by matching experimental data near 
threshold with an improved version of the original Nambu-Shrauner 
formula. They obtain M, ~ 9.5m, and the different determination can be 
ascribed to the different size of the finite pion mass effects (that in their 
case never exceed 10%). 

We are thus led to a quick comparison between our results and other 
theoretical calculations of low energy electroproduction. 


Dispersion Theory 


Clearly the dispersion relation approach, based on the solution of 
multipole integral equations remains the most complete theory. However, 
the practical determination of the multipoles and their dependence on 
the input parameters is not always easy and we believe that in the low 
energy region (below the first resonance) the present approach can be 
considered as a reasonable alternative to the dispersive scheme. In the 
photoproduction case the agreement between the two descriptions 
appears to be quite encouraging, as shown in Table 5, where our predic- 
tions are compared with the dispersive calculation of Donnachie and 
coworkers [31]. As far as electroproduction is concerned, we want to 
stress that the dispersive approach and the present one exhibit an explicit 
dependence on different quantities. For instance in the Born approxima- 
tion of dispersion relations the pion form factor is usually included; it 
does not appear here’®, but on the other hand G,(t) and G,(f) must 
be supplied. 


18 Although there is the photoelectric pion term, accounted for by the pion pole 
in D(t). 
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Table5. y and z at x? =O, and in the evaluation of Berends, Donnachie and Weaver. 
y=2M,,+M,-;z=E,,-M, +; G0) = — 1.23; M, = 7m, 


Ey(MeV) y? a Sal 29 Yodw. bw. dw.  Zbbw. 
160 -0.136 20 0169 -0007 -0.150 =0 0.165  —0.007 
180 -0259 20 031 -0013 -0306 =0 0.299 0.013 
200 -039 ~0 0470 -0019 -049 =0 0405 0.017 
220 -0489 20 0.608 -0.025 —0.590 0.004 0.496 0.019 


Current Algebra 


There are a few recent papers on the subject [32], where a more complete 
bibliography can also be found. The common idea, whose root is already 
in the Nambu’s papers, is to separate the amplitude into the contribution 
coming from the Born terms (pion photoelectric diagram included) and 
into the continuum: T= Tp.. + Troon. The evaluation of Tp.. does not 
offer any difficulty, of course, and it can be done with m, + 0. In particular 
as m,—0, Tr... reproduces, for photoproduction, the Kroll-Ruderman 
theorem. For T,,,, one can try to use the current algebra constraints, for 
instance in the form of Eq. (3.4) of p. 130 but these are not enough to 
determine completely all invariant amplitudes +°. We can say in general 
that no explicit estimate of the corrections alternative to our approach 
has been, at our knowledge, performed. 


Appendix 
The Cluster Decomposition 


The cluster decomposition provides the systematic procedure to identify 
the various contributions to our completeness sums. The starting point 
is that every matrix element can be separated into a connected and a 
disconnected contribution. Indeed, if we consider an amplitude of the 


form M=<plAlpı,P2-:- Prs B> (A.1) 


where A is a certain operator, p,p,...p, are momenta of identical 
particles and ß represents all other particles appearing in the final state, 


we have M = Mom + aise l (A.2) 


where i 
M*** = <p| p> <O|A|p2... P> 


+ <plp2> <0lAlpı. DE 
and the + signs reflect the statistics obeyed by the particle labelled with p. 


(A.3) 


19 An exception seems to be the work by Scadron and Silbar, where in photoproduction 
not only M{+® but also M{~, Mz,4 are claimed to be determined by current algebra. 
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An explicit form for M°" and M““ can be obtained using the 
asymptotic creation and destruction operators. The relevant formula is 


<plAlny = COl [asu (P), Alln> 
+ COJA alpin? - 

The first term in the r.h.s. corresponds to the fully connected con- 
tribution, as is easily seen using the LSZ reduction technique. The second 
term represents the disconnected part, since the operator a(p) acts as a 
destruction operator on the state |n) = |p, ...p,, B>, leading automatically 


to Eq. (A.3)?°. In our applications to current algebra we consider only 
single particle external states: 


<P2|A|P1> = <p21Alpı 2°" + <P2|P1> £0410) 


<p2|A|pi>°™" = <O|[[a(p2), A], a* (pı)]10> 
and, moreover, A = A, A}. 
Separating out the fully disconnected part, and expanding the double 
commutator, we get: 


M321 =<P2|A2A1|Pr™" 
= Ol[a2, Az] [A1, af 110% + <0|42 [a2, [Ai a7 J] 10> (A.6) 
+ <O|[[a., 42], a; | A,l0y+<0|[A,,a; ][a,, 4,]]9>, 


where the + signs apply respectively to external bosons and fermions. 
At this point we use a completeness, and find: 


M,,=M3,+M3,+M?), (A.7) 
M31 = $ <Olla,, 42] |”) <n|[A;, 47110) 


(AA) 


(A.5) 


A.8 
= 2 €P2| A [ayena <n| A, PA ‘ecm ( ) 
M3,=% <OllLa;, 42], ar] |m> <m|A,|0) 
+ = (0) A, Im’) <m’|[a,, [A,, a7 ])|0> 
4 (A.9) 
=) <p2|A2|m, pX)" <m|A,|0> 
+ Y<OlAs|m’> <m, pa1A ID)", 
M3 = + %, COLAg, af WD <i Lap, A,110> 
(A.10) 


2 <O| Ag II, p> <1, p214,10> . 


20 Since |p> is a one particle state, it is inessential to specify whether it is ingoing 
or outgoing. 
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In order to understand the meaning of the different terms, let us 
suppose explicitly |p; >, |p2> to be nucleon states, and A, , pseudoscalar 
operators with zero baryonic number (for example, axial charges). In 
this case we see that: 


1. The states |n) have baryonic number +1 (the lowest one is the 
nucleon state), and MY, corresponds to direct graphs, namely to the 
transition from a nucleon state to an intermediate one and again to the 
final nucleon state. 


2. The states |m), \m’> have zero baryonic number, and M®, represents 
the “semidisconnected graphs” (or “mass singularities”), where, e.g. a 
state |m) is created from the vacuum and the reaction m+ p,—A,+D), 
follows. The lowest contribution to Eq. (A.9) will come from the single 
pion state. Under certain conditions, this will be the dominant contribu- 
tion to the completeness sum. The results of soft pion theory can thus be 
derived from axial charge-commutators, under the assumption of pion 
dominance. 


3. The states |!) have baryon number — 1 and M4! is given by the 
contributions of “Z-graphs”, corresponding to the creation from the 
vacuum of the state |l,p,), with subsequent annihilation of the state 
Il, pı> (the simplest effect is the creation and destruction of a nucleon- 
antinucleon pair). 
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